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ABSTRACT 

V 

Decisionmaking  organizations  are  modeled  as  asynchronous  concurrent 
systems,  using  Timed  Petri  Nets,  The  modeling  allows  for  evaluating  the 
time-related  performances,  with  respect  to  the  following  measures:  The 
maximum  throughput  rate,  defined  as  the  maximum  processing  rate  achievable 
by  the  system,  and  the  execution  schedule,  which  determines  the  earliest 
instants  at  which  the  different  operations  can  occur  in  the  process.  These 
measures  of  performance  are  analyzed  and  expressed  as  a  function  of  the 
resource  and  time  constraints  that  are  specific  to  the  organization.  The 
characterization  obtained  makes  it  possible  to  compare  different 
organizational  forms  and  to  modify  existing  designs  so  as  to  improve 
performance. 
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CHAPTER  I 


INTRODUCTION 


1.1  BACKGROUND 

In  most  decisionmaking  organizations  (DM0),  and  especially  in  military 
organizations  supported  by  C*  (command,  control,  and  communications) 
systems,  timeliness  is  of  critical  importance.  The  ability  of  an 
organization  to  perform  the  tasks  and  transmit  information  in  a  timely 
manner  is  indeed  a  determinant  factor  of  effectiveness.  There  are  two 
types  of  contraints  which  affect  the  time  performance  of  a  DM0.  The  first 
type  is  related  to  the  internal  organizational  structure  that  determines 
how  the  various  operations  do  occur  in  the  process:  some  tasks  are  indeed 
processed  one  after  the  other,  characteristic  of  sequential  activities,  and 
others  are  processed  independently,  which  characterizes  concurrent 
activities.  The  sequential  and  concurrent  events  are  precisely  determined 
by  the  communication  and  execution  procotols  among  the  individual 
organization  members,  as  defined  in  [1] .  The  second  type  of  constraints 
consists  of  time  and  resource  constraints.  More  precisely,  the  time 
constraints  derive  from  the  task  execution  times,  i.e.,  the  amount  of  time 
necessary  to  perform  each  task.  In  addition,  the  organization  controls 
specified  resources,  which  are  generally  limited.  Depending  on  which  of 
the  resources  are  free  at  a  given  Instant,  oertain  activities  oan  take 
place  and  others  must  be  delayed:  Hence,  the  resource  limitations  also 
constrain  the  processing  of  inputs. 

This  set  of  constraints  makes  the  decisionmaking  prooess  ooour 
asynchronously  and  concurrently  in  real-time.  Indeed,  the  coordination 
among  the  concurrent  tasks  occurs  asynchronously  and  not  at  set  times.  At 
this  point,  the  Petri  Net  formalism,  introduced  in  [2]  for  modeling  the 
DM0,  provides  a  convenient  tool  for  analyzing  the  behavior  of  systems  that 
exhibit  asynchronous  and  concurrent  properties.  The  Petri  Net  framework 


will  therefore  be  used  in  this  thesis  to  evaluate  the  performance  of  an 
organization,  with  respect  to  time-related  measures. 

1.2  GOALS 

In  earlier  work  [3],  the  evaluation  of  a  particular  measure  of 
performance  (HOP),  the  response  time  of  the  DMO,  has  been  carried  out, 
using  the  Petri  Net  representation.  This  MOP  characterizes  the  time 
interval  between  the  momemnt  one  external  input  or  task  is  received  by  the 
organization  and  the  moment  a  response  can  be  made.  The  response  time  so 
computed  is  a  static  measure,  in  so  far  as  the  processing  is  supposed  to 
take  place  for  a  single  input  arriving  at  any  instant  of  time. 

In  the  analysis  developed  here,  we  want  to  Investigate  the  dynamic 
behavior  of  the  DMO,  for  successively  arriving  inputs.  More  precisely,  the 
objective  of  this  research  task  is  to  evaluate  the  performance  of  a  DMO, 
with  respect  to  the  following  time-related  measures: 

(a)  Maximum  Throughput  Rate  of  the  system.  Clearly  enough,  if  we 
assume  that  external  inputs  are  arriving  continuously  at  a  rate 
which  is  low  enough,  the  DMO  will  be  able  to  handle  all  inputs  as 
soon  as  they  arrive.  Then,  the  rate  at  which  inputs  are  being 
processed,  which  characterizes  the  throughput  rate  of  the 
organization,  will  precisely  correspond  to  the  arrival  rate  of 
inputs.  However,  beyond  a  certain  arrival  rate  of  inputs,  the 
DMO  will  be  overloaded:  inputs  will  queue  at  the  entry  of  the 
system  and  this  queue  will  grow  to  infinity  over  time,  meaning 
that  the  DMO  will  never  catob  up  to  all  the  inputs.  This  bound 
precisely  determines  the  maximum  throughput  rate  of  the  system, 
i.e.,  the  maximum  rate  at  which  inputs  can  be  processed.  This 
MOP  characterizes  therefore  the  best  performance  that  can  be 
achieved  by  the  organization,  with  respect  to  the  processing 


(b)  The  Execution  Schedule  of  the  organization.  Let  us  assume  that 
the  processing  starts  at  x=0  and  that  it  occurs  repetitively  at 


its  maximum  throughput  rate. 


We  want  to  determine  the  earliest 


Instants  of  time  at  which  the  various  tasks  can  occur  in  the 
repetitive  process.  In  our  analysis,  the  DMO  will  be  allowed  to 
process  simultaneously  more  than  one  input:  we  will  therefore 
obtain  dynaml c  measures  of  performance.  In  particular,  if  we 
assume  that  several  inputs  arrive  simultaneously,  it  will  be 
possible  to  determine  the  response  time,  which,  in  that  case, 
corresponds  to  the  time  interval  between  the  moment  inputs  were 
received  and  the  moment  a  response  to  each  of  the  Inputs  was 
made. 

These  MOPs,  as  described  above,  are  important  in  the  performance 
evaluation  of  an  organization.  Suppose,  for  instance,  that  the  goal  of  a 
C*  system  is  to  detect,  track,  and  identify  threats  and  allocate  weapons  to 
them.  If  several  threats  arrive  simultaneously  (or  within  a  short  period 
of  time)  it  will  be  possible,  from  the  execution  schedule,  to  evaluate  the 
ability  of  the  system  to  respond  to  all  these  threats  in  a  timely  manner. 
Likewise,  if  a  large  number  of  threats  arrive  successively  with  a  rate 
that  is  above  the  maximum  throughput  rate  of  the  system,  we  know  that  the 
system  will  be  overloaded. 

1 .3  APPROACH 

The  performance  evaluation  is  based  on  the  Petri  Net  modeling  of 
decisionmaking  organizations,  as  introduced  in  [2].  However,  some 
extensions  are  introduced,  in  order  to  account  for  the  time  and  resource 
constraints.  The  resource  constraints  actually  occur  for  both  individual 
decisionmakers  (DMs)  and  the  organization  as  a  whole.  In  our  approach,  the 
resource  limitation  existing  for  individual  DMs  allows  for  modeling  the 
limited  capacity  of  information-processing  that  characterizes  a  human  DM, 
as  described  in  the  information  theoretic  framework  by  the  bounded 
rationality  constraint  [4],  More  precisely,  according  to  the  analysis 
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carried  out  in  IS] ,  the  resource  constraint  is  derived  from  the  limited 
’ability  of  the  human  short-term  memory  to  handle  a  certain  amount  of 
information  at  the  same  time.  Therefore,  the  limitation  of  resources  is 
expressed  here  as  a  bound  on  the  total  number  of  inputs  that  a  DM  is  able 
to  process  simultaneously.  Likewise,  the  amount  of  information  being 
transmitted  between  the  DMs  is  necessarily  limited  by  the  storage  capacity 
of  the  organization.  We  should  recall  that  the  process  ocours 
asynchronously,  which  implies,  generally,  that  the  information  exchanged 
between  the  DMs  has  to  be  stored  temporarily.  Exactly  as  for  individual 
DMs,  this  capacity  constraint  is  modeled  as  a  bound  on  the  total  number  of 
Inputs  that  the  overall  organization  can  handle  at  the  same  time. 

As  already  stated,  the  time  constraints  correspond  to  the  task 
execution  times,  i.e.,  the  amount  of  time  it  takes  to  process  each  task. 
Consequently,  arbitrary  processing  times  are  assigned  to  each 
organizational  task,  according  to  the  Timed  Petri  Net  model  developed  by 
Ramchandani  [6].  The  performance  evaluation  of  the  DMO  is  therefore 
carried  out  using  the  properties  of  Timed  Petri  Nets. 

In  this  work,  the  maximum  throughput  rate  is  expressed  directly  as  a 
function  of  both  the  task  processing  times  and  the  resources  available, 
given  the  structure  of  the  organization.  In  addition,  an  algorithm  ia 
developed,  that  determines  the  execution  schedule  of  the  organization,  when 
the  processing  of  inputs  occurs  at  its  maximum  throughput  rate.  These  MOPs 
characterize  in  fact  the  best  dynamic  performance  achievable  by  the 
organization,  as  determined  by  the  maximum  rate  at  which  inputs  can  be 
processed  and  the  earliest  instants  of  time  at  which  the  various  tasks  oan 
be  exeouted.  These  measures  provide,  therefore,  a  useful  tool  to  evaluate 
and  compare  different  organizational  forms.  In  particular,  the  precise 
characterization  of  the  resource  and  time  constraints  makes  it  possible  to 
investigate  which  of  the  constraints  should  be  modified  in  the  aotual 
design,  so  as  to  improve  the  time  performance  of  the  organization. 


1.4  THE  THESIS  IN  OUTLINE 


The  thesis  is  organized  as  follows.  Chapter  2  is  a  review  of  Petri 
Net  theory  and  presents  the  extended  Timed  Petri  Nets  that  are  used  in  the 
analysis.  In  Chapter  3.  the  model  of  the  DMO  using  Timed  Petri  Nets  is 
discussed.  Included  in  the  modeling  are  the  resource  and  time  constraints, 
as  described  previously.  In  Chapter  4,  the  maximum  throughput  rate  is 
determined  as  a  function  of  the  time  and  resource  constraints.  The 
analysis  of  the  execution  schedule  is  carried  out  in  Chapter  3  and  a  method 
is  developed  to  obtain  the  precise  instants  at  which  the  various  operations 
take  place  in  the  process.  In  Chapter  6,  an  application  is  developed  for  a 
five  member  organization  and  the  corresponding  measures  of  performance  are 
computed.  Finally,  conclusions,  as  well  as  suggestions  for  further 
research  are  presented  in  Chapter  7 . 
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CHAPTER  II 


REVIEW  OF  PETRI  NET  THEORY 


Petri  Net  theory  forms  the  framework  for  the  model  of  the 
decisionmaking  organization  analyzed  in  this  thesis.  Accordingly,  the 
basic  definitions  are  first  reviewed  in  this  chapter,  as  they  are  developed 
in  [7].  The  important  concepts  of  S  and  T-invariants  [8]  are  also 
presented  and  will  be  used  later  to  develop  algorithms.  In  addition,  the 
model  of  Timed  Petri  Nets  [6]  is  described,  since  it  is  the  approach  taken 
in  this  thesis  to  carry  out  the  performance  evaluation  of  the  DMO. 

2.1  BASIC  DEFINITIONS 

2.1.1.  Structure  of  a  Petri  Net 


A  Petri  Net  (denoted  by  PN)  is  a  quadruple 
PN  =  (P,  T.  I,  0) 


where : 

P  =  (p1,...,pn)  is  a  finite  set  of  places 
T  =  is  a  finite  set  of  transitions 

I  is  a  mapping:  P  x  T - •>  (0,1}  corresponding  to  the  set  of  directed 

edges  from  places  to  transitions 

0  is  a  mapping:  T  x  P - >  {0,1}  corresponding  to  the  set  of  direoted 

edges  from  transitions  to  places. 

Figure  2.1  shows  a  simple  example  of  PN. 


Figure  2.1  Example  of  a  Petri  Net 


We  denote  by: 


t*  ■  {p  e  P  /  0(t,p)  =  1}  the  set  of  all  output  places  of  transition  t. 
't  =  {p  c  P  /  I(p,t)  =  1}  the  set  of  all  Input  places  of  transition  t. 
p*  and  'p  are  defined  for  each  place  p  in  a  similar  way. 

In  the  example  above:  't%  =  (p^p,)*  =  (pa.p4). 

2.1.2  Incidence  Matrix 

The  Incidence  Matrix  (also  called  Flow  Matrix),  denoted  by  C, 
characterizes  the  structure  of  a  PN  In  the  following  way:  The  columns  of 
the  matrix  correspond  to  the  transitions  of  the  net  and  the  rows,  to  the 
places.  C  Is  therefore,  a  n  x  m  matrix,  such  that: 


-1  if  Pj^  is  an  Input  place  of  tj 
+1  if  p^  is  an  output  place  of  tj 
0  if  there  is  no  arc  between  pA  and  tj 


We  can  easily  verify  that: 


C, ,  =  0(t„  p,)  -  I(p,,t  J 


(2.1) 


The  Petri  Net  of  Figure  2.1  has.  for  example,  the  following  inoidence 
matrix: 


2.1.3  Marked  Petri  Net 

We  call  marking  of  a  PN,  denoted  by  M,  a  mapping:  P - >  10,1,2,...,} 

(set  of  integers)  which  assigns  a  non-negative  integer  number  of  so  called 
tokens  to  each  place  of  the  net. 

Figure  2.2  shows  an  example  of  marking  for  the  PN  considered  in 
Section  2.1.1.  Dots  represent  tokens.  The  marking  is  usually  denoted  as  a 
n-dimensional  integer  vector,  with  each  element  corresponding  to  one  of  the 
places.  For  the  example  of  Figure  2.2: 

M  =  (2,  0.  1,  0)T 


Figure  2.2  Example  of  a  Marked  Petri  Net 
2.1.4  Rules  of  Operation 

(i)  A  transition  t  is  enabled  by  a  marking  M#  if  every  input  place 


M1^)  =  M*(pi)  +  Cjj 


for  j  and  all  1*1,2 


n 


or,  in  a  more  compact  form: 

M1  =  M®  +  C  Xj  (2.3) 

where  Xj  is  a  m-dimensional  vector  with  all  components  equal  to  zero  except 
the  j-th  one,  which  equals  1: 

Xj  =  [0,0 . 1 . 0]T 


We  can  now  consider  a  sequence  of  transition  firings,  denoted  by: 


a 

s 


which  means  that  transition  tj4  fires  first,  then  transition  tja  and  so  on 
until  transition  tja,  We  denote  by 

the  marking  reached  from  the  initial  marking  M*  by  firing  the  sequence  a&. 
Let  us  construct  the  m-dlmensional  integer  vector: 


V  [v  v 


-n  1 
m 


where  nj  ( j*l ,2, . . . ,m)  denotes  the  number  of  oocurences  of  transition  tj  in 
the  sequence  og.  Then  relation  (2.3)  is  generalized  in  the  form  [9]: 


Ms  =  M°  +  C  N 


(2.4) 
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The  algebraic  equation  (2.4)  allows  for  computing  directly  the  new  marking 
reached  by  any  sequence  of  transition  firings.  Let  us  remark  however,  that 
some  information  is  lost  when  using  the  vector  Na,  since  the  order  of  the 
firing  sequence  is  not  specified. 

2.1.5  Forward  Harking  Class 

Given  an  initial  marking  M#  of  the  net,  we  call  forward  marking  class, 
denoted  by  M°,  the  set  of  all  possible  reachable  markings.  In  other  words, 
a  marking  M  belongs  to  M°  if  there  exists  a  firing  sequence  o,  such  that: 


In  the  example  shown  in  Figure  2.2,  the  forward  marking  class  consists  of 
the  following  markings: 


(2,  0,  1.  0)T 

(1,  1,  1,  0)T 

(0,  0.  2.  1)T 

(0,  2,  1,  0)T 


(Initial  marking) 

(after  firing  of  tx) 

(after  firing  of  tx) 

(after  two  successive  firings  of  tx) 


2.1.6  Liveness  and  Boundeness 

Liveness :  A  marking  M8  is  live  if,  for  any  transition  t,  and  for 
every  reachable  marking  M  (i.e.,  M  s  M#),  there  exists  a  firing  sequence 
from  M  which  fires  t.  In  other  words,  this  property  guarantees  the  firing 
process  to  be  deadlock  free. 

The  example  of  Figure  2.2  is  trivially  not  live:  from  M*,  it  is  not 
possible  to  fire  any  more  transitions.  Figure  2.3  shows  a  simple  example 
of  a  live  Petri  Net. 


V v  v .r- .%  ,v\vlv*vv  va’//  \*  v  -  ”  w  v  o  o  ov\’ 

■  / •  O  ■*.*  O  s  '  ***»  •'«  *,fc*^s*\,**<* * 


Figure  2.3  Example  of  a  Live  Petri  Net 

Boundedness :  A  marking  M°  is  bounded,  if  there  exists  a  positive 
integer  N  such  that,  for  every  reachable  marking  M,  the  number  of  tokens  in 
each  place  is  bounded  by  N.  If  N  equals  one,  the  marking  is  said  to  be 
safe. 

The  example  of  Figure  2.2  is  trivially  bounded  (by  2).  However,  the 
example  above  (Figure  2.3)  is  not  bounded.  Each  time  that  tx  fires,  one 
token  is  added  in  pB.  Since  t4  can  fire  an  infinite  number  of  times,  the 
marking  of  p,  is  not  bounded. 

Remark:  There  exists  an  extension  to  the  definition  of  liveness  (resp. 
boundedness),  which  is  called  structural  liveness  (resp.  structural 
boundedness)  [10].  In  that  case  the  properties  are  characteristics  of  the 
PN,  independently  of  the  initial  marking.  Here  sure  the  definitions: 

(1)  A  PN  is  (structurally)  live  if  there  exists  at  least  one  initial 
marking  which  is  live. 

(2)  A  PN  is  (structurally)  bounded  if  any  Initial  marking  is  bounded. 
2.1.7  Application  of  the  Petri  Net  Model 


In  systems,  especially  those  which  exhibit  asynchronous  and  concurrent 
properties,  as  it  is  the  case  for  the  decision  making  process  studied  in 
this  thesis.  To  model  the  dynamic  behavior  of  a  system,  the  execution  of  a 
process  (or  task)  is  represented  by  the  firing  of  the  corresponding 
transition.  The  flow  of  tokens  represents  the  flow  of  information  in  the 
process  and  the  marking  of  the  net  at  any  instant  corresponds  to  a 
particular  state  of  the  system.  The  forward  marking  class  determines  the 
set  of  all  possible  states,  given  an  initial  state  of  the  system.  The  two 
properties  defined  above,  namely  liveness  and  boundedness,  characterize  a 
well  defined  system.  Indeed,  they  guarantee  that  first,  the  process  is 
deadlock  free  (liveness)  and  second,  that  information  does  not  accumulate 
at  some  steps  of  the  process  (boundedness).  It  becomes  possible, 
therefore,  to  analyze  the  steady-state  of  the  process. 

In  the  next  sections  of  this  Chapter,  we  are  going  to  present  some 
extensions  to  the  basic  Petri  Net  model,  that  will  be  used  later  in  our 
development. 

2.2  S  AND  T- INVARIANTS 

S  and  T-invariants  are  important  concepts  in  Petri  Net  theory  [8]  and 
will  be  used  later  to  describe  algorithms.  This  is  the  reason  why  they  are 
presented  in  this  Chapter.  S  and  T-invariants  are  dual  in  the  following 
sense:  Let  us  consider  a  Petri  Net  in  which  we  replace  all  the  places  by 
transitions  and  conversely.  The  Petri  Net  obtained  is  oalled  the  dual  [11] 
of  the  original  Petri  Net  and  it  is  easily  proved  that  the  incidence  matrix 
of  the  new  PN  is  the  tranpose  of  the  incidence  matrix  of  the  original  PN. 
We  will  show  later  that  the  S-invariants  of  a  PN  are  the  T-invariants  of 
the  dual  and  conversely.  This  is  the  reason  why  we  will  focus  in  this 
section  on  S-invariants  only.  The  properties  of  T-invariants  will  be 
obtained  directly  from  the  analysis  of  S-invariants . 


2.2.1  Definitions 


We  reeall  that  C  represents  the  incidence  matrix  of  a  PN,  as  defined 
in  Section  2.1.2.  C  is  a  n  x  m  matrix,  where  n  is  the  number  of  plaoes  and 
m  the  number  of  transitions. 

Definition:  An  n-dimensional  positive  integer  vector  X  is  called  an  S- 

invariant  if  and  only  if  XT  C  =  0. 

Definition:  The  set  of  places  whose  corresponding  components  in  X  are 

strictly  positive  is  called  the  support  of  X  and  is  denoted  by  1 1 X I  I . 

Definition:  The  support  1 1 X 1 1  of  an  S-invariant,  X,  is  3aid  to  be  minimal 
if  and  only  if  it  does  not  contain  the  support  of  any  other  invariant,  but 
Itself  and  the  empty  set. 

2.2.2  Example  of  S-Invariants 

Let  us  consider  again  the  example  of  Figure  2.1.  We  recall  that  the 
Incidence  Matrix  is: 


An  S-invariant:  X  =  [xx,  xa,  x, ,  x4]T  (where  xx,  x#,  x,  and  x4  are  positive 
or  null  integer)  satifies: 

T 

X  C  =  0 

which  yields  the  following  system  of  equations: 
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Therefore,  the  general  form  of  an  S-lnvariant  is: 

X  =  lx,  x,  y,  zlT 

such  that: 

y  +  z  =  2  x 

In  particular,  it  turns  out  that: 

Xx  -  [1  1  2  0]T 

X  -  [1  1  1  1]T 
2 

X,  -  [1  1  0  2]T 

are  all  S-invariants,  whose  supports  are  respectively: 

llxJI  *  P,*  P,J 

I  lX, 1 1  =  tP*'  P4'  P,'  P4} 

IIMI  -  tPa'  P,'  P4} 

Now,  the  support  of  any  non-null  S-invariant,  X  «=  [x,  x,  y,  z] ,  oontalns 
necessarily  px  and  ps.  Otherwise  we  would  have  x«0  and  therefore  y  ■  *  *  0 
(since  y  and  z  are  positive  integers  verifying:  y  +  z  ■  2  x).  In  addition, 
y  and  z  cannot  be  both  null  simultaneously,  which  implies  that  the  support 
of  any  non-null  S-invarlant  also  oontalns  at  least  p,  or  p4.  Accordingly, 


IKII  and  llx.ll  are  minimal  support  of  S-lnvariants.  By  construction, 
both  are  in  fact  the  only  minimal  supports,  since  we  have  just  proved  that 
the  support  of  any  non-null  S- invariant  neoessarlly  contains  (px,  p#,  pt) 
(i.e.,  |  |Xx  1 1 )  or  lp4,  pa,  p4)  (i.e.,  ||X,||). 


2.2.3  Properties  of  S-invariants 


The  fundamental  property  of  an  S-invariant,  which  Justifies  "a 
posteriori*  its  name,  is  the  following  one: 


Theorem  2.1:  X  is  an  S-invariant  iff  for  any  initial  marking  M°  and  for 
any  reachable  marking  M  (i.e.,  M  e  M  ), 


T  To 

X  M  -  X1  M 


(2.5) 


The  proof  is  straightforward  from  Equation  (2.4)  and  from  the  definition  of 


am  S-invariant,  i.e.,  X1  C  =  0. 


Relation  (2.5)  is  very  similar  to  an  equation  of  convervation:  it 
implies  indeed  the  conservation  of  tokens  for  the  places  belonging  to  the 
support  II X 1 1  and  weighted  by  certain  coefficients. 


For  instance,  in  the  example  of  Section  2.2.2,  the  three  S-invariants 
yield  the  following  equations  of  conservation: 


M(px)  +  M(pa)  +  2M(p|)  =  M  (p4)  +  M  (p#)  +  2M  (p$) 


(for  X,  -  [1  1  2  0]1) 


M(px)  +  M(p#)  +  2  M(p4)  =  M  (pj  +  M  (pa)  +  2  M  (p4) 


(for  Za  -  [1  1  0  21' 


}  M(Pl)  =  J  M#  (p± ) 


(for  X,  «  [1  1  1  l]1) 


I 

I 

Ki’ 


This  last  equation  is  especially  easy  to  interpret:  there  is  neither 
creation  nor  loss  of  tokens  in  the  net. 

The  next  theorem  is  the  basis  of  algorithmic  procedures  to  determine 
all  the  S-invariants  of  a  Petri  Net: 

Theorem  2.2:  (Memmi-Sifakis)  [8]:  Let  I#  be  the  support  of  an  invariant 
and  lx,....Ik  be  the  minimal  supports  contained  in  I,.  then: 


-  for  every  invariant  X  such  that  1 1 X 1 1  =  I#,  there  exists  positive 
rational  coefficients  (1-1. k)  such  that: 

k 

x-5  IjXj 
1=1 


where  XA  is  an  invariant  whose  support  is  Ii. 

Corollary  2.2:  Let  IX,...,IS  be  all  the  miniaal  supports  of  S-invariants 
of  a  PN  (they  are  necessarily  finite,  sinoe  the  number  of  plaoes  is 
finite).  Let  X^^  be  an  invariant  whose  support  is  IA,  i.e.,  11x^1  ■  (we 
can  choose  actually  the  minimal  one.  as  determined  by  the  partial  ordering 
2  between  integer  vectors).  Then  the  set  ^.....Xg)  forms  a  minimal  set 
of  generators  of  all  invariants,  i.e., 

-for  every  invariant  X,  there  exists  positive  rational  coefficients 
such  that 

s 

x '} 

i=l 
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-if  one  element  is  removed  from  the  set  -  let  us  say  X^  for  example  - 
the  set  {Xl,...Xs}  -  {Xk}  is  not  generator. 

In  the  example  of  Section  2.2.2,  Ix  *  {px,  p,.  p, }  and  I,  -  {px,  p4,  p4) 
are  the  minimal  supports.  Xx  =  [1  1  2  0]T  and  X,  ■  [1  1  0  2]T  form  the 
minimal  set  of  generators  and  we  can  check  for  example  that: 

Xa  =  [1  1  1  1]T  =  i  X4  +  ^  X, 

i  •  e  • ,  Xj  2 

2.2.4  T-Invariants 

As  stated  earlier,  T- invariants  can  be  considered  as  duals  of  S- 
invariants.  A  T-invariant  is  therefore  a  m-dimensional  positive  integer 
vector  Y  such  that 

C  Y  -  0  (2.6) 

The  set  of  transitions  whose  corresponding  components  in  Y  are  strictly 
positive  is  called  the  support  of  Y  and  is  denoted  by  INI-  The 
definition  of  the  minimal  support  is  the  same  as  for  ^invariants,  but  it 
refers  now  to  a  set  of  transitions.  Likewise,  Theorem  2.2  holds  also  for 
T-invariants. 

From  equation  (2.6),  a  T-invariant  can  be  Interpreted  as  a  positive  or 
null  integer  assignment  to  each  transition  of  the  net  suah  that,  at  every 
place,  the  sum  of  integers  assigned  to  its  input  transitions,  equals  the 
sum  of  integers  assigned  to  its  output  transitions. 

2.3  PETRI  NET  THEORY  AND  GRAPH  THEORY 

Formally,  a  Petri  Net  can  be  defined  as  a  bipartite  directed  graph. 


whose  nodes  are  places  (also  called  oircle  nodes)  and  transitions  (also 
called  bar  nodes)  and  whose  arcs  correspond  to  the  directed  edges  between 
places  and  transitions.  We  give  first  some  definitions  that  come  directly 
from  Graph  Theory. 

2.3.1  Definition  from  Graph  Theory 

Connectivity:  A  Petri  Net  is  said  to  be  connected  if  and  only  if  there 
exists  a  path  (not  necessarily  directed)  from  any  node  (place  or 
transition)  to  any  other  node . 

Strong  Connectivity:  A  Petri  Net  is  said  to  be  strongly  connected  if  and 
only  if  there  exists  a  directed  path  from  any  node  to  any  other  node. 

Simple  connectivity  implies  only  that  a  PN  cannot  be  split  in  two  or 
more  Petri  Nets  that  are  independent.  All  the  examples  of  PNs  presented 
previously  are  simply  connected.  Strong  connectivity  is  however  a  much 
more  stringent  characteristic.  The  example  of  Figure  2.1  is  obviously  not 
strongly  connected:  there  is,  for  example,  no  directed  path  between  p,  and 
p4.  Figure  2.4  shows  an  example  of  a  strongly  connected  net. 


Figure  2.4  Example  of  a  Strongly  Connected  let 

Directed  Elementary  Circuit:  A  directed  elementary  cirouit  is  a  directed 
path  from  one  node  back  to  Itself  such  that  none  of  the  nodes  are  repeated. 
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Naturally,  the  condition  that  none  of  the  nodes  should  be  repeated 
refers  to  the  fact  that  the  circuit  is  elementary.  A  directed  path  from  a 
node  baok  to  Itself  is,  generally  speaking,  a  directed  circuit  (not 
necessarily  elementary).  However,  we  will  always  oonaider,  in  the  sequel, 
only  the  elementary  circuits  of  a  PN.  Tor  simplification,  what  will  be 
called  later  directed  circuits,  shall  always  refer  to  elementary  circuits. 

In  the  example  of  Fig.  2.4,  the  sequence  px  =  (px  tx  px  tx  p,  t,) 
constitutes  a  directed  elementary  circuit.  However,  the  sequence 
p  =  (px  tx  px  tx  p,  t,  px  tx  p,  t,)  is  a  directed  circuit,  which  is  not 
elementary. 

2.3.2  Event-Graph:  Definition  and  Properties 

Event-Graphs  constitute  a  special  class  of  Petri  Nets,  whose 
properties  are  especially  appropriate  for  analyzing  the  dynamic  behavior  of 
a  system.  In  this  section,  we  present  the  relevant  definitions  and 
properties  that  will  be  used  later  in  this  thesis. 

Definition:  An  Event-Graph  [9]  (also  known  as  Marked  Graph  [12])  is  a 

connected  Petri  Net,  in  which  each  place  has  exactly  one  input  and  one 
output  transition. 


This  means  that  tokens  are  always  generated,  at  a  given  place,  by  a 
predefined  transition  (its  only  input  transition)  and  oonsumed  by  a 
predefined  transition  (its  only  output  transition).  Neither  the  example  of 
Figure  2.1,  nor  the  example  of  Figure  2.4  are  Event-Graphs.  Indeed,  plaoe 
px  has,  in  both  oases,  two  output  transitions.  Figure  2.5  is  an  example  of 
an  Event-Graph.  We  should  point  out  that,  in  this  particular  example,  the 
input  transition  (i.e.,  transition  tx)  has  no  input  places.  This  is 
allowable  in  Petri  Net  theory  and  simply  means  that  the  corresponding 
transition  is  always  enabled.  Otherwise  stated,  transition  tx  oan  fire  any 
number  of  times  at  any  instant. 


Figure  2.5  Example  of  an  Event-Graph 


The  following  two  theorems,  due  to  Commoner  and  Holt  [12],  describe 
useful  properties  of  Event-Graphs: 

Theorem  2.3:  In  an  Event-Graph,  the  number  of  tokens  in  any  directed 
elementary  cirouit  (called  the  token  content)  remains  invariant  by 
transition  firings. 

Theorem  2.4:  A  marking  of  an  Event-Graph  is  live  if  and  only  if  the  token 
content  of  every  directed  elementary  circuit  is  strictly  positive. 

For  instance,  the  marking  of  the  event-graph  shown  in  Figure  2.6  is 
live:  The  token  content  of  the  unique  oircuit  (ts  p,  t,  p,  t4  p4)  is  one. 


P* 


Figure  2.6  Example  of  a  Live  Event-Graph 

The  other  interesting  result  relates  the  boundedness  property  to  the 
condition  of  strong  connectivity  in  the  following  way: 


Tbeorem  2.5:  An  Event-Graph  is  bounded  if  and  only  if  it  is  strongly 
connected. 

In  order  to  prove  this  theorem,  we  use  a  result  obtained  by  Sifakis 
[10],  stating  that  a  Petri  net,  which  is  (structurally)  live  and  bounded, 
is  necessarily  strongly  connected.  Let  us  now  prove  Theorem  2.5. 

(i)  Let  us  assume  that  the  Event-Graph  is  bounded .  Prom  Tbeorem  2.4, 

we  know  that  an  Event-Graph  is  always  (structurally)  live: 
Indeed,  it  is  sufficient  to  choose  an  initial  marking,  such  that 
the  token  content  of  every  directed  elementary  circuit  is 
strictly  positive.  Given  the  result  stated  above,  the  Event- 

Graph  is  therefore  strongly  connected,  since  it  is  live  and 
bounded . 

(ii)  Let  us  assume  now  that  the  Event-Graph  is  strongly  connected  and 

let  p±  (resp.  tj)  be  any  place  (resp.  any  transition).  There 
exists  a  directed  path  from  pA  to  tj,  that  we  denote  by  and 

a  directed  path  from  tj  to  pA,  that  we  denote  Now,  the 

sequence  of  nodes  belonging  to  the  union  of  the  two  paths, 
wij+<Tji'  determines  a  circuit  in  the  net.  The  token  content  is 
invariant  in  the  circuit  (Theorem  2.3)  and  the  marking  of  pi  is 
therefore  bounded  by  this  token  oontent.  Hence  the  Event-Graph 
is  bounded.  Q.E.D. 

We  can  check,  for  instance,  that  the  PN  of  Figure  2.6  is  not  bounded: 
tx  can  fire  any  number  of  times,  as  said  previously,  so  that  the  marking  of 
pt  is  never  bounded.  Indeed,  the  PN  is  not  strongly  oonnected. 

It  turns  out  that  Theorem  2.4  and  2.5  provide  a  very  easy  way  to  oheck 
the  two  important  properties  -  liveness  and  boundedness  -  that  are 
characteristics  of  a  well-defined  system,  if  the  system  is  modeled  as  an 
Event-Graph . 


»j<l 


2.3.3  S-Invariants  of  Event-Gn_, — 


Before  characterizing  the  S-invariants  of  an  Event-Graph,  we  need  a 
further  definition:  S^oomponents  [10] . 

Definition:  Let  X  be  any  S- invariant  of  a  PN  and  llxll  its  support.  We 
recall  that  llxll  is  a  set  of  places.  We  call  S-ooaponent  the  unique 
subnet  whose  set  of  plaoes  is  precisely  llxll  and  whose  set  of  transitions 
consists  of  all  the  transitions  connected  to  the  plaoes  of  1 1 X 1 1 .  In  other 
words,  the  corresponding  S-component  is  the  subnet  PNa  *  (P8,  Ts,  Ig,  0S) 
such  that: 


P  =11x11 

S 

T  =  U  (P’U'P) 

pepa 

(l.e.,  T„  contains  all  the  input  and  output  transitions  of  the  plaoes 

3 

of  Ps).  Ia  (res.  0g)  is  the  restriction  of  I  (resp.  0)  to  Pg  u  Tg. 

Let  us  take  for  example  the  PN  shown  in  Figure  2.1.  The  corresponding 
S-component  of  the  S-invariant  Xt  =  [1  1  2  0]T  is  the  subnet  shown  in 
Figure  2.7. 

Given  the  definition  of  a  minimal  support,  an  S-component  is  said  to 
be  minimal  if  it  corresponds  to  an  S-invarl&nt  whose  support  is  minimal, 
which  means  also  that  it  does  not  contain  any  other  S-component  but  Itself 
and  the  empty  set. 


Figure  2.7  An  Example  of  an  S-oomponent 
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The  following  new  result,  shows  the  relationship  between  the  directed 
circuits  and  the  S-invariants  of  an  Event  Graph. 

Theorea  2.6:  The  minimal  S-oooponents  of  an  Event-Graph  are  exactly  the 
directed  circuits. 

Proof  of  Theorea  2.6 

(i)  Let  us  prove  first  that  a  directed  circuit  is  a  minimal  S-oomponent. 
Assume  then  that  p  is  any  directed  circuit  and  let  us  call  Pp  the  set 
of  places  of  the  circuit,  Tp  the  set  of  transitions  of  the  circuit. 
We  have  trivially: 

T  *’  U  (p'U'P)  (2.7) 

p  peP 

P 

because  each  place  has  exactly  one  input  and  one  output  transition. 
From  Theorem  2.3,  we  know  that,  for  any  Initial  marking  M*  of  the  net 
and  for  any  reachable  marking  M: 

5  M(p>  =  y  M°(p)  (2*8) 

peP  peP 

P  P 

Let  Xp  be  the  n-dimenslonal  vector  whose  components  corresponding  to 
the  places  of  Pp  are  equal  to  1  and  to  other  places  equal  to  0. 
Equation  (2.8)  can  be  written: 

XT  M  -  XT  M*  (2.9) 

P  P 

From  Theorem  2.1,  this  means  that  Xp  is  an  S-invariant.  Because  of 
relation  (2.7),  the  corresponding  S-component  is  the  same  as  the 
directed  circuit  p.  Let  us  prove  now  that  the  support  of  Xp,  i.e.. 


Pp,  is  necessarily  minimal.  If  we  oonsider  the  circuit  p  as  a  subnet 
(which  is  allowable  since  it  is  an  S-component) ,  we  can  construct 
its  Incidence  Matrix  Cp,  which,  after  an  appropriate  ordering 
of  plaoes  and  transitions,  has  the  form  (the  columns  correspond 
to  the  transition,  tp  ,  tp2,...,tn  ,  and  the  rows  to  the  plaoes, 
Po  »  Po  » • • *Po  ) • 


.  0  +1 


+1  -1 


+1  -1 
0  +1 


0  0  0 


+1  -1 


where  p  corresponds  to  the  circuit: 


P  =  (p  t  p  .  .  .  .  p  t  ) 

Pi  Px  P»  Pk  Pk 


Let  Xk  =  [ax  .  .  .  ak]T  be  any  positive  integer  vector.  Now  Xk  is 
^invariant  of  the  subnet  iff: 


X*  C  =  0 

k  p 


which  yields  immediately: 


Otherwise  stated,  all  the  S-invariants  of  the  subnet  have  the  same 
support,  which  is  actually  the  set  of  plaoes  of  the  oirouit,  i.e,  Pp. 
Hence  Pp  is  a  minimal  support  of  the  Event-Graph:  If  it  were  not,  it 
would  strictly  contain  another  minimal  support,  which  would  also  be  a 
support  of  invariant  for  the  corresponding  subnet.  However,  we  have 


•  •  .  •  S  .  J*  f  J*  «*■  T  ’  J*  *'  ,*  *  ■  a  .*  *  -  *.  *  w  m  "  V*  .  •  *  .  “  .  «  .  '  „  -  N  **  A  ^  . 

V - s  •.  •.  s  v  .  N  ••  V  •-  vA  A  - 


showed  that  Pp  cannot  strictly  contain  another  support  of  invariant. 


Q.E.D. 


(ii)  Let  us  prove  now  that  a  minimal  S-ooaponent  is  a  directed  circuit. 
Assume  then  that  PN_  is  a  minimal  S-component  corresponding  to  the  Sr- 
invariant  Xg  and  call  Ps  (resp.  T„)  its  set  of  places  (resp.  its  set 


S  *8' 

of  transitions).  By  construction,  PNS  is  an  Event-Graph.  For  any 

0 


initial  marking  M  and  for  any  reachable  marking  M,  we  deduce 
immediately  from  Theorem  2.1  that,  for  all  the  places  pi  belonging  to 


M  (p.)  < 


T  • 
X  M 
s 


where  xA  denotes  the  i-th  component  of  vector  Xs.  Hence,  PNS  is  a 
bounded  net.  From  Theorem  2.5,  we  deduce  that  PNg  is  a  strongly 
connected  Event-Graph.  We  can  find  therefore  at  least  one  directed 
circuit,  p,  in  PNS:  we  choose  any  two  nodes  and  consider  the  directed 
path  that  goes  from  the  first  node  to  the  second  one  and  vice-versa. 
However,  we  know  from  (i)  that  p  is  in  that  case  a  minimal  S-component 


of  PNg  and  consequently  of  the  original  Event-Graph  PN.  Since  PNS  is 


by  definition  a  minimal  S-component,  and  since  p  is  lnoluded  in  PN_, 

9 


PNS  is  necessarily  identical  to  p:  Hence,  PNg  is  a  directed  circuit. 

Q.E.D. 


Theorem  2.6  will  be  used  to  determine  all  the  circuits  of  an  Event- 
Graph,  from  an  algorithm  that  determines  all  the  minimal  support  S- 
invariants  of  a  net  [13], 


The  last  section  of  this  chapter  concerns  the  Timed  Petri  Net  model, 
as  it  will  be  used  to  represent  Decision  Making  Organizations.  Indeed, 
Timed  Petri  Nets  provide  an  appropriate  model  for  describing  the  dynamic 
behavior  of  systems,  especially  those  which  exhibit  asynchronous  and 
concurrent  properties. 
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2.4  TIMED  PETRI  NETS 


The  strength  of  the  Petri  Net  formalism,  as  described  In  the  previous 
sections,  lies  In  the  rather  siaple  representation  It  provides  for  modeling 
complex  systems,  in  which  processes  are  both  concurrent  and  sequential. 
However,  the  major  weakness  of  ordinary  Petri  Net  models  comes  from  the 
lack  of  tools  to  analyze  the  real-time  process.  Indeed,  no  assumptions  are 
made  regarding  the  length  of  time  it  takes  to  complete  the  different 
processing  operations.  Tokens  move  in  the  system  according  to  the 
transition  firings,  which  are  assumed  instantaneous. 

To  handle  this  shortcoming,  several  types  of  Timed  Petri  Nets  -  in 
which  the  notion  of  time  is  explicitly  introduced  -  have  been  developed  in 
the  literature.  The  main  differences  between  the  models  come  actually  from 
the  type  of  processing  times  considered:  deterministic  or  random. 
Ramchandanl  [6]  was  the  first  to  introduce  a  model  of  Timed  Petri  Net  in 
which  the  processing  times  were  assumed  deterministic  and  assigned  to  the 
transition  of  the  net.  Sifakis  [14]  also  considered  deterministic 
processing  times,  but  assigned  to  the  places.  He  has  proved,  however,  that 
both  models  are  equivalent.  Stochastic  Timed  Petri  Nets,  corresponding  to 
the  case  where  times  are  assumed  random  have  been  extensively  analyzed  by 
Wiley  [15].  In  this  thesis,  the  type  of  TPN  (Timed  Petri  Net)  introduced 
by  Ramchandanl,  will  be  used  for  modeling  the  DMO.  In  this  section,  the 
related  definitions  are  presented. 

2.4.1  Definitions 

A  Timed  Petri  Net  (TPN)  is  a  pair  (PN,  p),  where  PN  is  a  Petri  Net  and 
p  is  a  firing  time  function  that  assigns  a  positive  rational  number  to  each 

transition  of  the  net:  p  :  T  - >  R+  (R+  denotes  the  set  of  positive 

rational  numbers).  In  a  Timed  Petri  Net,  each  transition  t  takes  the  time 
p(t)  to  fire.  The  reason  why  we  require  the  firing  times  to  be  rational 
(and  not  merely  real)  is  that  we  can  discretize  the  processing  times  in 
uni ts  of  time  and  precisely  describe  the  state  of  the  process  at  each 


instant  of  time 


The  rules  of  operation  of  a  TPN  remains  the  same  as  an  ordinary  PN, 
except  that  the  firing  of  a  transition  consists  of  the  three  following 
phases: 

(1)  The  firing  initiation,  that  can  occur  whenever  the  transition  is 

enabled.  When  the  firing  of  a  transition  is  initiated,  one  token  from 

each  of  its  input  places  is  removed, 

(2)  The  firing  execution  which  occurs  during  the  processing  time  p(t). 

(3)  The  firing  termination,  which  occurs  at  the  end  of  the  execution.  At 

this  moment  one  token  is  added  in  each  of  the  output  places  of  the 

transition. 

The  three  phases  of  a  transition  firing  can  be  visualized  by  imagining 
every  transition  as  consisting  of  two  transitions  and  an  intermediate 
place,  as  shown  in  Figure  2.8.  The  firing  time  of  the  transition  t  can  now 
be  asociated  with  the  place  pt  in  the  following  manner:  When  transition  t 
initiates,  tb  fires  instantaneously,  a  token  is  removed  from  each  input 
place  of  tb  and  a  token  is  deposited  on  place  pt.  This  token  stays  in  pfc 
for  the  interval  p(t),  the  firing  time  of  t.  At  the  end  of  this  interval, 
transition  tg  fires,  corresponding  to  the  termination  of  the  transition  t. 
The  substitution  described  in  Figure  2.8  is  actually  the  transformation 
used  by  Sifakis  [14]  to  prove  the  equivalence  of  both  types  of  TPN  Ci.e., 
when  processing  times  are  either  associated  to  transitions  or  places). 

Remark:  With  this  definition  of  a  TPN,  as  given  above,  it  should  be  clear 
that,  at  time  x,  the  state  of  the  Petri  Net  is  not  completely  described  by 
its  marking  M(x):  some  transitions  are  being  processed  at  that  tine,  which 
implies  that  some  tokens  are  not  taken  into  account  in  M(x).  If  we 
substitute  for  each  transition  the  subnet  model  described  in  Figure  2.8, 
the  tokens  being  in  the  execution  phase  would  be  in  the  corresponding 


Figure  2.8  Timed  Petri  Net  Equivalence 

place  pt.  This  remains  however  a  pure  matter  of  convention,  which  does  not 
affect  the  conservation  of  tokens  in  the  net.  In  fact,  in  modeling  the 
dynamic  behavior  of  the  system  and,  more  precisely,  analyzing  how  the 
sequence  of  transition  firings  occurs  in  real  time,  only  the  "free"  tokens 
deserve  attention;  these  tokens  determine  whioh  transitions  are  enabled 
and,  therefore,  which  transitions  can  fire  at  a  given  instant  of  time  x. 

In  order  to  have  a  complete  description  of  the  state  of  the  system  at 
any  instant,  Ramchandani  has  introduced  the  following  two  functions,  which 
will  be  often  referred  to,  in  this  thesis: 

-  Ij(t),  which  denotes  the  number  of  initiations  of  trams! ti on  tj  in 
the  interval  of  time  [0,x] . 

-  Tj(t),  which  denotes  the  number  of  terminations  of  transition  tj  in 
the  interval  [0,t], 

It  becomes  easy  to  determine  the  marking  of  the  net  at  any  instamt, 
using  these  two  functions.  Let  us  write  the  incidence  matrix: 


where  C+  characterizes  the  directed  arcs  that  go  from  transitions  to 
places,  i.e.. 


■  0(V  »i> 


and  C~  characterizes  the  directed  arcs  that  go  from  places  to  transitions. 


i.e. : 


cIj  -  I(W 


Let  M^(t)  denote  the  marking  of  place  p^^  (i=l,...,n)  at  the  instant  of  time 
t.  Then: 


m  m 

V'1  -  <  *  1  <4  V°  -  I  Cl3  !)<'>  (2-,0> 

J-l  J=1 

This  equation  describes  the  marking  of  the  net  at  any  instant,  given  the 
functions  Tj(v)  and  Ij(t)  (J=l,...,m)  and  we  will  make  use  of  it  in  later 
sections. 

In  the  next  section,  we  are  going  to  apply  the  model  of  TPN  to  the 
special  class  of  Petri  Nets,  namely,  Event-Graphs.  The  resulting  Timed 
Event  Graph  will  actually  be  the  model  used  to  carry  out  the  performance 
evaluation  of  the  DMO. 

2.4.2  Timed  Event-Graph 

A  Timed  Event-Graph  (TEG)  is  simply  a  model  of  TPN,  as  presented  in 
the  last  section,  i.e.,  a  pair  (PN,p)  where  PN  is  in  this  case  an  Event- 
Graph  .  The  model  of  TEG  used  in  this  thesis  corresponds,  however,  to  the 
one  defined  by  Chretienne  and  Carller  [16],  which  includes  an  additional 
condition  in  the  firing  process  of  a  transition.  In  the  definition  of  a 


TPN,  a  firing  can  be  initiated  as  soon  as  the  corresponding  transition  is 
enabled.  This  implies,  in  particular,  that  firings  can  be  initiated  even 
if  the  transition  is  already  executing.  In  other  words,  no  assumptions  are 
made  regarding  the  capacity  of  place  pt  (Figure  2.8),  which  therefore  can 
contain  more  than  one  token  at  a  time.  This  is  not  allowable  in  the 
modeling  of  the  decisionmaking  process  developed  in  this  thesis.  To  the 
extent  that  tokens  represent  information  messages  and  transitions  model  the 
different  processing  stages,  no  more  than  one  information  message  can  be 
processed  in  a  certain  stage  at  a  time.  In  other  words,  inputs  are 
processed  one  by  one  at  any  stage.  This  constraint  implies,  in  fact,  that 
a  transition  is  not  allowed  to  initiate,  if  it  is  already  executing. 

This  contraint  can  actually  be  modeled  in  the  Petri  Net  representation 
by  adding  a  self-loop  to  each  transition  of  the  net,  as  shown  in  Fig.  2.9. 


Figure  2.9  Representation  of  a  Self-loop 

The  fact  that  there  is  only  one  token  in  the  extra  place  pg  accounts 
for  this  constraint.  The  same  kind  of  transformation,  as  used  previously, 
describes  clearly  what  happens  in  the  firing  process  (Fig.  2.10).  Onoe  t,. 
has  fired,  one  token  is  in  pt  and  place  ps  is  empty.  Therefore,  transition 
tb  cannot  fire  again  (since  it  is  not  enabled)  until  te  has  fired, 
which  occurs  at  the  end  of  the  interval  }i(t)  (i.e.,  the  firing  time  of 
transition  t) . 

Although  self-loops  allow  for  a  clear  representation  of  the  firing 
constraint  in  the  Petri  Net  model,  they  are  not  really  convenient  for  two 


• v"  v  v 

»  *.  »  /«  -  j»  ■  *  *  »  >  ■ 


Figure  2.10  Model  of  Transition  Firing 


reasons:  First,  they  do  not  change  the  structure  of  the  system,  but 
increase  its  complexity,  since  the  number  of  places  is  considerably 
greater.  Second,  they  cannot  be  included  in  the  incidence  matrix,  because 
the  structure  of  the  incidence  matrix  only  permits  to  represent  the 
presence  (or  absence)  of  one  directed  arc  between  any  two  nodes.  For  these 
reasons,  self-loops  will  not  be  represented  in  our  Petri  Net  model  of  DMOs. 
They  will  be,  however,  implicit  and  it  will  be  important  to  keep  them  in 
mind,  when  analyzing  the  real-time  process.  Indeed,  we  will  see  that  the 
circuits  of  the  net  play  a  orucial  role  in  the  time-related  performance 
measures  of  the  system.  At  this  point,  each  transition  should  also  be 
considered  as  a  circuit  containing  one  token,  precisely  because  of  the 
implicit  self-loop. 
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CHAPTER  III 


PETRI  NET  DESIGN  OP  DECISIONMAKING  ORGANIZATIONS 

In  this  chapter,  we  develop  some  extensions  to  the  Petri  Net 
representation  of  the  DMO,  as  previously  introduced  in  [2].  Included  in 
our  modeling  are  the  resource  and  time  constraints  of  the  DMO.  The 
resource  constraints  concern  both  individual  DMs  and  the  organization  as  a 
whole.  This  type  of  constraint  is  modeled  as  a  limitation  on  the  number  of 
inputs  that  can  be  bandied  at  the  same  time.  Time  constraints  occur  from 
the  various  task  processing  times.  Firing  times  sire  therefore  assigned  to 
the  transitions,  so  that  the  representation  corresponds  to  a  Timed  Petri 
Net  model. 

3.1  AGGREGATED  MODEL  OF  THE  INTERACTING  DECISIONMAKER 

The  modeling  of  decisionmaking  organizations  using  Petri  Nets  has  been 
introduced  in  [21.  Figure  3.1  shows  the  aggregated  Petri  Net  model  of  the 
single  interacting  decisionmaker.  The  decision  process  occurs  in  four 

stages:  Situation  Assessment  (SA),  Information  Fusion  (IF),  Command 

Interpretation  (Cl)  and  Response  Selection  (RS).  Incoming  inputs  -  either 

from  the  environment  or  from  other  DMs  -  are  processed  at  the  first  stage 
to  produce  the  situation  assessment.  This  information  is  then  fused  at  the 
IF  stage  with  situation  assessments  oommunicated  by  other  DMs.  The 

resulting  information  is  combined  with  commands  received  by  the  DM  in  the 
Cl  stage,  so  as  to  select  a  response  in  the  RS  stage.  Each  of  the  four 

stages,  which  corresponds  to  a  particular  task  performed  by  the  DM,  is 

modeled  by  a  transition.  The  firing  of  a  transition  represents  the 

execution  of  the  corresponding  stage  in  the  process.  According  to  the 
rules  of  operation  of  a  Petri  Net,  a  transition  can  fire  when  it  is 

enabled,  i.e.,  when  there  is  at  least  one  token  In  each  of  its  input 
places.  When  a  transition  fires,  one  token  is  removed  from  each  of  its 
input  places  and  one  token  is  added  in  each  of  its  output  places. 


*.  s 


Figure  3.1  Model  of  the  Interacting  Declalonaaker 


Assuming  that  tokens  represent  Information  messages,  the  flow  of  tokens 
models  the  flow  of  Information  In  the  process.  Places  reoelve  and  transmit 
the  information  that  is  exchanged  between  the  different  processing  stages. 

If  the  role  of  transitions  Is  dearly  defined,  since  they  represent 

the  various  processing  tasks,  the  role  of  places  needs  to  be,  however, 

further  specified.  From  Figure  3.1,  we  can  notice  that  places  carry  in 
fact  two  types  of  information: 

-  places  like  px,  pf,  p4,  pT  and  p,  carry  the  information  being 
exchanged  between  the  DM  and  the  environment  or  other 
organizational  members. 

-  places  like  pa,  p,  and  p4  carry  the  Information  that  is 

internally  processed  by  the  DM. 

The  first  type  of  places  allows  for  modelling  the  interactions  between 
the  DM  and  the  environment  or  other  DMs:  they  are  either  inputs  to  the  DM 
(like  px,  pT,  pa)  or  outputs  of  the  DM  (like  pf,  p(),  but  in  no  case  can 
they  be  both  at  the  same  time.  If  we  eonslder  these  places  for  all  the  DMs 
Involved  in  the  organization,  they  determine  in  faot  the  organisational 
structure,  l.e.,  the  Internal  interactions  between  the  DMs  and  the 

interactions  between  the  DMs  and  the  environment.  This  is  the  reason  why 
we  will  call  them  structural  places,  in  contrast  to  the  other  places. 
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The  second  type  of  places,  together  with  the  four  tasks  defined  above, 
model  the  internal  structure  of  the  interacting  DM.  These  places  carry  the 
information  that  is  produced  at  one  stage  of  the  decisionmaker  process  and 
used  at  the  next  stage.  In  contrast  to  the  structural  places,  these  are 
both  input  and  output  places  of  the  same  DM.  At  this  point,  we  should  note 
that  the  processing  of  a  specific  input  takes  place  in  an  asynchronous 
manner,  i.e.,  delays  generally  occur  between  the  different  processing 
stages,  precisely  because  of  the  interactions  with  the  other  DMs  or  the 
environment.  For  example,  once  the  situation  assessment  stage  is 
completed,  one  token  is  in  pa,  but  the  information  fusion  cannot  take  place 
before  one  token  is  also  in  pT,  which  is  the  information  transmitted  from 
the  SA  stage  of  another  DM.  Such  delays  imply  that  information  has  to  be 
stored  temporarily  in  places  pa,  p,  or  p4  until  the  processing  of  the 
corresponding  stage  can  occur.  This  is  the  reason  why  we  call  these  places 
memory  places:  they  model  the  Internal  (short-term)  memory  of  the  DM, 
where  the  Information  being  processed  has  to  be  stored  temporarily  until 
the  different  processing  stages  are  completed. 

3.2  MODEL  OF  THE  INTERACTING  DECISIONMAKER  WITH  LIMITED  RESOURCES 

The  model  of  the  Interacting  Decisionmaker  does  not  take  into  acoount 
the  limited  capacity  for  information  processing,  that  characterizes  the 
human  DM.  Indeed,  as  long  as  information  messages  are  present,  i.e., 
tokens  are  available  in  place  px,  the  processing  can  start,  i.e.,  the 
transition  corresponding  to  the  SA  stage  can  fire.  However,  the 
information  processing  of  human  DMs  is  subject  to  the  bounded  rationlity 
constraint,  as  defined  in  14] .  In  the  information- theoretio  approach,  the 
amount  of  information  processed  is  measured  by  the  total  aotlvity,  G,  of 
the  DM,  which  also  characterizes  the  DM's  workload.  It  is  assumed  that 
there  exists  an  upper  bound  Gr,  above  which  the  DM  becomes  overloaded  and 
his  performance  degrades: 


When  the  analysis  is  carried  out  for  the  steady-state  process,  the  above 
constraint  takes  the  form: 

Cl  FfG  (3.2) 

r 

where  F  is  the  processing  rate  constraint  that  characterizes  the  human  DM, 
and  1/t  (resp.  x)  is  the  average  arrival  rate  of  inputs  (resp.  average 

interarrival  time).  This  constraint  implies  that  the  DM  must  process 

inputs  at  a  rate  at  least  equal  to  the  rate  at  which  they  arrive. 

In  our  work,  we  deal  with  time-related  performance  measures  and  do  not 
address  the  accuracy  of  the  response,  which  is  based  on  the  comparison 
between  the  actual  organization's  response  and  the  ideal  or  desired 
response  [1].  In  particular,  the  way  a  DM  reacts  to  Information  overload 
(which  occurs  when  G  >  F  t)  and  the  extent  to  which  it  affects  his 

performance  are  not  a  matter  of  concern  here.  This  remains  however 
allowable,  in  so  far  as  we  also  include  in  our  modelling  the  actual 

processing  constraint,  that  can  be  expressed  by  writing  inequality  (3.2)  in 
another  form: 

1/r  <  (3.3) 

Using  relation  (3.3),  the  bounded  rationality  limitation  turns  out  to  be, 
in  that  case,  a  constraint  on  the  allowable  rate  of  incoming  inputs,  i.e., 
on  the  maximum  rate  of  inputs  that  can  be  bandied  by  the  DM,  without  being 
overloaded.  In  other  words,  rather  then  assuming  that  Inputs  are  prooessed 
at  a  rate  at  least  equal  to  the  rate  with  which  they  arrive,  and  derive  a 
bound  on  the  total  activity  allowable,  G,  (given  F),  we  assume  that  no 
information  overload  does  occur  in  the  process  and  derive  the  corresponding 
constraint  on  the  allowable  rate  of  incoming  inputs. 

Let  us  see  now  how  it  is  possible  to  model  this  constraint  using  the 
Petri  Net  framework.  In  fact,  the  limited  processing  capabilities  of  a 


human  DM  come  from  the  limited  (cognitive)  resources  available  to  perform 
the  various  processing  tasks.  In  particular,  the  bounded  rationality 
constraint  is  very  much  related  to  the  limited  capacity  of  the  human  short¬ 
term  memory,  defined  as  the  memory  in  which  the  information  is  held 
temporarily  (in  contrast  to  the  long-term  memory,  in  which  the  information 
is  stored  permanently).  Indeed,  this  bound  means  that  a  DM  cannot  handle 
properly  too  much  information  at  the  same  time.  Interestingly  enough,  the 
analysis  of  the  human  short-term  memory  carried  out  in  [5]  has  shown  that  a 


maximum  of  six  or  seven  units  of  information  can  only  be  held  in  this 


memory  (typically  for  a  few  seconds),  without  any  loss  of  information. 
This  is  of  much  importance  here,  because  the  DM  has  precisely  to  handle 
different  information  messages  during  the  time  necessary  to  complete  the 
various  processing  tasks.  This  is  the  reason  why  we  have  identified,  in 
the  previous  section,  the  role  of  places  pa,  pa  and  p4  (Figure  3.1)  as 
short-term  memory:  these  places  contain  the  Information  messages 
(represented  by  tokens),  that  are  being  processed  internally  by  the  DM. 
Quite  naturally,  the  limited  capacity  of  the  short-term  memory  can 
therefore  be  modeled,  using  the  Petri  Net  formalism,  as  a  capacity 
constraint  on  the  corresponding  places  pa,  pf  and  p4.  Accordingly,  we 
extend  the  model  of  Interacting  Decisionmaker  as  shown  in  Figure  3.2. 


Figure  3.2  Model  of  Interacting  Decisionmaker  with  Limited  Resources 
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Now  the  added  place  R1#  whose  Input  transition  corresponds  to  the  RS 
stage  and  output  transition  to  the  SA  stage,  allows  for  modeling  the 
Information-processing  constraint.  Ve  shall  call  it  the  reaouroe  plaoe, 
because  the  number  of  tokens  put  initially  in  this  place  represents  the 
resources  available  for  processing.  Indeed,  the  directed  path: 

P  *  pa  t,  p,  t,  p4  t4) 

determines  a  circuit  in  the  net,  in  which  the  token  content  remains 
invariant  by  transition  firings  (Chapter  2,  Theorem  2.3),  i.e.,  for  any 
reachable  marking  M: 

4  4 

M(Rx)  +  ^  M(PjL)  =  M#(Ri)  +  J  M*(Pi)  =  H°(Ri)  (3.4) 

1=2  i=2 

It  is  assumed  here  that  places  px,  p,  and  p4  contain  initially  no 
tokens,  whloh  means  that,  initially,  there  are  no  information  messages 
being  processed.  Now,  the  constraint  on  the  memory  capacity  is  trivially 
satisfied.  Assuming  that  there  are  n  resources  available  (i.e.,  M*(RX)  =  n) 
which  may  represent  the  amount  of  memory  space  available  (in  other  words, 
the  cognitive  resources),  we  deduce  indeed  from  the  marking  equation  (3.4), 
that,  at  any  reachable  state  of  the  process: 

H(pa)  +  M(pf)  +  M(p4)  <  n 

which  precisely  models  the  short-term  memory  limitation,  as  described 
earlier . 

At  this  point,  it  is  possible  to  extend  further  our  analysis  of  the 
model,  especially  with  respect  to  the  processing  rate  constraint. 
Let  us  specify  now  how  the  processing  of  inputs  does  occur,  bolding 
the  same  assumptions  about  the  initial  marking  of  each  place,  i.e. 
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The  processing  of  any  new  input  starts  by  the  firing  of  transition  t4, 
(corresponding  to  the  SA  stage),  which  consumes  therefore  one  token  (i.e., 
one  resource)  from  place  R4 .  Likewise,  when  the  processing  is  completed, 
transition  t4  (corresponding  to  the  RS  stage)  fires,  which  produces  one 
token  back  in  R4.  In  other  words,  one  resource  is  engaged  at  the  beginning 
of  the  process  and  is  released  at  the  end  of  the  process.  Clearly  enough, 
the  DM  can  only  process  at  most  n  different  inputs  at  the  same  time.  If  n 
inputs  are  currently  being  processed,  the  resource  place  R4  is  empty  (which 
follows  directly  from  the  balance  equation  (3.4))  and  transition  t4  is 
consequently  not  enabled.  Otherwise  stated,  the  resource  constraint  turns 
out  to  bound  the  number  of  inputs  that  the  DM  can  process  simultaneously. 

To  analyze  deeper  how  it  affects  the  processing  rate,  it  is  necessary 
to  introduce  the  processing  times  of  the  four  stages  SA,  IF,  Cl,  and  RS, 
that  will  be  denoted  respectively  by  | ix,  p4,  p,  and  p4.  Let  us  point  out 
that  a  complete  description  of  the  model,  with  firing  times  associated  to 
the  transitions,  will  be  presented  later  in  this  chapter.  Our  intention 
here  is  only  to  make  clear  how  the  resource  limitation  actually  bounds  the 
processing  rate.  Now,  if  the  decision  process  were  fully  synchronous, 
which  means  that  no  delays  would  occur  between  the  different  processing 
stages,  it  would  take  the  amount  of  time 

»*.  -  +  N,  +  N,  + 

to  complete  the  processing  of  any  input.  Since  the  DM  can  handle  at  most  n 
inputs  at  the  same  time,  the  processing  rate  is  necessarily  bounded  by: 


There  is,  however,  am  additional  constraint  to  include:  the  ezeoution  of 
any  processing  stage  can  only  take  place  for  only  one  input  at  a  tine  as 
discussed  in  Chapter  2,  Section  2.4.2.  Therefore,  the  actual  bound  is 
determined  by: 

p  =  min  (f,  —  ,  —  ,  —  ,  —  )  (3.5) 

u  u  u  u 

rl  *2  *»  r« 

We  will  show  in  the  next  chapter  that  p  determines  precisely  the  maximum 
rate  of  information  processing  that  characterizes  the  DM.  This  rate 
constraint  is  derived  from  the  limited  capabilities  of  a  human  DM,  whloh 
Include  both  the  resource  limitations  (as  specified  by  n)  and  the 

processing  time-constraints  (as  specified  by  the  amount  of  time  required  to 
perform  the  various  tasks),  p  is  in  fact  similar  to  the  processing  rate 
constraint  F,  introduced  in  the  information  theoretic  framework. 

Remark:  We  have  extended  the  Petri  Net  model  of  the  Interacting 

Decisionmaker,  so  as  to  include  the  (cognitive)  resource  oonstraints, 
derived  from  the  limited  capacity  of  the  human  short-term  memory.  However, 
the  model  is  also  suitable  for  any  kind  of  resouroe  limitations  that 

oonstrain  the  processing  of  inputs.  Suppose,  for  example,  that  the  same 
specific  resource  is  needed  to  complete  the  processing  of  any  input. 
Putting  initially  only  one  token  in  the  resource  plaoe  Rx  allows  for 
modeling  the  corresponding  constraint. 

Moreover,  as  explained  earlier,  the  resource  limitation  turns  out  to 
be  a  constraint  on  the  number  of  inputs  that  can  be  bandied  at  the  same 
time.  Therefore,  the  resources  available  to  the  DM  may  have  in  reality 
various  forms,  but  it  is  only  necessary  to  take  into  account,  in  the 

modeling,  the  most  stringent  resource  limitation.  If,  for  instance,  the  DM 


has  the  capacity  to  handle  several  inputs  simultaneously,  but  needs  a 
unique  specified  resource,  the  corresponding  processing  constraint  will  be 
modeled  by  putting  only  one  token  in  Rx,  as  explained  above. 


3.3  AGGREGATED  MODEL  OF  THE  DM0 

After  having  analyzed  the  model  of  the  simple  Interacting  Decision¬ 
maker,  we  present  in  this  section  the  aggregated  Petri  Net  model  for  the 
overall  organization.  At  first,  we  deal  with  the  interactions  between  the 
organization  and  the  environment. 

3.3.1  Model  of  the  DM0  Interacting  with  the  Environment 

The  interactions  between  a  DM0  and  the  external  environment  has  been 
described  in  [1].  The  first  processing  stage  of  the  DM0  consists  of  the 
partitioning  of  the  external  signals  or  messages  into  a  set  of  inputs  that 
are  assigned  to  different  organization  members.  This  particular  type  of 
allocation  has  been  addressed  in  [17].  Assuming,  for  example,  that  DM1,  2, 
and  3  are  the  DMs  that  interact  directly  with  the  environment,  the 
corresponding  Petri  Net  representation  is  shown  on  Figure  3.3. 


Figure  3.3  Model  of  the  Partitioning  Operation 
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The  place  ps  represents  the  source  of  signals  or  messages,  and  the 
transition  tx  models  the  partitioning  operation.  Since  it  is  the  first 
processing  stage,  tx  will  be  called  the  input  transition  of  the  process. 
Each  time  that  tx  fires,  one  token  is  sent  in  the  plaoes  px,  p,,  and  p,, 
which  are  the  input  places  of  DM1,  2  and  3.  It  should  be  clear  that  this 
model  implies  an  overall  synchronization  between  the  individual  inputs 
received  by  each  of  the  DMs.  In  the  following  stages,  however,  the 
processing  of  the  inputs  by  each  DM  becomes  asynchronous  and  concurrent,  as 
emphasized  in  the  introduction. 

Likewise,  the  Petri  Net  representation  of  the  DMO  has  an  output 
transition,  which  characterizes  the  last  processing  stage.  The  output 
place  of  this  transition  contains  the  organizational  responses.  For  the 
moment,  we  will  denote  this  transition  by  tffl,  where  m  is  the  total  number 
of  transitions  in  the  net.  pr  will  denote  this  output  place.  Finally,  the 
aggregated  representation  of  the  DMO  interacting  with  the  environment  is 
shown  in  Figure  3.4. 


Figure  3.4  Model  of  DMO  Interacting  with  the  Environment 

Interestingly  enough,  we  can  imagine  the  overall  process  as  taking 
place  through  a  sequence  of  three  major  stages.  Each  of  the  stages  would 
correspond  to  a  particular  action,  as  defined  in  [3],  as  follows: 


-  the  input  action,  corresponding  to  the  partitioning  operation  of 


the  external  messages 


-  the  processing  of  the  Information,  carried  out  by  the  different 
DMs  Involved  in  the  organization) 

-  the  output  action,  corresponding  to  the  production  of  the  overall 
response  that  characterizes  the  organizational  response. 

3.3.2  Model  of  the  DM0  with  Limited  Resources 


The  analysis  of  the  resource  constraint,  carried  out  in  Section  3.2 
for  the  Interacting  Decisionmaker,  can  be  applied  in  fact  to  the 

organization  as  a  whole.  Indeed,  the  resources  used  by  the  overall 
organization  may  have  various  forms,  but  there  exists  always  at  least  a 

processing  constraint  that  comes  from  the  limited  capacity  of  the 

structural  places.  This  constraint  is  very  similar  to  the  one  existing  for 
the  memory  places  of  the  individual  DMs. 

As  described  in  Section  3.1,  the  structural  places  reoeive  and 
transmit  the  information  between  the  DMs  and  between  the  DMs  and  the 

environment;  Now,  since  the  process  is  asynchronous,  information  is  to  be 
stored  temporarily  in  these  places,  exactly  as  what  ooours  for  the  internal 
decision  process  of  a  DM.  These  places  oarry  in  fact  the  role  of  buffer 
storage  within  the  organization,  where  some  information  is  stored 
temporarily.  It  is  therefore  important  to  make  sure  that  at  no  instant 
does  the  amount  of  information  stored  exceed  the  buffer  capacity  of  the 
system,  i.e.,  the  capacity  of  the  structural  places. 

Let  us  assume,  for  instance,  that  the  arrival  rate  of  external  inputs 
exceeds  the  maximum  processing  rate  of  any  one  of  the  DMs  (defined  as  p  in 
Section  3.2).  With  the  present  Petri  net  model  of  the  DM0,  it  follows  that 
tx  will  fire  at  the  same  rate  with  which  inputs  arrive:  Tokens  will 
therefore  necessarily  accumulate  in  the  system  and,  in  particular,  the 
token  content  of  some  structural  places  of  the  net  will  eventually  grow  to 
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infinity  over  time:  in  such  case,  we  will  say  that  the  DMO  is  overloaded 


To  handle  this  oontraint,  we  modify  the  Petri  net  model  of  the  DMO  by 
adding  an  extra  place,  exactly  as  we  did  for  the  single  Interacting 
Decisionmaker.  The  resulting  model  is  shown  on  Figure  3.5.  R#  is  the 
resource  of  the  overall  organization:  the  number  of  tokens  put  initially 
in  this  place  bounds  the  number  of  Inputs  that  the  DMO  can  process  at  the 
same  time.  Indeed,  each  time  that  a  new  input  is  processed  (i.e.,  whenever 
tx  fires),  a  token  is  removed  from  R#  and  comes  back  again  in  R#  once  the 
processing  is  completed  (i.e,  when  tm  fires).  Assuming  that  R,  contains 
initially  n  tckens  and  that  n  inputs  are  currently  being  processed,  then  R, 
is  empty  and  no  more  inputs  can  be  processed,  since  tx  is  not  enabled. 


Let  us  point  out  that  the  same  remark,  made  in  Section  3.2,  applies 
here:  it  is  only  worth  taking  into  account  the  type  of  resouroe  limitation 
which  is  the  most  stringent.  If,  for  instance,  the  same  speolfio  resouroe 
is  needed  to  complete  the  processing  of  any  organizational  input,  putting 
one  token  in  R0  will  account  for  this  processing  constraint. 

In  the  last  two  sections,  we  have  dealt  with  a  specific  type  of 
processing  constraints  derived  from  the  limitation  of  resources.  In  the 
next  section,  we  will  address  time  constraints,  as  determined  by  the 
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various  task  processing  times 


3.4  MODELING  TIME  CONSTRAINTS 

Apart  from  the  resource  constraints,  which  were  treated  in  Section  3.2 
and  3.3,  the  other  type  of  processing  constraints  is  time  constraints, 
which  arise  from  the  task  execution  times.  As  already  introduced  in 
Section  3.2,  the  task  processing  times  are  included  in  the  model  by 
assigning  to  each  transition  a  corresponding  firing  time.  At  this  point, 
however,  the  model  of  Timed  Petri  Net  described  in  Chapter  2,  Section  2.4 
where  the  firing  times  are  assumed  deterministic,  cannot  be  applied 
directly  to  the  aggregated  model  of  the  DM0  developed  so  far.  The  problem 
arises  from  switches,  whose  function  is  more  complex  than  ordinary 
transitions. 

In  the  aggregated  model  of  the  Interacting  Decisionmaker  described  in 
Section  3.1,  all  four  processing  stages  SA,  IF,  Cl  and  RS  were  represented 
by  transitions.  In  reality,  the  representation  of  both  the  SA  and  RS 
stages  includes  a  decision  node  or  switch  as  defined  in  [1],  which  is  a 
more  complex  form  of  transition.  An  n-decision  switch  is  in  fact  a  subnet 
that  contains  n  (simple)  transitions,  as  shown  on  Figure  3.6. 


Figure  3.6  Petri  Net  Model  of  an  n-Decislon  Switch 


Each  transition  of  a  switch  ( tgl, tsl, . . . , tsn)  represents  an  alter- 
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native  process  or  algorithm.  The  selection  of  one  algorithm,  i.e.,  the 
firing  of  one  of  these  transitions,  occurs  for  every  input  (i.e.,  token) 
processed.  It  turns  out  that  each  algorithm  may  have  a  different 
processing  time,  which  implies  that  the  processing  time  of  the  SA  (or  IS) 
stage  is  not  necessarily  unique,  but  depends  on  the  algorithm  selected.  In 
fact,  if  in  our  model  we  represented  explicitly  all  the  transitions  of  a 
switch,  it  would  be  possible  to  assign  to  each  transition  of  the  net  a 
deterministic  firing  time,  according  to  the  Timed  Petri  Net  model  presented 
in  Chapter  2. 

However,  as  it  will  be  clear  later,  the  performance  analysis  oarrled 
out  in  this  thesis  depends  crucially  on  the  fact  that  the  DMO  can  be 
modeled  as  an  Event-Graph,  i.e.,  each  place  has  exactly  one  input  and  one 
output  transition;  this  condition  would  be  violated,  if  the  representation 
shown  in  Figure  3.6  is  used.  This  is  why  the  aggregated  model  is  very 
convienent  for  our  purpose.  Now,  to  handle  the  problem  of  swltohes,  we 
will  introduce  the  following  extension  to  the  Timed  Petri  Net  model:  The 
firing  time  of  a  switch  will  be  considered  as  a  discrete  random  variable, 
having  a  finite  number  of  possible  values  (corresponding  to  the  processing 
times  of  the  different  algorithms),  according  to  a  certain  probability 
distribution.  Naturally,  the  probability  distribution  will  depend  on  the 
decision-rule  that  determines  the  selection  among  the  different  algorithms. 
Because  tokens  are  indistinguishable  in  a  Petri  Net,  the  probability 
distribution  will  be  assumed  Independent  from  one  token  to  the  next  (i.e., 
independent  from  one  input  to  the  next) . 

In  the  next  section,  we  summarize  the  key  points  that  have  been 
developed  so  far.  In  particular,  we  show  that  the  model  of  the  DMO,  which 
includes  resource  and  time  constraints,  corresponds  to  a  strongly  connected 
Timed  Event-Graph  (following  the  definitions  given  in  Chapter  2).  We  will 
illustrate  the  development  with  a  two  decisionmaker  organization. 
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3.5  MODEL  OF  THE  DM0  AS  A  STRONGLY  CONNECTED  TIMED  EVENT-GRAPH 

Let  us  recall  that  an  Event-Graph  is  a  Petri  Net  such  that  each  place 
has  exactly  one  input  and  one  output  transition.  The  strong  connectivity 
means  that  there  exists  a  directed  path  between  any  two  nodes  (places  or 
transitions).  Now.  our  aggregated  model  of  the  DM0  is  dearly  an  Event- 
Graph  since,  as  seen  in  Section  3.2  and  3.3,  each  place  of  the  net  — 
whether  it  is  a  structural,  memory  or  resource  place  —  has  indeed  one 
input  and  one  output  transition.  For  the  memory  and  resource  places,  this 
is  due  directly  to  the  internal  structure  of  the  Interacting  Decisionmaker 
(Figure  3.2).  It  is  also  clear  for  the  structural  places,  because  the 
information  messages  (i.e.,  tokens)  are  transmitted  from  a  certain 
processing  stage  to  another  stage  (both  modeled  by  transitions),  as 
determined  by  the  execution  protocol  of  the  organization. 

Naturally,  we  are  interested  in  analyzing  the  performances  of  the 
admissible  organizational  forms,  as  formulated  in  [18] .  Two  constraints 
should  in  particular  be  satisfied: 


v. 


(1)  the  information  structure  should  have  no  loops 

(2)  a  directed  path  should  exist  from  the  input  node  to  every  node 
and  a  directed  path  should  exist  from  any  node  to  the  output 
node. 

At  this  point,  it  is  quite  important  to  clarify  what  the  first 
contraint  implies  in  our  model.  This  constraint  means  that  the  information 
structure  should  be  acycllcal,  i.e.,  there  should  be  no  Information  loops. 
Of  course,  in  our  Petri  Net  model  of  the  DM0,  there  are  loops  (or 
circuits),  because  of  the  resource  places.  But  the  loops  that  exist 
because  of  the  resource  places,  i.e.,  the  loops  that  include  a  least  one  of 
the  resource  places,  do  not  actually  violate  the  constraint.  As  far  as  our 
model  is  concerned,  the  precise  formulation  of  this  constraint  is  the 
following : 


>•3 


(1)  Only  circuits  which  contain  one  or  more  of  the  resource  places 
are  allowed.  In  other  words,  if  we  delete  from  the  Petri  Net 
all  the  resource  places,  the  resulting  net  should  have  no 
circuits  (the  net  so  obtained  models  the  information  structure, 
as  developed  in  [1]). 

Now,  the  second  constraint  makes  sure  that  our  Petri  Net  model  is 
strongly  connected,  precisely  because  of  the  resource  place  of  the  overall 
organization,  R#.  As  we  have  seen,  this  resource  place  connects  the  output 
transition  to  the  input  transition  of  the  net  (Figure  3.5).  Let  us  now 
oonsider  any  two  nodes  of  the  net.  There  exists  a  directed  path  from  the 
first  node  to  the  output  transition  and  therefore  to  the  input  transition 
(via  the  resource  place).  There  exists  also  a  directed  path  from  the  input 
transition  to  the  second  node.  By  appending  these  two  paths,  we  have 
constructed  a  directed  path  between  the  two  nodes. 

Finally,  given  the  processing  times  assigned  to  the  transitions,  as 
described  in  Section  3.4,  the  Petri  net  model  of  the  DMO  turns  out  to  be  a 
strongly  connected  Timed  Event-Graph.  We  illustrate  the  model  for  two 
different  structures  of  a  two  member  organization,  as  shown  in  Figures  3.7 
and  3.8.  The  example  is  derived  from  the  analysis  carried  out  in  [18]. 
Transitions  t,,  t,,  t4,  and  t,  (reap,  t,,  tT,  t#,  and  tt)  correspond 

respectively  to  the  SA,  IF,  Cl,  and  RS  stages  of  DM1  (resp.  DM2).  In  both 
cases,  DM2  transmits  the  result  of  his  SA  stage  to  mi,  which  constitutes 
information  sharing  among  them.  However,  in  structure  (a),  information  is 
transmitted  from  the  RS  stage  of  DM1  to  the  IF  stage  of  DM2  and  this 
interaction  is  of  the  result  sharing  type,  as  defined  in  [18].  In 
structure  (b)  the  same  information  is  transmitted  to  the  Cl  stage  of  DM2, 
which  corresponds  in  that  case  to  a  command .  This  implies  that  there  is, 
in  that  case,  a  hierarchical  relationship  between  the  two  DMs.  Naturally, 
we  have  added  to  the  Petri  Net  model  all  the  resource  places.,  i.e.,  R#  for 
the  organization  and  Rx  (resp.  R,)  for  DM1  (resp.  DM2),  in  accordance  with 
the  analysis  carried  out  in  this  chapter.  Both  models  will  be  used  to 
illustrate  the  performance  analysis  that  will  be  developed  in  the  next 
chapter. 


CHAPTER  IV 


MAXIMUM  THROUGHPUT  RATE  OF  THE  DMO 


Using  the  model  of  the  DMO  developed  in  the  previous  chapter,  we 
present  in  this  chaper  a  method  for  computing  the  maximum  throughput  rate 
of  the  DMO.  This  MOP  characterizes  the  maximum  rate  of  processing  of  the 
overall  organization.  This  measure  is  especially  important  because  it 
bounds  the  allowable  rate  of  external  inputs  that  can  be  handled  by  the 
DMO.  Given  the  internal  structure  of  the  system,  the  maximum  throughput 
rate  is  characterized  as  a  function  of  both  the  time  constraints  (i.e.,  the 
various  task  processing  times)  and  the  resource  contraints.  It  will 
therefore  be  possible  to  investigate  which  constraints  should  be  relaxed  in 
order  to  improve  this  performance.  Additionaly,  different  organizational 
structures  can  be  compared  with  respect  to  their  maximum  throughput  rate. 
These  types  of  performance  anal/sis  will  be  illustrated  in  this  chapter 
using  the  example  of  the  two  decisionmaker  organization  presented  in 
Chapter  3. 

4.1  PROPERTIES  OF  THE  PETRI  NET  MODEL 

Before  computing  the  throughput  rate,  it  is  first  necessary  to  verify 
that  the  Petri  Net  model  of  the  DMO  satisfies  the  properties  of  a  well- 
defined  system,  namely  liveness  and  boundedness.  As  stated  in  Chapter  2, 
Section  2.1.6,  these  properties  guarantee  the  system  to  be  deadlock  free 
(liveness)  and  that  information  does  not  accumulate  in  the  system 
(boundedness),  in  other  words  that  the  flow  of  outputs  equals  the  flow  of 
inputs. 

Let  us  recall  that  an  Event-Graph  is  live,  iff  the  token  content  of 
every  directed  circuit  is  strictly  positive  (Theorem  2.4)  and  it  is 
bounded,  iff  the  net  is  strongly  connected  (Theorem  2.5).  Now,  we  have 
assumed  in  thge  previous  chapter,  that  the  model  of  the  DMO  is  admissible. 


which  implies,  that  there  should  not  exist  any  directed  circuit  which  does 
not  contain  any  of  the  resource  places  (otherwise  the  information  structure 
would  not  be  acyclioal).  Therefore,  every  directed  oircuit  contains  at 
least  one  resource  place  and,  consequently,  its  token  oontent  is  strictly 
positive.  Since  the  model  is  also  strongly  connected,  as  seen  in  Chapter 
3,  Seotion  3.5,  the  system  is  indeed  live  and  bounded. 

The  liveness  and  boundedness  properties  are  very  important  for 
analyzing  the  steady-state  process,  because  these  characteristics  imply 
also  that  the  system  is  consistent.  Actually,  an  Event-Graph  is  said  to  be 
consistent,  if  there  exists  a  cyclic  firing  sequence,  i.e.,  a  sequence 
which  fires  each  transition  once  and  brings  the  marking  back  to  its  initial 
state.  Thus,  if  a  system  is  consistent,  it  goes  back  to  its  initial  state 
after  each  cycle  and  then  repeats  itself.  If  the  system  is  inconsistent, 
it  can  be  shown  that  either  it  produces  an  infinite  number  of  tokens  (i.e., 
needs  infinite  resources)  or  consumes  tokens  and  eventually  comes  to  a 
stop.  Consistency  is,  therefore,  a  critical  property  of  systems  which 
should  function  continuously,  with  a  finite  amount  of  resources,  as  it  is 
the  case  for  the  model  of  the  DHO  studied  here.  As  proved  in  [6],  a  live 
and  bounded  Event-Graph  is  consistent.  Having  assured  that  our  Petri  Net 
model  of  the  DMO  is  consistent,  we  can  now  analyze  the  throughput  rate  of 
the  corresponding  system. 

4.2  COMPUTATION  OF  THE  MAXIMUM  THROUGHPUT  RATE 

Once  we  are  sure  that  the  model  of  the  DMO  under  consideration 
corresponds  to  a  live  and  consistent  system,  and  therefore  each  transition 
can  fire  repeatedly,  we  would  like  to  determine  the  maximum  rate  the 
transitions  can  actually  fire.  Indeed,  this  will  yield  the  maximum 
throughput  rate.  Of  course,  this  supposes  that  we  take  explicitly  into 
account  the  transition  firing  times.  At  this  point,  our  analysis  will  be 
based  on  the  model  of  Timed  Petri  Nets  and  more  precisely  on  Timed  Event- 
Graph  described  in  Chapter  2,  Section  2.4.  Some  of  the  results  obtained 
here  have  been  originally  studied  by  Ramchandani  [6],  but  in  a  more 
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restrictive  case.  In  bis  work,  it  was  assumed  that  the  steady-state 
process  was  deterministic  (i.e.,  all  the  transition  firing  times  were 


deterministic)  and 


periodic,  which  means  that  every  transition 


fires  at  regular  intervals  of  time,  with  the  same  period.  No  such 
assumptions  are  made  in  our  analysis,  tfe  can  handle  both  oases:  firing 
times  that  are  deterministic  or  discrete  random  variables  (with  a  finite 


set  of  possible  values),  according  to  the  operating  rules  of  switches,  as 
defined  in  Chapter  3,  Section  3.4.  In  addition,  we  will  deal  in  our 
analysis  with  average  firing  rates,  rather  than  fixed  rates.  Indeed,  as  it 
will  be  shown  in  the  next  chapter,  the  steady-state  process  is  generally 
not  strongly  periodic  even  for  the  deterministic  case:  the  periodicity  is 
of  a  more  complex  type,  that  will  be  described  later  as  the  K-periodicity 
Introduced  by  Chretienne  [16],  Finally,  the  only  assumption  made  in  our 
analysis  is  that  inputs  are  supposed  to  queue  at  the  entry  of  the  system  so 


that  the  process  occurs  repetitively  over  time. 


4.2.1.  Average  Cycle  Time  of  Transitions 


In  this  section,  we  present  the  first  important  and  rather  intuitive 
result  that  the  average  firing  rate  is  the  same  for  all  the  transitions. 
The  average  rate  considered  here  is  very  precisely  defined.  We  assume  that 
the  processing  starts  at  x  =  0  and  then  occurs  repetitively,  as  said 
earlier.  We  then  denote  by  S°  the  Instant  of  time  at  which  the  transition 
tA  initiates  its  n~th  firing,  corresponding  to  the  n-th  occurence  of  this 
particular  processing  task.  Naturally,  (S°  -  sj-x)  represents  the  time 
that  has  elapsed  between  the  (n-l)-th  occurence  and  the  n-th  occurence  of 
transition  tA.  This  can  also  be  interpreted  as  a  cycle  time,  since  the 
system  is  assumed  to  function  cyclically,  each  cyle  corresponding  in  fact 
to  the  complete  processing  of  one  input.  The  definition  of  the  average 


cycle  time  is  therefore  as  follows: 


Definition:  The  average  cycle  time,  denoted  by  0^,  of  transition  t^  is 


defined  as: 


0 


i 


lim 

n-»+® 


(4.1) 


This  is  naturally  an  average  measure,  since  the  time  interval  between  two 
successive  occurences  of  transition  t^  is  averaged  over  the  total  number  of 
occurences  (from  the  initial  instant).  Relation  (4.1)  can  actually  be 
written  in  a  more  simple  form: 


e 


i 


lim 

n->+® 


i  i 

lim  —  -  lim  — 

n-»  +®  n  n-»  +®  n 


and  being  a  constant  value: 


lim  — =■  =  0 
_*  ,  n 
n->+® 


so  that: 


0 


i 


lim 
n  ->+• 


(4.2) 


This  will  be  the  expression  used  later  in  the  developments.  Naturally, 


represents  the  average  firing  rate  of  transition  t ^ 
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We  are  now  going  to  prove  the  following  result: 


Theorem  4.1:  All  the  transitions  have  the  same  average  cycle  time  and 
consequently  the  same  average  firing  rate. 

As  it  will  be  shown,  this  result  holds  because  it  is  possible  to  find, 
for  any  two  transitions,  a  directed  circuit  that  Includes  them.  The  result 
is  therefore  crucially  dependent  on  the  fact  that  the  Petri  net  model  of 
the  DMO  corresponds  to  a  strongly  connected  Event-Graph,  as  stated  in 
Chapter  3,  Section  3.5. 

Proof :  Let  t^  and  tj  be  any  two  transitions.  There  exists  at  least  one 

directed  circuit  that  contains  both  transitions,  because  the  net  is 
strongly  connected:  Indeed,  we  can  find  a  directed  path  from  tA  to  tj  and 
from  tj  to  t^.  Let  us  point  out  that  the  directed  circuit  is  not 
necessarily  elementary  (as  defined  in  Chapter  2,  Section  2.3.1),  but  this 
would  not  change  the  nature  of  the  proof  given  here:  In  that  case  we  could 
split  the  circuit  into  a  set  of  elementary  circuits.  Without  loss  in 
generality,  we  can  therefore  denote  by 


the  corresponding  circuit,  where: 


t 


and 


t 


J 


The  circuit  is  shown  on  Figure  4.1.  Because  the  net  is  an  Event-Graph, 
each  place  of  the  circuit  has  exactly  one  input  and  one  output  transition 
and  the  token  content  of  the  circuit  is  invariant  by  transition  firings 
(Theorem  2.3,  Chapter  2,  Section  2.3.2).  Let  m£  denote  the  initial  marking 
of  place  pk  and  M^j  (resp.  M^)  the  number  of  tokens  in  the  places 
belonging  to  the  path  from  tA  and  tj  (resp.  from  tj  to  t^),  i.e.. 


k-1  r 

o  V  o  o  v  o 

Mu  "  2  Hi  and  Mji  ■  1  Hi 

P-1  p  p-k  p 

Let  us  recall  from  Chapter  2,  Section  2.4.1,  that  we  denote  by: 

IA(t)  (reap.  I j (x) ) ,  the  number  of  initiations  of  transition  t±  (resp. 
tj)  in  the  interval  of  time  [0,t], 

T^(t)  (resp.  Tj(t)),  the  number  of  terminations  of  transition  t^ 
(resp.  tj)  in  the  interval  [0,t]. 

tir 


tj  = 


til  =  ti 


Figure  4.1  Directed  Circuit  Containing  tA  and  tj 

The  processing  is  assumed  to  start  at  t=0  (no  transitions  are 
therefore  executing  at  t=0)  so  that,  at  any  instant  t,  we  have  trivially: 


T^t)  <  I  (t) 

VT>  i  VT) 


(4.3) 


Now,  since  in  the  firing  process  the  initiation  of  a  transition  •consumes* 
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one  token  from  its  input  places  and  the  termination  "produces*  one  token  in 


its  output  places,  we  have  trivially,  at  any  instant  x. 

for  the  transitions 

tA  and  tj  considered: 

hu)  <  ♦  Tj«> 

(4.4) 

Z)M  1  +  V° 

(4.5) 

Using  inequality  (4.3),  (4.4)  and  (4.5)  yields: 

V’* i  V'1 

(4.6) 

2  V'1 

(4.7) 

Consider  the  instant  of  time 


where  n  is  any  positive  integer  (this  is  the  instant  at  which  the 
(n+M^J-th  firing  initiation  of  transition  tA  occurs).  By  definition,  at 
this  instant: 


I 


i 


U)  -  °  *  Mji 


Replacing  in  (4.5)  yields: 

V'1 2  - 


Otherwise  stated,  the  number  of  initiations  of  transition  tj  in  [0,t]  is  at 
least  equal  to  n,  which  implies: 


Likewise,  considering  the  instant  of  time: 


n-Mjj 

*  =  Si  J 


where  n  is  assumed  greater  than  M^j,  then: 


Ii(t)  *  n  - 


and  replacing  in  (4.7)  yields: 


n  >  Ij(x) 


which  implies  this  time: 


Sn  >  t  =  s'1  MiJ 
Sj  “  T  Si 


Finally,  for  any  positive  integer  n,  such  that  n  >  M?j,  we  have: 


“-"Ij  <  n  <  “+"3l 

S1  -  Sj  *  S1 


Hence,  dividing  by  n: 
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ii'i 

[IV 


B 
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IV 
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•i  V : 


n  -  M, 


sn-M^ 

Si 


n  -  M, 


n+H  . 

S,  Ji  n  +  M„ 

<  _ _JA 

n  *  Mji  “ 


By  taking  the  limits  as  n  goes  to  infinity: 


n->+® 


n-H  n+M  . 

S.  ij  S.  S,  ji  n+M., 

-i—  <11.  -J- <  llm  -i—  - ^ 

n-M^j  n—>+“  n->+«  n+M^j  u 


which  means  exactly: 


e  <  e  <  e 

i  j  i 


therefore:  =  0j.  Q.E.D. 


In  the  next  sections,  we  will  denote  by  0  the  unique  average  cycle 


time  and 


f -i 


will  denote  the  average  firing  rate  of  all  the  transitions.  $ 
characterizes  the  average  processing  rate  of  the  overall  system  and  for 


that  reason,  will  be  called  the  throughput  rate. 


4.2.2.  Average  Circuit  Processing  Time 


In  the  previous  section,  we  have  proved  in  that  any  two  transitions  of 
a  directed  circuit  have  the  same  average  cycle  time.  It  is  now  possible  to 
determine  what  would  this  average  cycle  time  be,  assuming  that  the  circuit 
were  by  itself,  i.e.,  could  function  independently  of  the  remaining  system. 


» 


i 
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Let  us  consider  a  directed  circuit,  as  shown  in  Figure  4.2,  that  we 
denote,  without  loss  of  generality  by  p  =  (tx  pa  ta  ...  t^  p^).  Let  us 
assume  now  that  the  circuit  has  only  one  token  in  it,  which  is  initially  in 
place  p^.  Let  p(p)  denote  the  sum  of  the  transition  firing  times,  i.e., 

k 

p(p)  =  ^ 

i=l 


where  =  pC^)  is  the  firing  time  of  t^,  as  defined  in  Chapter  2,  Section 
2.4.1.  We  will  begin  by  assuming  that  the  firing  times  are  deterministic 
and  will  then  extend  the  analysis  to  the  random  case.  Now,  in  the  steady 
state  process,  the  token  fires  every  transition  in  the  circuit  in  turn  and 
reappears  in  p^  every  p(p)  seconds,  assuming  that  no  time  is  allowed  to 
elapse  between  a  transition  being  enabled  and  being  initiated.  Quite 
naturally,  we  shall  call  p(p)  the  circuit  processing  time  corresponding  to 
the  amount  of  time  it  takes  to  complete  the  processing  operations  of 
the  circuit,  under  the  previous  assumption.  In  that  case,  every 
transition  fires  at  intervals  of  p(p)  seconds:  therefore,  the  cycle  time  of 
any  transition  of  the  circuit  is  precisely  p(p),  i.e.,  6  -  p(p). 

tk 

< 

t2 

Figure  4.2  Directed  Circuit 


Now,  suppose  that  the  circuit  has  n  tokens  instead  of  1,  assuming  that 


all  tbe  tokens  are  Identical  (which  is  always  implicit  in  ordinary  Petri 
Nets).  Then,  the  average  cycle  time  of  every  transition  of  the  circuit 
becomes : 


6  =  . 
n 

Clearly  enough, 

p(p) 

n 

represents  the  average  amount  of  time  between  a  token  moving  from  p^  and  a 
token  reappearing  in  p^.  Extending  the  previous  analysis,  we  shall  call 
this  ratio  the  average  circuit  processing  time.  At  this  point,  it  is 
interesting  to  recall  that  the  token  content  of  a  circuit  represents,  in 
our  model,  the  resources  available  for  processing  the  oorrespoding 
operations  of  the  circuit.  Hence,  the  average  circuit  processing  time  is 
determined  as  the  sum  of  the  task  processing  times  divided  by  the  number  of 
resources  available.  Since  the  resources  (i.e.,  the  token  content)  are 
specified  by  the  initial  marking  of  the  places  of  the  circuit,  the  precise 
definition  is  finally  the  following  one: 

Definition:  Let  p  be  any  directed  circuit  that  we  denote  (without  loss  of 
generality)  by  p  =  (tx  pt  ...  tA  pA  ...  t^  pfc) .  The  average  circuit 
processing  time,  denoted  by  a(p),  is  defined  as: 


H±  =  H(ti>  denotes  the  firing  time  of  transition  tA 


denotes  the  initial  marking  of  place  p^,  and 

M°(p)  characterizes  the  token  content  of  the  circuit.  Furthermore, 
the  quantity 


1 

o(p) 

will  be  called  the  average  circuit  processing  rate. 

Now,  for  the  random  process,  corresponding  to  the  case  where  the 
transition  firing  times  may  be  discrete  random  variables  (with  a  finite  set 
of  values)  with  a  fixed  probability  distribution,  this  definition  i3  easily 
extended  by  taking  the  expected  firing  time  for  each  transition.  For 
instance,  let  us  suppose  that  the  firing  time  of  transition  tA  is  a  random 
variable  that  can  take  the  values 


‘■“i,.'  "l.. . •‘i.U1 


according  to  the  probability  distribution 


*Yi,i'  Yi,*' *  *  "Yi,u^ 


where  u  denotes  the  number  of  possible  processing  times  for  transition  tA. 
The  probability  distribution  satisfies,  as  usual,  the  equality 


J-u 


1=1 


(4.9) 


For  simplicity,  we  assume  that  the  processing  times  of  all  the  other 
transitions  but  t^  are  deterministic.  The  results  can  be  extended 
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trivially  to  the  case  where  more  than  one  transition  has  a  no  in¬ 
deterministic  processing  tine.  We  also  assume,  for  clarity,  that  the  token 
oontent  of  the  circuit  is  one.  Let  us  denote  then  by  the  expected 
firing  time  of  t^,  i.e.. 


j-u 

:i*  5  »i.j  “i.j 
J-l 


(4.10) 


It  turns  out  that  the  circuit  processing  time,  i.e.,  the  amount  of  time  it 
takes  for  the  token  to  complete  the  processing  operation  of  the  circuit, 
takes  the  value 

' *■  >■* 

whenever  the  firing  time  of  transition  t^  takes  the  value:  j. 
If  we  consider  a  large  number  of  repetitions  of  the  circuit,  i.e.,  assuming 
that  the  token  runs  cyclically  in  the  circuit,  the  expected  circuit 
processing  time  will  be: 


5(p>  -  }  <p,  *  »,  ♦•••*  Pi,j  ♦•••♦  **> 

j-i 

which  can  be  written  simply,  given  (4.9)  and  (4.10) 
p(p)  =  Pt  +  P,  +•••  +  +...+ 

Naturally,  in  the  general  case,  we  would  obtain 

p(p)  =  Px  +  P,  +...+  +...+  pk 


I 


assuming  that  all  the  transition  firing  times  are  discrete  random 
variables.  Likewise,  if  we  assume  that  there  are  n  tokens,  instead  of  one, 
given  the  assumption  that  the  probability  distribution  is  independent  from 
one  token  to  the  next,  as  stated  in  Chapter  3,  Section  3.4,  the  average 
circuit  processing  time  would  be: 

a(p)  =  (4.11) 

Hence,  in  the  non-determinlstic  case,  the  average  ciroult  processing  time, 
as  defined  by  (4.11)  determines  the  average  cycle  time  of  all  the 
transitions.  Otherwise  stated,  (4.8)  is  simply  extended  by  taking,  in  that 
case,  the  expected  firing  time  of  each  transition. 

In  brief,  we  have  showed  that,  under  the  assumption  that  the  circuit 
is  by  itself,  and  that  no  time  would  elapse  between  a  transition  being 
enabled  and  being  initiated,  the  average  circuit  processing  time,  as 
defined  previously,  determines  the  average  cycle  time  of  every  transition 
of  the  circuit. 

4.2.3  Maximum  Average  Circuit  Processing  Time 

In  the  previous  section,  we  have  considered  a  circuit  of  the  net,  as 
if  it  were  Isolated  from  the  remaining  system  and  we  have  determined  in 
that  case  what  would  be  the  average  cycle  time  of  each  transition  of  the 
circuit.  The  fundamental  assumption  was  that  there  was  no  time  delays 
between  the  different  processing  operations.  Considering  now  the  system  as 
a  whole,  it  is  clear  that  the  different  circuits  are  in  fact 
interconnected,  as  shown  for  instance  in  Figure  4.3. 

It  turns  out  that  the  Interconnected  circuits  will  affect  each  others 
processing  time.  For  instance,  in  the  example  of  Figure  4.3,  transition  tx 
belongs  to  both  circuit  px  and  p,.  Clearly  enough,  the  average  cycle  time 
of  this  transition  cannot  to  lower  than  the  maximum  of  a(px)  and  a(px). 


O  “  11111  ^a(pi)  *****  a(pp) j 

will  be  called  the  minimum  average  circuit  processing  rate.  At  this  point, 
it  is  emphasized  that  the  minimum  is  taken  over  all  the  direoted  circuits 
of  the  net. 

We  have  now  the  following  result: 

Theorem  4.2:  The  maximum  average  circuit  processing  time  is  a  lower  bound 
of  the  average  cycle  time,  i.e.,  B  1  a. 

Let  us  first  point  out  that  this  result  holds  only  for  the  type  of 
Petri  Nets  used  here  to  modeling  the  DHO,  i.e.,  strongly  connected  Event- 
Graphs.  As  already  stated  in  the  introduction  of  this  chapter,  this  result 
was  first  obtained  by  Ranchandami  [6]  but  with  restrictive  assumptions:  In 
his  work,  the  transition  firing  times  were  assumed  deterministic  and, 
moreover,  the  steady-state  process  was  assumed  strongly  periodic  (i.e., 
transitions  fire  at  regular  intervals  of  time).  We  extend  the  result  to 
the  non-deterministic  case,  where  firing  times  may  be  discrete  random 
variables.  In  addition,  we  do  not  make  any  assumption  regarding  the 
periodicity  of  the  deterministic  system,  which  will  be  shown  in  the  next 
chapter  to  be  of  a  more  complex  type  (as  already  mentioned  before,  the 
deterministic  process  is  in  fact  K-periodic). 

In  order  to  prove  Theorem  4.2,  it  is  sufficient  to  prove  that,  for  any 
directed  circuit  p  of  the  net,  we  have  necessarily:  B  1  a  (p).  This  is 
actually  the  proof  given  below. 


Proof:  Let  p  be  any  directed  circuit  that  we  denote  (without  loss  of 

generality)  p  =  (t4  px  ta  ...  tk  pk).  Let  us  consider,  for  instance,  the 
sequence  S*  (n  *  1,2,3,...)  corresponding  to  the  instant  at  which  the 
firings  of  transition  tt  are  initiated  successively.  We  first  assume  that 
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all  the  transition  firing  times  are  deterministic 
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Now,  let  us  consider  any  couple  of  transitions  (t^,  ti+l)  of  the 
clrouit  that  are  connected  by  place  p^  as  shown  in  Figure  4.2  (p^  is  the 
output  place  of  t^  and  the  input  place  of  t^+1  in  the  clrouit).  The 
inequalities  (4.4)  -  (4.5)  established  previously  for  (t^,tj)  can  be 
applied  here  for  (t^»ti+1),  so  that,  at  any  instant  x: 


W0  -  <  *  v*> 


(4.13) 


where  denotes  the  initial  marking  (at  x=0)  of  place  pA- 


Let  us  consider  then  the  instant: 


x  =  S^“M1  +  m1 


where  n  is  any  positive  integer,  such  that  n  >  M*  and  ^  is  the  firing  time 
of  tA .  x  is  the  instant  of  the  (n-M^)-th  firing  termination  of  transition 
t^,  hence  at  this  instant: 


T±(x)  =  n  -  M 


Using  relation  (4.13)  yields: 


Ii+»(x)  -  a 


which  implies: 


(4.14) 


I* 


i 


Now,  by  iterating  (4.14)  from  i=l  to  i-k,  we  obtain: 


S 


n 

1 


>  gn-M#(p) 


k 


i=l 


where 


k 

M°(p)  =  }  M® 
i=l 


Finally,  for  any  positive  integer  n,  such  that  n  >  M°(p): 

gn  >  gn-M  (p)  +  p(p)  (4.15) 

where  M#(p)  denotes  the  token  content  of  the  circuit  and  p(p)  the  circuit 
processing  time. 

Let  us  now  take  a  place  of  the  circuit  which  has  a  non-null  initial 
marking  (i.e.,  contains  initially  at  least  one  token).  We  will  assume, 
without  loss  in  generality,  that  this  is  the  case  for  place  pt.  If  we 
consider  px  as  the  origin  for  the  circuit,  we  shall  say  that  the  oirouit 
has  been  repeated  exactly  K  times  (K=l,2,3...)  when  all  the  tokens  have 
gone  through  px  exactly  K  times  (i.e.,  all  the  tokens  have  been  prooessed  K 
times  through  the  circuit).  In  that  case,  the  number,  n,  of  firing 
occurences  of  transition  t4  is  determined  by: 
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n  =  K  M  (p)  +  r 


0  <  r  <  M*  (p) 


where  r  denotes  the  number  of  tokens  that  are  initially  not  in  place  px. 
Indeed,  all  these  tokens  must  be  fired  by  transition  tx  once  to  reach  place 
px  (the  origin)  for  the  first  time. 


By  iterating  (4.15)  for  n  *  r,  H*(p)  +  r  ,...,  K  M#(p)  +  r,  we  obtain: 


s“  2  s[  +  K  p(p) 


(4.16) 


Now,  dividing  by  n,  we  obtain: 


n  K  M#(p)  +  r  K  M*(p)  +  r 


Finally,  when  n  goes  to  infinity: 


SQ  Sr  _ 

lim  —  2  lim  - - — - -  +  lim  - - - 

n-)+<D  n  K->+®  KM  (p)  +  r  K-*+®  KM  (p)  +  r 


=  nki 


M  (p) 


which  means,  given  (4.2)  and  (4.8): 


6  =  lim  —  2  o(p)  =  ^  f  P--- 


n->  +®  n  M# ( p ) 


Q.E.D. 


In  case  where  the  firing  times  are  discrete  random  variables,  it  is 


uwmmamam®*** 
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not  difficult  to  extend  the  proof.  Suppose,  for  instance,  that  the  firing 
time  of  t^  is  a  random  variable  that  can  take  the  values 


**!,»**  *  * '^i.u3 


according  to  the  probability  distribution 


tYi,x'  Yi.a'"*'ri.u3 


In  addition,  all  the  other  processing  times  are  deterministic,  exactly  as 
we  have  assumed  in  the  previous  section.  Recall  that  the  circuit 
processing  time  takes  the  value 


^(P)  l-i.J  I*, 


with  probability  when  pi( t± )  =  p^j.  Clearly  enough,  (4.16)  becomes 


j-1 


(4.17) 


where  Kj  denotes  the  number  of  times  that  the  circuit  processing  time  was 
equal  to  ni  j (p )  for  the  K  occurences  of  the  circuit.  By  the  definition  of 
the  probability  distribution: 


lim 
K— )  +°° 


*1 

•=**■  =  y,  ,  for  j  =  1,2,. ...u. 

2.  t  J 


From  (4.17),  we  deduce: 


warn* 


When  □  goes  to  infinity,  we  obtain  the  same  result  given  in  (4.11): 


0  =  lim 
n->°° 


Ihl. 

M#(p) 
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M  (p) 


a(p) 


Remark:  In  the  determination  of  the  maximum  average  circuit  processing 

time  (as  defined  by  (4.12)),  it  is  necessary  to  include  all  the  transition 
processing  times.  Indeed,  in  our  model  of  the  DMO,  each  transition 
contains  in  fact  a  self-loop  with  one  token,  as  discussed  in  Chapter  2, 
Section  2.4.2.  Therefore,  each  transition  constitutes  by  itself  a  circuit, 
whose  token  content  is  one.  In  that  case,  the  average  circuit  processing 
time  is  trivially  the  processing  time  of  the  transition.  In  order  to 
differentiate  later  the  self-loops  from  the  other  directed  circuits,  we 
will  refer  to  the  former  as  the  trivial  circuits. 


4.2.4  Maximum  Throughput  Rate 

In  the  previous  section,  we  have  found  that  the  maximum  average 
circuit  processing  time,  as  determined  by  (4.12),  is  a  lower  bound  of  the 
average  cycle  time,  0.  Naturally,  this  result  can  also  be  stated  this  way: 

The  throughput  rate  of  the  system,  f  =1/0,  as  defined  in  Section 
4.2.1,  is  bounded  from  above  by  the  minimum  average  circuit  processing 
rate,  i.e.,  by  1/a 

The  problem  before  us  now  is  to  know  whether  or  not  this  bound  is 
actually  achievable.  In  other  words,  we  have  to  determine  whether  the 
minimum  average  circuit  processing  rate  characterizes  the  maximum 


throughput  rate  of  the  system.  It  is  clear  that  if  the  arrival  rate  of 
external  inputs  is  low  enough,  the  DMO  will  be  able  to  handle  all  the 
arriving  inputs  and  the  throughput  rate  will  be  precisely  the  rate  at  which 
inputs  arrive.  However,  there  is  a  rate  above  which  the  DMO  will  be 
overloaded,  in  the  sense  that  inputs  will  queue  at  the  entry  of  the  system 
and  this  queue  will  grow  to  infinity  over  time.  This  rate  characterizes 
the  maximum  throughput  rate  of  the  DMO,  i.e.,  the  maximum  rate  at  which 
inputs  can  be  processed  in  the  organization.  Prom  the  analysis  carried  out 
so  far,  we  know  that  the  maximum  throughput  rate  is  necessarily  lower  or 
equal  to  the  minimum  average  circuit  processing  rate,  as  determine  by  1/a. 
Indeed,  the  average  rate  at  which  transitions  can  fire,  i.e.,  0,  cannot 
exceed  1/a  (Theorem  4.2):  Therefore,  if  the  arrival  rate  of  inputs  exceeds 
1/a,  we  are  sure  that  the  queue  of  inputs  will  grow  to  infinity  over  time. 
As  we  are  going  to  see,  1/a  characterizes  in  fact  the  maximum  throughput 
rate  for  the  deterministic  system,  i.e.,  when  all  the  transition  firing 
times  are  deterministic.  It  will  be  shown  that  it  is  generally  not  true 


for  the  non-deterministic  case. 

4. 2. 4.1  Deterministic  System 

For  the  deterministic  system,  the  following  result  holds: 

Theorem  4.3:  The  minimum  average  circuit  processing  rate  determines  the 
maximum  throughput  rate. 

As  stated  earlier,  we  have  proved  that  the  minimum  average  circuit 
processing  rate,  as  determined  by  1/a  (relation  (4.12)),  bounds  the  maximum 
throughput  rate.  We  will  not  prove  here  that  1/a  corresponds  effectively 
to  the  maximum  throughput  rate.  This  result  will  be  derived  from  the 
analysis  carried  out  in  Chapter  5,  which  concerns  the  execution  schedule  of 
the  system.  In  Chapter  5,  we  will  determine  the  precise  instants  of  time 
at  which  transitions  fire,  assuming  that  the  process  occurs  repetitively 
and  that  the  firings  are  initiated  as  soon  as  the  transitions  are  enabled. 
From  the  execution  schedule  obtained,  it  will  be  shown  that  the  (average) 


84 


rate  at  which  the  transition  fires  is  determined  by  the  minimum  average 
circuit  processing  rate. 


Theorem  4.3  allows  the  maximum  throughput  rate  to  be  computed  very 
easily,  once  we  have  determined  all  the  directed  circuits  of  the  net  (which 
only  depend  on  the  structure  of  the  Petri  Net).  In  order  to  determine  all 
the  directed  circuits,  we  have  used  an  algorithm  developed  by  Alalwan  and 
Toudlc  [19]  that  determines  the  S-lnvarlants  (Chapter  2,  Section  2.2)  of  an 
ordinary  Petri  Net.  This  algorithm  is  described  in  Appendix  A  of  the 
thesis. 

Let  us  analyze  further  how  we  can  relate  the  evaluation  of  the  maximum 
throughput  rate  to  the  resource  and  time  constraints.  To  that  extent,  we 
need  a  further  definition  characterizing  the  critical  cirouits. 

Definition:  A  circuit  p^  is  said  to  be  critical,  if  its  average  processing 
time  is  maximal  over  all  the  directed  circuits,  i.e., 

a  =  max  (a(pi)  ,...,  aip^))  =  o(p^) 

The  following  corollary  of  Theorem  4.3  is  straightforward: 

Corollary  4.3:  The  average  processing  rate  of  the  critical  circuits 
determines  the  maximum  throughput  rate  (for  the  deterministic  system). 

This  corollary  is  very  interesting  because  it  provides  a  way  for 
determining  which  resource  and  time  constraints  are  actually  binding,  i.e., 
bound  the  throughput  rate.  Indeed,  these  will  be  the  task  processing  times 
(resp.  resources)  that  correspond  to  the  transitions  (resp.  token  content) 
of  the  crl tical  circuits.  Hence,  the  problem  of  modifying  the  right 
constraints  so  as  to  improve  the  maximum  throughput  rate  beoomes 
straightforward.  In  addition,  by  determining  the  critical  circuits  of 
different  organizational  structures,  a  comparative  study  of  performance  can 


be  achieved,  with  respect  to  the  throughput  rate.  This  type  of  performance 
evaluation  will  be  carried  out  In  the  last  section  of  this  chapter  for  the 
example  of  the  two  DM  organization. 

4. 2. 4. 2  Non-Deterministic  System 

Let  us  recall  that  the  characteristic  of  the  non-determlnlstlc  system 
studied  here  is  the  following:  The  transition  firing  times  are  discrete 
random  variables,  that  can  take  a  finite  number  of  possible  values, 
according  to  a  fixed  probability  distribution  (i.e.,  independent  from  one 
token  to  the  next).  As  discussed  in  Chapter  3,  Section  3.4,  these 
assumptions  allow  for  modeling  the  decision  nodes  or  switches .  Recall  also 
that  the  average  circuit  processing  time  is  computed  by  taking  the  expected 
(or  mean)  firing  time  of  each  transition.  In  that  case,  it  turns  out  that, 
in  general,  the  maximum  throughput  rate  is  lower  than  the  minimum  average 
circuit  processing  rate,  as  defined  by  relation  (4.12).  Let  us  see  why  on 
a  simple  example. 

Let  us  consider  the  Petri  Net  shown  on  Figure  4.4,  which  is  a  strongly 
connected  Event-Graph  with  the  following  two  circuits: 

Pa  -  Pt  *a  pa) 

Pa  =  P,  ta  P4> 

Now,  suppose  that: 

-  the  token  content  of  both  circuit  is  one 

-  the  firing  time  of  tB  is  deterministic:  |i#=l. 

-  the  set  of  possible  firing  times  for  transition  tx  is  =  {1,3} 
according  to  the  probability  distribution:  {y1>a  =  0.5,  ra>*  =  0.5} 
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-  the  set  of  possible  firing  times  for  ta  is  pa=  {2,4},  according  to 

the  probability  distribution:  {T*fl.  ■  0.5,  =  0.5}. 

There  are  four  deterministic  cases  that  correspond  to  the  randan 
system.  In  each  case,  we  can  easily  compute  the  best  cycle  time 
achievable,  as  determined  by  the  maximum  average  circuit  processing  time 
(Theorem  2.3).  These  are: 

-  for  pa  =  1,  pa  =  2  then  ax  =  a(pa)  *  pa  +  p#  *  3 

-  for  pa  -  1,  pa  -  4  then  aa  =  a(pa)  =  Pa  +  P,  •  5 

-  for  px  =  3,  pa  =  2  then  aa  =  a(px)  =  px  +  p0  =  4 

-  for  px  =  3,  pa  =  4  then  a4  *  a(pa)  =  pa  +  p,  *  5 

Given  the  probability  distribution,  the  average  cycle  time  of  the  system 
(taken  over  a  large  number  of  repetitions  of  the  process)  is  clearly: 

0  15  +  +  Tl,***,!®*  +  Yx,»Ta,xa4 

=  0.25  (3+5+4+5)  -=  4.25 

As  said  earlier,  the  maximum  throughput  rate  corresponds  to  the  average 
rate  at  which  transitions  fire,  assuming  that  the  system  works  oontlnously 
and  that  the  firings  are  initiated  as  soon  as  the  transitions  are  enabled. 


By  construction,  f  =  1/4.25  determines  the  maximum  throughput  rate  in  that 
case.  If  we  compute  now  the  maximum  average  circuit  processing  time  of  the 
random  system,  by  taking  the  expected  firing  time  of  tx  and  ts,  we  obtain: 


'i 


otpj)  =  jlj  +  p#  =  0. 5(1+3)  +1*3 
o(p,)  *  }ia  +  p0  =  0. 5(2+4)  +  1  =  4 
a  =  max  (a(px),  o(p,))  =  4 


and  therefore  the  minimum  average  circuit  processing  rate  equals: 


1  1 

_  s  — 

a  4 


Finally,  it  turns  out  that: 


i  = 


l 

4.25 


1. 
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This  example  proves  that,  in  oontrast  to  the  deterministic  system,  the 
maximum  throughput  rate  can  be  lower  then  the  minimum  average  circuit 
processing  rate,  as  determined  by  1/a.  The  reason  should  be  clear  from 
this  example:  the  critical  circuits  are  not  necessarily  the  same  among  the 
deterministic  systems  obtained  by  setting  a  deterministic  value  (among 
those  possible)  for  each  transition  firing  time.  In  our  example,  when  px=l 
and  p,=2,  pa  is  the  critical  cirouit  and  when  px  =  3,  ps  =  2,  the  critical 
circuit  is  px.  By  averaging  the  corresponding  circuit  processing  times 
over  the  probability  distribution,  we  do  not  necessarily  obtain  a(p),  where 
p  is  the  critical  circuit  obtained  when  the  different  cirouit  processing 
times  are  computed  by  taking  the  expected  firing  times  of  each  transition. 


In  the  non-determinlstlc  case,  we  can  obtain  a  lower  and  an  upper 
bound  for  tbe  maximum  throughput  rate.  From  tbe  analysis  above,  tbe  upper 
bound  is  given  by  1/a,  tbe  minimum  average  circuit  processing  rate.  If  tbe 
arrival  rate  of  inputs  exceeds  1/a,  we  are  sure  that  tbe  DM0  will  be 
overloaded,  i.e.,  tbe  queue  of  inputs  will  grow  to  infinity  over  time. 
Now,  tbe  lower  bound  is  trivially  obtained  by  assigning  to  each  transition 
its  worst,  i.e.,  longest,  firing  time  possible  and  computing  tbe  minimum 
average  circuit  processing  rate  for  tbe  corresponding  deterministic  system. 
For  arrival  rates  of  inputs  that  do  not  exceed  this  lower  bound,  we  are 
sure  that  tbe  DM0  will  be  able  to  handle  all  inputs.  Interestingly  enough, 
if  tbe  upper  and  lower  bound  turned  out  to  be  tbe  same,  whicb  is  tbe  case, 
when  in  our  model  of  the  DM0  the  critical  circuits  do  not  contain  switches, 
we  bave  tben  completely  characterized  tbe  maximum  throughput  rate.  Let  us 
point  out  that,  as  for  tbe  deterministic  system,  tbe  lower  and  upper  bound 
so  determined  are  directly  expressed  as  a  function  of  both  the  time  and 
resource  constraints.  Therefore,  the  performance  evaluation,  based  on  the 
circuits  that  are  critical  can  also  be  applied  here. 

4.3  APPLICATION  TO  THE  TWO  MEMBER  ORGANIZATION 

In  this  section,  we  illustrate  the  results  obtained  so  far  using  the 
example  of  the  two  member  organization  presented  in  Chapter  3,  Section  3.5. 
We  recall  that  the  first  structure,  sbown  on  Fig.  3.7,  includes  an 
interaction  of  the  result  sharing  type  between  the  two  DMs,  while  the 
second  one,  shown  on  Fig.  3.8,  implies  a  hierarchical  relationship  between 
them. 

4.3.1  Analysis  of  the  Processing  Rate  Constraints 


In  order  to  compute  the  maximum  throughput  rate,  it  is  first  necessary 
to  determine  all  the  directed  circuits  of  the  net.  The  algorithm  developed 
by  Alawain  and  Toudic  [19],  that  is  described  in  Appendix  A,  has  been  used. 
In  this  simple  example,  it  is  however  easy  to  check  the  directed  circuits 
graphically.  Given  the  labeling  of  transitions  and  places  shown  in  Figs. 


3.7  and  3.8,  there  are  three  directed  circuits  that  have  the  same  structure 
in  both  cases.  These  are: 

Pi  “  P*  t,  p,  t4  p4  t,  Rj 

P*  *  P «  t ,  pT  t,  p,  t,  R») 

P*  *  ^i  P t  P«  *7  Pi  P,  t,  p14  t1#  R#) 

Additionally,  for  structure  (a).  Fig.  3.7, 

P*a  ~  P»  *■»  P*  ^4  P<  P xo  P7  P» 

P*a  ~  ^1  Pi  ^1  P*  Pi  ^4  P4  Pi#  ^7  P7  P#  P11  t*i#  R# ) 

P#a  =  *fci  p*  tt  P»  fc*  t4  P4  tf  Pi«  fc7  P7  t#  P#  fc»  Pll  fci#  *#> 

and  for  structure  (b),  Fig.  3.8, 

P 4 b  "  p,  t,  p,  t4  p4  t,  p1#  t,  p,  t,  R4) 

P*b  -  ^1  Pi  ^a  P»  P»  ^4  P4  P10  P#  P11  ^10  R#) 

P#b  =  (ti  P*  P#  t,  p,  t4  p4  tf  plB  t,  p,  t#  plx  tl0  R# ) 

If  we  denote  by  ^  the  processing  time  of  transition  tA  and  by  n#  (resp. 

nx,  na)  the  amount  of  resources  (l.e.,  the  number  of  tokens  that  is 
initially  in  the  corresponding  resource  place)  available  for  processing  the 
organizational  inputs  (resp.  the  inputs  to  DM1  and  DM2),  then  the  average 
circuit  processing  times  are: 


We  recall  that,  in  order  to  determine  the  maximum  throughput  rate.  It  Is 
also  necessary  to  take  into  account  the  processing  times  of  the  individual 
transitions,  which  constitute  the  trivial  circuits  of  the  net.  Let  us 


max 


u  =  max  11. 

*“  1-1 . 10  1 


Then,  the  maximum  average  circuit  processing  time  is,  for  structure  (a): 


a  =  max  (aCpj),  a(p,),  a(pB),  a(p4_),  o(p,  ),  a(p4  ),  p  ) 


and,  for  structure  (b): 


ab  =  max  (a(p1),  a(p4),  a(pt),  o(p4b),  a(pfb),  a(p4b),  nmax) 


However,  we  have  seen  in  Chapter  3,  Section  3.2,  (relation  (3.5))  that: 


/  1  1  1  1  1  \ 
1  (a(Pl)  ’  H,  '  M,'  *,) 


describes  the  maximum  processing  rate  that  characterizes  DM1.  Identically, 
for  DM2,  it  is: 


=  min 


/— 1  i_  i_\ 

\a<P.)'  •».  »*.’  S/ 


Moreover,  it  turns  out  that: 


+  M4  +  M, 

a(P«a)  =  o(p3)  + - n - 


Hence  a(p4a)  2  a(pt),  whatever  the  processing  times  and  the  resources  may 
be.  Similarly,  we  have  always: 
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°<P<b>  -  a(p4) 


We  deduce  finally  that  the  minimum  average  circuit  processing  rate  is 
determined  for  structure  (a),  by: 


i 


a 


=  min 


'*■  ““-.a* 


a(fisa) 


1  1 
a(p«a>  '  *** 


and  for  structure  (b): 


* 


b 


min 


o(P.b) 


o(P«b) 


-L.\ 

\of 


Now,  the  expressions  above  are  especially  interesting  to  interpret.  It 
turns  out  that  the  throughput  rate  of  the  DMO  is  bounded  by  three  types  of 
constraints: 

(i)  The  constraint  that  comes  from  the  maximum  processing  rate  of 
each  DM,  as  defined  by  px  and  p%. 

(ii)  The  constraints  due  to  the  structure  of  the  organization,  as 
determined  by  the  interactions  between  the  DMs  and  between  the 
DMs  and  the  environment.  In  the  example,  these  constraints  are 
derived  from  the  average  processing  rates  of  the  circuits  p4a, 
p<a  and  p<a  (resp.  p4b,  plb  and  p§b).  We  will  refer  to  these 
constraints  as  structural  constraints. 

(iii)  The  time  constraint  px,(resp.  p1#),  specific  to  the  first  (resp. 
last)  processing  stage,  which  corresponds  to  the  allocation  of 
Inputs  between  the  DMs  (resp.  the  production  of  the 
organizational  response).  These  tasks  are  Indeed  different  from 
the  four  processing  stages  that  characterize  the  DM,  i.e.,  the 
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SA,  IF,  Cl,  and  RS  stages,  as  described  in  Chapter  3,  Section 
3.1. 


However,  it  should  be  clear  that  all  these  constraints  are  actually 
strongly  Interrelated;  For  example,  the  resources  available  to  DM2  (l.e., 
na)  are  used  to  compute  pa  but  also  to  compute  the  average  processing  rate 
of  the  circuit  p4a.  Likewise,  the  average  processing  rate  of  the  cirouits 
pJa  and  p4a  depends  on  p4  and  p10.  In  particular,  the  critical  cirouits 
which,  as  said  earlier,  determine  the  maximum  throughput  rate,  may  be  at 
the  same  time  in  the  three  categories  described  above.  In  other  words, 
there  is  not  necessarily  only  one  type  of  constraint  which  is  the  binding 
constraint,  given  the  values  of  the  various  processing  times  and  resources 
available. 

The  differentiation  is  nevertheless  very  useful  in  carrying  out  the 
performance  evaluation,  especially  when  the  point  is  to  compare  different 
protocols  of  operations  for  the  same  organization.  Indeed,  in  that  case, 
it  is  sufficient  to  look  at  the  structural  constraint  only.  For  Instance, 
assuming  that  the  processing  times  and  resources  available  are  identical 
for  both  organizations  (structure  (a)  and  <b) ) ,  it  turns  out  that  the 
following  inequalities  are  always  verified: 

a(p4b)  <  a(p4a) 

a(pJb)  <  a(p,a) 

a(p4b)  <  a(p4a) 


Therefore,  §b  2  ia. 


In  fact,  the  equality  between  #b  and  fa  is  only  verified  if  the 
binding  constraint  comes  from  the  processing  rate  of  one  of  the  DMs 
(i.e.,  fa  =  pa  or  fa  =  pa)  or  from  the  processing  time  of  the  input  or 


output  task,  i.e.: 


In  any  case,  it  turns  out  that  the  performance  of  the  hierarchical 
structure  (structure  (b))  is  always  at  least  as  good  as  that  of  the  result 
sharing  type  (structure  (a)).  We  are  now  going  to  analyze  further  how  the 
performance  of  the  organization  is  related  to  the  different  types  of 
constraints  previously  described,  by  assuming  various  numerical  values  for 
the  resource  and  time  constraints.  To  that  extent,  we  will  consider  only 
one  organizational  form,  structure  (b)  of  the  two  member  organization. 


4.3.2  Evaluation  of  the  Maximum  Throu 


jut  Rate 


In  this  section,  we  consider  the  two  member  organization,  with 
structure  (b),  as  shown  in  Figure  3.8.  For  simplicity,  all  the  processing 
times  are  assumed  deterministic.  Indeed,  our  objective  here  is  to  make 

clear  how  to  carry  out  the  performance  evaluation  depending  on  the  type  of 

binding  constraints.  To  that  extent,  and  according  to  the  development 
presented  in  the  previous  section,  we  will  consider  three  cases:  The  first 
one  is  an  example  where  the  throughput  rate  is  bounded  by  the  maximum 

processing  rate  of  one  DM.  In  the  second  one,  the  binding  processing  rate 

constraint  comes  from  a  structural  constraint.  In  the  last  one,  the 
processing  time  of  the  partitioning  operation  (which  corresponds  to  px)  is 
the  binding  constraint. 


Remark:  For  clarity,  the  circuits  previously  denoted  by  p4b,  p<b  and  p4b 
will  be  simply  denoted  by  p4,  pf,  and  p4. 

4. 3. 2.1  Example  1 

Let  us  assume  that  it  takes  one  unit  of  time  to  complete  any  one  of 


the  tasks,  except  for  the  SA  stage  of  DM1,  which  takes  four  units  of  time, 
i.e.,  pj;  =  1  for  i  =  1,3, 4,..., 10  and  p,  *  4.  In  addition,  DM1  and  DM2  can 
only  process  one  input  at  a  time,  i.e., 

n»  =  n4  =  1 

and  the  overall  organization  can  handle  at  most  eight  inputs  at  a  time, 
i.e.. 


n0  =  8 


In  that  case,  the  average  processing  times  of  the  different  circuits  are: 


a(pi) 
a(pa ) 

a(p,) 

a(p4) 

a(p,) 

a(Pt) 


3_ 

4 


=  7 

=  4 

=  £ 

8  = 

=  6 

=  _U 
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and 


max 


max 


i=li 


,10 


“i.  =  “a  =  4 


The  maximum  throughput  rate  is  therefore  determined  by: 


*  ,  /I  1  4  1  8  .  .  1 

f  min  (?  ,  4  ,  3  .  6  ,  n  ,  U  -  7 


The  maximum  processing  rate  of  each  DM  is: 
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min  ( 


o(px) 
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min 


a(p,) 
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The  maximum  processing  rate  of  DM1  is  clearly  the  binding  constraint  to  the 
throughput  rate  of  the  overall  organization.  Indeed,  pa  is  the  unique 
critical  circuit  of  the  net. 


Now,  there  are  only  two  possible  alternatives  so  as  to  improve  the 
performance  of  the  organization:  either  increase  the  resources  available  to 
DM1  or  reduce  the  processing  time  of  one  its  processing  stages  (or  both), 
which  are  indeed  the  only  ways  to  reduce  a(px).  For  Instance,  it  would  be 
possible  to  improve  the  performance,  if  DM1  was  better  trained,  so  that  he 
could  handle  more  information  at  the  same  time  (which  means  that  he  has 
more  cognitive  resources)  or  execute  the  tasks  faster  (or  both).  In  any 
case,  the  only  way  to  improve  the  organizational  performance  is  by 
improving  the  individual  performance  of  DM1 . 


4. 3. 2. 2  Example  2 


Let  us  keep  the  same  values  as  before,  except  that  pa  equals  now  one 
unit  of  time  like  the  other  task  processing  times.  In  this  example: 

3 

a(pt)  =  a(p#)  =  4,  o(pt)  =  j  ,  a(p4)  =  6  ,  a(pf)  =  o(p4)  =  1 

p  =1 

max 


Hence : 


& 


-•vv 


,  <  ,1  4  1  ,  .  1 

{  *  min  {4  *  3  »  g  »  “  6 


and  the  maximum  processing  rate  of  each  DM  is  the  same: 


K  -  -  T 


The  maximum  throughput  rate  is  determined  by  a  structural  constraint,  sinoe 
p4  is  now  the  unique  critical  circuit.  Given  the  structure  of  p4.  the  only 
way  to  improve  the  throughput  rate  is,  in  that  case,  by  either  increasing 
the  resources  available  to  DM2  (i .e. ,  na)  or  reducing  the  processing  times 
of  the  corresponding  tasks  of  the  circuits  (or  both).  Interestingly 
enough,  it  will  have  no  effect  to  either  increase  the  resources  of  DM1 
(i.e.,  nx)  or  the  processing  time  of  its  SA  stage  (i.e.,  pa),  in  oontrast 
to  the  previous  case. 


4. 3. 2. 3  Example  3 


At  last,  let  us  assume  that: 


-ns-: 


p±  =  1  for  i  =  2, 


,10 


px  -  2  (processing  time  of  the  input  transition,  corresponding 
to  the  partitioning  operation) 

n,  «  8,  nx  =  n4  =  4  (Both  DMs  can  now  handle  four  inputs 

simultaneously) 


m 


then,  in  that  case: 


il] 


y.ij 

’*1 


$1 
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t(px)  ■  o(pt)  ■  1  ,  a(p4)  “  y  ,  o(p4)  -  y  ,  (p4)  -  a(p4)  “  y 


-  2 


Hence: 


.  4  ,,  *  2  8  1,  _  1 
x  *  nln  *  ij  x  ^  i  2  2 


In  that  case,  the  processing  time  of  t1,  corresponding  to  the 
partitioning  operation,  determines  the  maximum  throughput  rate.  The  way  to 
improve  the  performance  of  the  organization  is  very  simple  in  that  case, 
sinoe  the  only  alternative  possible  is  to  reduoe  the  amount  of  time  it 
takes  to  allocate  the  external  inputs  to  the  DMs  interacting  directly  with 


the  environment. 


4 .4  CONCLUSION 


We  have  showed  in  this  chapter  that  the  maximum  throughput  rate  is 
bounded  from  above  by  the  minimum  average  circuit  processing  rate.  In 
addition,  this  bound  characterizes  the  maximum  throughput  rate  for 
deterministic  systems.  In  the  next  chapter,  we  will  analyze  the  execution 
schedule,  i.e.,  the  instant  of  times  at  which  the  various  operations  take 


place  in  the  process. 


CHAPTER  T 


ANALYSIS  OF  EXECUTION  SCHEDULES 


In  the  previous  chapter,  we  have  carried  out  the  computation  of  the 
maximum  throughput  rate  as  a  function  of  both  the  time  and  resource 
constraints,  specific  to  the  DMO.  This  throughput  rate  bounds  the  rate 
with  which  inputs  can  arrive  without  overloading  the  organization.  In  this 
chapter,  we  want  to  determine  the  execution  schedule  (also  oalled  firing 
schedule,  by  reference  to  the  transition  firings),  l.e.,  the  Instant  of 
times  at  which  the  various  operations  occur  in  the  prooess,  for  allowable 
rates  of  incoming  Inputs.  It  is  emphasized  that  the  firing  schedule 
computed  here  will  characterize  the  dynamic  behavior  of  the  DMO:  starting 
from  an  Initial  state,  the  prooess  will  be  assumed  to  ooour  repetitively, 
which  supposes  that  there  are  always  inputs  available  at  the  entry  of  the 
system.  By  doing  so,  we  will  obtain  the  best  performance,  with  respeot  to 
real-time  processing,  that  can  be  achieved  by  the  organization. 

Naturally,  the  precise  firing  schedule  oan  only  be  obtained  for  the 
deterministic  process,  in  which  all  the  processing  times  are  deterministic. 
For  the  non-deterministic  case,  where  processing  times  may  be  discrete 
random  variables  (as  defined  for  switches),  it  will  be  possible  to  obtain 
the  earliest  (resp.  latest)  firing  sohedule  by  taking  the  shortest  (resp. 
longest)  firing  time  of  each  transition. 

5.1  ANALYSIS  OF  SEQUENTIAL  AND  CONCURRENT  OPERATIONS 

Before  modeling  the  dynamic  behavior  of  the  DMO,  taking  into  aooount 
the  time  and  resource  constraints,  it  is  first  neoessary  to  analyze  further 
bow  the  different  operations  take  place  in  the  prooess.  Indeed,  the 
precise  sequence  of  operations  in  the  prooess  is  completely  determined  by 
the  structure  of  the  net.  In  other  words,  the  structure  of  the  Petri  Net 
characterizes  the  causality  of  the  system,  i.e.,  the  faot  that  certain 


transitions  oust  fire  before  others,  corresponding  to  the  sequential  tasks, 
while  others  can  fire  independently,  corresponding  to  the  concurrent  tasks. 


Apart  from  the  task  processing  times  and  the  resources  available,  this 
partial  order  in  the  different  operations,  as  determined  by  the  structure 
of  the  net,  is  clearly  a  determinant  factor  of  the  execution  schedule. 
This  is  the  reason  why  our  objective  here  is  to  find  a  convenient 
representation  that  could  account  for  the  causal  dependence  and 
independence  between  the  various  processing  tasks. 

5.1.1  Representation  of  the  Process  as  an  Occurence  Net 

In  order  to  model  the  precise  sequence  of  activities  in  the  process, 
we  first  consider  the  net  that  is  obtained  by  deleting  all  the  resource 
places,  including  the  resource  place  of  the  overall  organization  (R, )  and 
of  each  DM.  The  corresponding  net  for  the  example  of  the  two  member 
organization  described  in  Chapter  3,  Section  3.5  (structure  (a))  is  shown 
in  Figure  5.1 . 


Naturally,  the  net  so  constructed  is  exactly  the  Petri  Net  model  of 
the  DM0,  as  Introduced  in  [2],  which  represents  the  information  structure 
within  the  organization.  Let  us  recall  that,  as  long  as  the  model 
corresponds  to  an  admissible  organizational  structure,  we  are  sure  that 


this  net  will  contain  no  loops  (or  directed  oirouits) ,  as  discussed  in 
Chapter  3,  Section  3.5.  At  this  point,  a  Petri  Net  which  satisfies  the 
two  properties: 

(1)  it  Is  acyclic  (i.e..  contains  no  loops) 

(2)  each  place  has  at  most  one  input  transition  and  at  most  one 
output  transition 

is  called  an  Occurence  Net  [20] .  The  Petri  Net  obtained  is  therefore  a 
model  of  an  occurence  net.  The  second  property  is  trivially  satisfied, 
since  our  model  of  the  DM0  is  an  Event-Graph  (each  place  has  exactly  one 
input  and  one  output  transition).  If  we  identify  the  reachable  markings  as 
the  states  of  the  process  and  the  transitions  as  the  steps,  then  property  1 
implies  that  the  process  is  completed  in  a  finite  sequence  of  steps  and 
property  2  implies  that  the  passage  from  one  state  to  another  one  is 

uniquely  defined  (and  can  be  predicted) . 

5.1.2  Characterization  of  Lines  and  Slices 

We  call  marked  places  the  places  of  the  occurence  net  that  contain  one 
token.  We  assume  that  initially,  there  are  no  marked  places  in  the  net. 

Let  us  denote  by  Tx  =  {tx},  where  tx  is  the  input  transition  of  the 

occurence  net.  Assume  now  that  tx  fires  and  let  ?x  be  the  resulting  set  of 
marked  places,  which,  in  that  case,  are  also  the  output  places  of  t}.  In 
the  example  of  Figure  5.1,  Px  =  tpx,  p,J.  Let  T,  be  the  set  of  transitions 
that  are  enabled  by  Px  :  T,  *  {t,»  t,).  Assume  now  that  all  the 

transitions  of  T,  fire  and  let  P,  be  the  new  set  of  marked  plaoes: 
Pj  *  (pa,  P«»  p,).  Similarly  T,  will  be  the  set  of  transitions  that  are 
enabled  by  Pa  and  so  on.  By  iterating  this  procedure,  we  will  eventually 
construct  a  sequence  Tx,  TS,...,TS  such  that:  Ts  -  {tBJ,  where  tB  denotes 
the  output  transition. 

In  the  example  of  Figure  5.1: 
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Let  us  analyze  the  properties  of  the  sequence  (T^)  1-1,.., s  so  constructed. 
We  first  give  the  definition  of  a  slice  and  a  line  of  an  occurence  net,  as 
described  in  [21].  Let  <  denote  the  partial  order  relation,  defined  as 
follows.  Let  a  and  b  be  any  two  nodes  (place  or  transition)  of  the  net, 
then: 


a  <  b  iff  there  exists  a  directed  path  that  goes  from  a  to  b. 


Naturally  the  definition  is  consistent,  since  an  occurence  net  is 
acyclical  (and  therefore  there  could  not  exist  a  directed  path  from  a  to  b 
and  from  b  and  a).  The  definitions  of  a  line  and  a  slice  are: 


Line:  A  line  is  a  set  B  of  nodes  (places  or  transitions)  such  that: 


(1)  for  any  two  nodes  a  e  B,  b  s  B: 


a  <  b  or  b  <  a 


(2)  B  is  maximal:  there  does  not  exist  a  node  c  in  the  net  which 


satifies  (1)  but  does  not  belong  in  B. 


Slioe:  A  slice  is  a  set  D  of  nodes  (places  or  transitions)  such  that: 

(3)  for  any  two  nodes  a  e  D,  b  e  D: 


a  <  b  and  b  <  a 
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(4) 


D  Is  aarimal :  there  does  not  exist  a  node  c  in  the  net  which 


satifies  (3)  but  does  not  belong  in  D. 

Clearly,  the  elements  of  a  slice  are  unordered  (independent  with  respect  to 
<  ) :  the  elements  of  a  slice  characterize  the  concurrent  activities  of  the 
process.  In  contrast,  the  elements  of  a  line  are  strictly  ordered 
(causally  dependent):  these  elements  characterize  the  activities  that  take 
place  sequentially  in  the  process. 

In  our  analysis,  we  are  only  Interested  in  the  slices  that  are  set  of 
transitions,  since  we  want  to  determine  the  tasks  that  are  concurrent  in 
the  process.  By  construction,  each  T^,  as  previously  defined,  is  a  slice. 
Assume,  for  instance,  that  it  would  take  exactly  one  unit  of  time  for  each 
transition  to  fire.  In  that  case,  the  transitions  of  are  all  the 
transitions  that  fire  concurrently  at  the  i-th  instant.  In  fact,  the 
sequence  Tx,  Ta,...,Tg  represents  the  steps  of  the  process  with  the 
transitions  that  fire  concurrently  at  each  step.  The  order  of  these  slices 
determines  precisely  the  partial  order  between  the  various  processing 
tasks,  in  the  following  sense:  If  we  consider  a  line  (a  sequential 
subprocess),  then  each  transition  of  the  line  belongs  to  a  different  slice 
(a  line  cuts  necessarily  any  slice  at  most  once,  as  it  can  be  deduced 
straightforwardly  from  the  definitions)  and  the  order  of  the  corresponding 
slices  corresponds  to  the  order  in  which  the  transitions  fire  in  the 
sequential  subprocess.  This  property  of  Tx,...,Ta  is  fundamental  and  will 
be  used  in  the  next  section  for  determining  the  precise  firing  schedule. 

Remark :  In  our  model,  where  there  is  an  input  transition  (tx)  and  an 
output  transition  (tm),  the  lines  are  exactly  the  directed  paths  that  go 
from  tj  to  tm.  The  determination  of  the  lines  is  straightforward  from  the 
development  in  Chapter  4.  By  adding  the  resource  place  (R#)  of  the  overall 
organization  to  a  line,  we  obtain  a  directed  circuit.  Therefore,  once  we 
have  determined  all  the  directed  circuits  of  a  net  (so  as  to  compute  the 
maximum  throughput  rate),  we  can  obtain  all  the  lines  by  selecting  the 
circuits  that  contain  R#  as  unique  resource  place  and  deleting  it. 


5.2  EXECUTION  SCHEDULE  OF  THE  DMO 


In  this  section,  we  are  going  to  carry  out  the  computation  of  the 
firing  schedule  and  analyze  its  properties.  As  stated  previously,  it  is 
only  possible  to  predict  the  precise  firing  schedule  for  the  deterministic 
process.  Therefore,  in  this  section,  all  the  transition  firing  times  are 
therefore  assumed  deterministic.  We  will  see  later  how  to  use  the  results 
when  processing  times  are  discrete  random  variables  (corresponding  to  the 
switches)  so  as  to  obtain  the  earliest  execution  schedule  and  the  latest 
one.  He  also  assume  here  that  external  inputs  are  always  available  at  the 
entry  of  the  system:  In  other  words,  there  is  a  queue  of  inputs,  such  that 
the  process  occurs  repetitively.  Given  this  assumption,  we  will  be  able 
to  characterize  the  execution  schedule  as  a  function  of  both  the  time  and 
resource  constraints  specific  to  the  organization,  independently  of  the 
arrival  rate  of  inputs.  As  for  the  maximum  throughput  rate,  the  execution 
schedule  obtained  will  characterize  the  dynamic  performance  of  the  DMO. 
Let  us  point  out  that,  here  again,  the  analysis  is  based  on  the  fact  that 
the  Petri  Net  Model  corresponds  to  a  strongly  connected  Event-Graph. 

5.2.1  Assumptions  of  the  Model 

Let  us  recall  that  s£  denotes  the  instant  of  the  n-th  initiation 
firing  of  transition  tA  and  the  firing  time  of  t^.  From  now  on,  we  will 
use  the  sequence  (S^)  for  n=l,2,3, . . .  ,N  and  for  all  the  transitions 
(i=l,...,m)  to  characterize  the  execution  schedule  of  the  system.  For 
clarity,  we  will  denote  the  places  differently:  since  each  place  has 
exactly  one  input  and  one  output  transition,  it  is  possible  to  specify 
precisely  a  place  by  its  unique  input  and  unique  output  transition. 
Accordingly,  we  wi.’  1  denote  by  p^  the  place  whose  input  transition  is  t^ 
and  output  transition  tj.  Likewise,  (res.  M^j  (x))  will  denote  the 
marking  of  place  p^  initially  (resp.  at  instant  t). 

At  this  point,  it  is  emphasized  that  our  objective  here  is  to  analyze 
the  dynamic  behavior  of  the  DMO.  Otherwise  stated,  starting  from  an 


initial  state,  we  want  to  compute  the  sequence  of  initiation  firings  (S£), 
assuming  that  the  process  occurs  repetitively.  Quite  naturally,  we  will 
assume  that  the  processing  starts  at  x  =  0,  so  that  initially  =  0  for 
all  the  places  pAj,  except  fo  the  resource  places,  meaning  that  there  are, 
at  the  Initial  instant,  no  inputs  being  processed.  Naturally,  the  initial 
marking  of  the  resource  places  will  denote  the  resources  available  for 
processing  (and  this  Initial  marking  is  consequently  strictly  positive). 

5.2.2  Feasible  Firing  Schedules 


The  underlying  assumption  concerning  the  firing  schedule  that  we  are 
interested  in  computing,  is  that  the  firing  initiations  occur  as  soon  as 
the  corresponding  transitions  are  enabled.  Otherwise  stated,  the  tasks  are 
executed  as  soon  as  the  corresponding  information  is  available.  This  will 
give  us  indeed  the  best  time  performance  of  the  DM0,  since  the  instant  of 
times  at  which  the  various  tasks  are  executed  are  clearly  the  earliest 
ones.  However,  we  could  imagine  that  some  tasks  are  arbitrarily  delayed 
and,  in  particular,  that  the  firing  schedule  is  set  up  externally  for  the 
system.  At  this  point,  in  treating  (S^)  as  parameters  for  the  system,  the 
problem  necessarily  arises  as  to  whether  or  not  this  firing  schedule  is 
realizable  (or  feasible),  given  the  transition  firing  times  (time 

constraints)  and  the  initial  marking  (resource  constraints).  A  firing 
schedule  would  be  feasible,  if  every  transition  was  actually  enabled  at  the 
instant  the  firing  should  be  initiated.  If  it  were  not  the  case,  the 
firing  schedule  would  be  infeasible,  i.e.,  the  firing  initiation  of  at 
least  one  transition  is  required  to  occur  before  it  is  enabled. 

Ramchandani  [6]  has  given  a  precise  characterization  of  a  feasible 
firing  schedule,  in  case  of  an  Event-Graph: 

Theorem  5.1:  A  firing  schedule  (S^),  i=l,...,m  (where  m  is  the  number  of 

transitions)  is  feasible  if  and  only  if,  for  any  place  p^j  (whose  input 

transition  is,  by  definition,  t^  and  output  transition  t.)  and  for 

. . 
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srMij  2  ^ 

The  above  result  is  intuitively  clear.  If  we  consider  the  instant  of 
the  n-th  termination  firing  of  transition  ti#  i.e.,  the  instant 


then,  transition  tA  has  exactly  "produced"  n  tokens  in  place  p±j  (since 
has  fired  n  times ) .  Now 


n+M 


'.-sj 


ij 


is  the  instant  at  which  the  firing  of  transition  tj  is  initiated  for  the 
(n+M^j)-th  time.  At  this  instant,  transition  t^  has  "consumed"  (n+M^j) 
tokens  in  pAj.  Because  there  were  initially  only  M?j  tokens  in  pAj  (i.e., 
at  x=0),  x,  should  necessarily  occur  after  xx,  i.e.,  after  the  instant  at 
which  transition  tA  has  fired  n  times.  The  inequality 


t.  >  x 


corresponds  exactly  to  relation  (5.1).  It  is  noted  that  a  rigorous 
demonstration  of  Theorem  5.1  is  provided  in  Appendix  B. 

5.2.3  Determination  of  the  Execution  Schedule 

The  result  stated  in  Theorem  5.1  provides  an  easy  way  to  compute  the 
earliest  firing  schedule,  i.e.,  the  one  in  which  the  firings  are  initiated 
as  soon  as  the  transitions  are  enabled.  Let  tj  be  any  transition  and  let 
us  denote  by 
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the  input  places  of  tj.  Following  the  notation  defined  earlier, 

t,  ,  t.,  .  ..t,  correspond  to  the  input  transitions  of  each  of  these  places, 
1  2  r 

as  shown  in  Figure  5.2.  Applying  now  the  previous  result  to  each  of  these 
places  yields : 

0 

n+H:  , 

S.  kJ  >  s“  +  p  for  k  =  1,2, ... ,r  (5.2) 

J  Xk  *k 

At  this  point,  we  recall  that  our  model  of  the  DMO  as  a  Timed  Event-Graph 
includes  an  additional  constraint,  when  compared  to  the  model  of  Timed 
Petri  Nets  introduced  by  Ramchandani  [6].  As  discussed  in  Chapter  2, 
Section  2.4.2,  a  transition  is  not  enabled,  if  it  is  already  executing, 
which  can  be  modeled  by  adding  a  self-loop  to  each  transition.  The 
corresponding  place  should  be  denoted  by  Pjj  (which  means  that  its  input 
and  output  transition  is  tj)  with  an  initial  marking  Mjj  =  1.  Applying  the 
inequality  (5.1)  to  this  place  yields: 

■r1  -  s3  *  ^  <5-3’ 

This  is  therefore  the  additional  constraint  that  should  be  included  in  our 
modeling,  for  specifying  a  feasible  firing  schedule. 

Let  us  now  assume  that  n,  the  number  of  repetitions,  is  large  enough, 
so  that: 

n  >  M*  for  k  =  l,...,r.  (5.4) 

AkJ 

The  previous  Inequalities  (5.2)  and  (5.3)  can  be  written: 


Figure  5.2  Input  Places  and  Transitions  of  t. 


Sn  >  SQ  Mikj  + 


for  k  =  1,2,. . . ,  r . 


(5.5) 


*  *r  *  ^ 


(5.6) 


The  earliest  n-th  firing  initiation  of  transition  tj  is  then  determined  by: 


n-M®  , 


n  w 

S  .  =  max 


k=l ..... r 


u  nj  ^  - 

(S  ,  ,Sj  *  Kj) 


(5.7) 


By  construction,  such  a  firing  schedule  is  feasible  and  it  is  necessarily 
the  earliest  one.  The  instant,  as  defined  in  (5.7),  corresponds  to  the 
instant  at  which  the  transition  is  enabled  for  the  n-th  time.  As  expected. 


the  places  p^j  for  which  the  maximum  is  obtained  have  received  a  token  at 


the  latest  instant  of  time,  therefore  enabling  transition  tj. 


Now  if. 


v  v-* 


n  i  M, 


for  some  k  (k  e  {1,2, . . .,r}) 


the  equality  (5.7)  still  holds,  except  that,  in  the  determination  of  the 
maximum,  we  do  not  take  into  account  the  corresponding  values  of  k. 
Indeed,  it  means  that  the  corresponding  places  pA  j  have  still  a  marking 
strictly  positive  (from  the  initial  marking)  and  therefore  these  places  do 
not  disable  transition  t4. 


Relation  (5.7)  turns  out  to  be  an  implicit  equation,  which  permits  the 
computation  of  the  firing  schedule  by  iteration  on  n,  the  number  of  firing 
repetitions.  If  the  initial  marking  of  all  the  places  were  strictly 
positive  (i.e.,  >  0  for  all  the  places  P^j)*  S°  (for  jsl,2,...,m)  could 
be  easily  computed,  using  (5.7),  from  the  determination  of  (S^),  for 
k=l ,2, . . . ,n-l  and  1=1,2, ... ,m.  This  is  however  not  the  case  in  the  model, 
given  the  assumptions  made  in  Section  5.2.1;  As  already  stated,  only  the 
resource  places  have  a  non-null  initial  marking.  For  a  pair  of  transitions 
(tiftj)  connected  by  the  place  p±j  such  that  M^j  =  0,  we  obviously  need  to 
determine  first  S“  to  compute  S°,  using  (5.7).  This  condition  implies  that 
the  computation  of  S1?  for  the  various  transitions  (j  =  l,2,...,m)  should  be 

J 

carried  out  following  a  precise  order.  As  it  will  be  disoussed  now,  this 
order  is  characterized  by  the  slices  of  the  net,  as  introduced  in  Section 
5.1.2. 

Given  our  assumptions  that  the  processing  starts  at  t=0,  the  places 
Pjj  for  which  =  0,  as  stated  in  Section  5.2.1,  are  all  the  places  of 
the  net,  except  the  resource  places.  In  other  words,  these  places  are  the 
places  of  the  occurence  net  that  describes  the  precise  sequence  of 
activities  in  the  process,  as  constructed  in  Section  5.1.1.  Let  us 
consider  then  the  sequence  of  slices  Tx,...?s,  which  has  allowed  us  to 
characterize  the  sequential  and  concurrent  activities  of  the  process. 

Any  pair  of  transitions  ( t^  ,  tj),  (connected  by  place  pA j )  belongs  to 


a  line  (since  trivially  <  t  J  and  both  transitions  are  therefore 

elements  of  two  different  slices,  let  us  say  T&  and  T^.  As  we  have 
explained,  the  order  of  the  slices  reflects  the  partial  order  between  the 
transitions,  so  that  we  have  necessarily:  a  <  b.  It  turns  out  that,  if  we 
have  determined  the  instants  of  the  n-th  firing  initiation  for  all  the 
transitions  belonging  to  Tx, . . .  .Tj^,  then  we  will  be  able  to  derive,  using 
relation  (5.7),  the  corresponding  Instant  for  the  transitions  belonging  to 
Tk.  Now,  S“,  corresponding  to  Ta  =  {ta}  (input  transition)  is  never  a 
problem,  since  tx  has  a  unique  input  place,  which  is  the  resouroe  place, 
R#,  of  the  overall  organization.  Hence,  M^,  where  pffll  denotes  the 
resource  place  R0,  given  the  notation  defined  previously  (tm  being  the 

output  transition),  is  always  strictly  positive  and  a  direct  iteration  on  n 
generates  Sj,  i.e.,  using  (5.7): 

0 

=  max  (S  +  n  ,  1  +  n  ) 

i  m  mix 

The  order  in  which  the  computation  of  Sj,  for  j  =  l,...,m,  should  be 

carried  out,  is  determined  by  the  partial  order  defined  by  the  slices 
T1,Ta, . . . ,T3.  Let  us  point  out  that,  for  the  transitions  belonging  to  the 
same  slice,  the  order  has  strictly  no  importance  since  the  transitions  fire 
concurrently. 

In  brief,  relation  (5.7)  allows  the  determination  of  the  execution 
schedule  by  iteration  on  n,  which  corresponds  to  the  number  of  firing 

repetitions  from  the  initial  instant.  At  each  iteration,  however,  the 
instants  at  which  the  different  tasks  are  executed  should  be  determined, 
following  the  order  with  which  the  operations  occur  in  the  process,  i.e., 
the  partial  order  defined  by  the  slices.  A  detailed  description  of  the 
corresponding  algorithm  is  provided  in  Appendix  C. 

Finally,  we  have  been  able  to  determine  the  best  execution  schedule 
(i.e.,  the  earliest  one)  of  the  DM0,  as  a  function  of  both  the  time 
constraints  and  resource  constraints.  Indeed,  in  relation  (5.7),  are 


the  task  processing  times  and,  as  we  have  said,  the  initial  marking  of 
the  net  models  the  resource  constraints.  The  execution  schedule  so 
obtained  characterizes  the  dynamic  behavior  of  the  DMO:  starting  from  the 
Initial  state  at  x  -  0  (where  no  Inputs  are  processed  and  all  the  resources 
are  available),  the  process  is  assumed  to  occur  repetitively  (external 
inputs  are  assumed  queueing).  Naturally,  the  results  developed  here  are 
only  valid  for  the  deterministic  process,  as  we  have  said  in  the 
introduction.  Now,  in  the  case  where  processing  times  are  discrete  random 
variable,  it  is  not  possible  to  predict  precisely  the  instants  of  time  at 
which  the  various  operations  take  place.  However,  by  considering  on  the 
one  hand  the  shortest  processing  time  for  each  random  variable  and  on  the 
other  hand  the  longest  processing  time,  we  obtain  two  deterministic  models. 
Computing  the  firing  schedule  yields  in  the  first  case  the  earliest 
execution  schedule  possible  and  in  the  second  case  the  latest  one. 
Naturally,  the  real-time  firing  schedule  is  in  between  these  two  bounds. 
Accordingly,  we  can  say  that  the  earliest  execution  schedule  characterizes 
the  best  possible  performance  of  the  DMO  and  the  latest  one,  the  worst 
possible  performance.  In  addition,  we  can  easily  compute  the  probability 
of  each  occuring.  From  the  probability  distribution  assigned  to  the 
possible  processing  times,  we  know  the  probability  of  the  shortest  (resp. 
longest)  processing  time  for  each  transition.  By  computing  the  product  of 
these  probabilities  for  all  the  transitions,  we  obtain  precisely  the 
probability  that  the  best  (resp.  worst)  performance  occurs. 

5.2.4  Properties  of  the  Execution  Schedule 

In  this  section,  we  carry  out  the  analysis  of  the  firing  schedule  and 
investigate  the  properties  of  the  steady  state  process.  The  fundamental 
result  presented  here,  which  characterizes  the  periodicity  of  the  steady- 
state  process,  has  been  obtained  by  Chretienne  [16].  Let  us  first  give  the 
definition  of  the  K-periodicity  that  he  introduced  in  [16]. 


Definition:  A  sequence  SQ  (n  =  1,2,3,...)  is  said  to  be  K-per iodic,  with 
period  n,  iff  for  all  n  2  1: 


jn  =  K  q  +  r 
(o  <  r  <  K 


S 


n 


S  +  q  n 
r  * 


(5.8) 


If  we  Identify  the  sequence  SQ  as  the  Instant  of  times  at  which  the  same 
operation  occurs  successively,  it  means  that  the  K  successive  occurences 
repeat  periodically,  as  it  can  be  visualized  on  Figure  5.3. 
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Figure  5.3  K-Periodic  Sequence  with  period  n 

When  K  =  1,  it  corresponds  to  the  usual  periodicity  (also  called  strong 
periodicity  by  Ramchandani  [6]).  The  ratio  K/n  represents  the  (average) 
rate  (or  frequency)  with  which  the  operation  occurs.  We  will  denote  this 
rate  by: 


so  that  (5.8)  can  be  rewritten: 


S  =  S  +  q  5  (5.9) 

n  r  H  4 

The  fundamental  property  that  characterizes  the  dynamic  behavior  of  the 
system  is  the  following  one: 

Theorem  5.2  [16]:  The  firing  schedule  (as  determined  previously)  becomes 
K-perlodlc  after  a  finite  number  of  firing  repetitions,  which  means  that 


114 


m  (where  a 


there  exist  E  and  N0  such  that,  for  n  1  N#  and  for  j  -  1,2 
is  the  number  of  transitions): 

n-N#  =  Eq  +  r 

Sn  -  S1J'»+r  +  q  I-  (5.10) 

v  0  i  r  <  E  J  J  * 

The  rate  i  is  equal  to  the  minimum  average  circuit  processing  rate  as 
defined  in  Chapter  4,  Section  4. 2 .2.1. 

Let  us  point  out  first  that  this  Theorem  holds  only  for  the  type  of 
Petri  Nets  used  for  modeling  the  DM0,  i.e.,  strongly  connected  Event- 
Graphs,  and  for  deterministic  processing  times.  Let  us  now  comment  on  the 
result. 

This  result  completes  the  development  carried  out  in  the  previous 
chapter,  concerning  the  maximum  throughput  rate.  At  it  was  discussed 
there,  the  minimum  average  circuit  processing  rate  bounds  the  allowable 
rate  at  which  inputs  can  arrive,  in  the  sense  that  above  this  rate,  the 
queue  of  inputs  at  the  entry  of  the  system  will  necessarily  grow  to 
infinity  over  time.  This  theorem  makes  sure  that,  in  the  deterministic 
case,  this  rate  corresponds  effectively  to  the  maximum  throughput  rate  that 
can  be  achieved  by  the  DM0,  as  stated  in  Chapter  4,  Section  4. 2 .4.1. 

In  addition,  we  have  a  further  result:  We  know  that  this  throughput 
rate  is  actually  achieved  within  a  finite  amount  of  time  or,  in  other 
words,  after  a  finite  number  of  repetitions  of  the  process.  This  result  is 
not  surprising;  A  Petri  Net  which  is  live  and  bounded  (which  is  the  case  of 
our  model)  has  necessarily  a  finite  number  of  reachable  markings,  i.e.,  the 
forward  marking  class,  as  defined  in  Chapter  2,  Section  2.1.5  is  finite. 
Therefore,  there  is  at  least  one  state,  sQ,  (i.e.,  one  marking)  which 
occurs  cyclically  in  the  repetitive  process.  Furthermore,  in  case  of  an 
Event-Graph,  the  set  of  all  reachable  markings  from  any  Initial  marking  is 
uniquely  defined  and  can  be  predicted,  as  stated  in  Section  5.1.1. 
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Consequently,  all  the  states  that  take  plaoe  between  the  successive 
occurences  of  sQ  also  repeat  cyclically.  It  means  that,  once  the  state  s0 
is  reached  for  the  first  time,  the  reachable  states  repeat  periodically, 
determining  the  steady-state  process.  However,  as  Theorem  5.2  proves  it, 
the  period  in  real  time  (i.e.,  given  the  transition  firing  times)  turns  out 
to  be  of  a  more  complex  type  than  the  simple  (or  strong)  periodicity. 

From  the  determination  of  the  firing  schedule  (S°)  and  of  the  rate  i, 
computed  as  the  minimum  average  circuit  processing  rate,  it  is  now 
straightforward  to  deduce  E,  the  periodicity  factor.  However,  an 
interesting  result  further  obtained  in  [16],  states  that  if  all  the 
critical  ciroults,  as  defined  in  Chapter  4,  Section  4. 2. 2. 2,  contain  only 
one  token,  then  E  *  1.  In  that  case,  it  means  that  the  steady-state  process 
is  strongly  periodic,  i.e.,  each  transition  fires  at  the  same  regular 
Interval  of  time.  Otherwise  stated,  the  processing  rate  is  constant  (for 
the  steady-state).  Interestingly  enough,  for  the  model  of  the  DMO  analyzed 
here,  it  would  be  this  case  if: 

-  either  the  critical  circuits  contain  only  one  resource  place  and  this 
resource  place  has  initially  one  token  (meaning  that  the  corresponding 
DM  or  the  overall  organization  can  only  process  one  input  at  time) 

-  or  the  critical  circuits  correspond  to  individual  transitions  (recall 
that  a  transition  contain  in  fact  a  self  loop  with  one  token). 

5.3  MEASURES  OF  PERFORMANCE  FROM  THE  EXECUTION  SCHEDULE 

In  the  first  two  sections,  we  have  been  able  to  characterize  the 
dynamic  behavior  of  the  DMO,  having  first  determined  the  sequential  and 
concurrent  operations  of  the  process  and  then  computed  the  precise 
execution  schedule  for  any  number  of  repetitions  of  the  prooess.  In 
addition,  we  know  that  the  steady  state  process  is  E-periodic,  with  a 
processing  rate  f  which  corresponds  to  the  minimum  average  circuit 
processing  rate.  Recall  that  u.is  analysis  is  appplicable  under  the 


following  assumptions: 


(i)  the  process  is  deterministic,  which  implies  that  all  the  task 
processing  times  are  deterministic. 

(ii)  the  process  occurs  repetitively  from  the  initial  instant  t  ■  0. 

(ill)  The  tasks  are  executed  as  soon  as  the  corresponding  information 
is  available,  which  means  that,  in  the  Petri  Net  model,  the 
firing  initiations  occur  as  soon  as  the  transitions  are  enabled. 
The  execution  schedule  so  obtained  is  the  best  possible,  in  the 
sense  that  the  instants  of  time  at  which  the  various  tasks  are 
processed  are  the  earliest  ones. 

In  the  case  where  the  processing  times  are  discrete  random  variable, 
as  it  is  the  case  for  switches,  we  can  extend  this  analysis  by  considering 
the  deterministic  process  associated  with  the  shortest  processing  times  and 
the  one  associated  with  the  longest  processing  times  of  each  transition. 
In  the  first  case,  we  obtain  the  earliest  execution  schedule  possible  and 
in  the  second  case  the  latest  execution  schedule  possible.  Naturally,  the 
real-time  process  is  in  between  these  two  bounds. 

It  is  now  Interesting  to  analyze  qualitatively  the  results  obtained 
and  derive  related  measures  of  performance.  There  are  actually  two 
measures  that  are  particularly  interesting  to  analyze:  The  execution 
schedule  (Sx)  of  the  input  transition  and  the  execution  schedule  (S°)  of 
the  output  transition. 

5.3.1  Execution  Schedule  of  the  Input  Transition 

Each  time  that  the  firing  of  the  input  transition  tx  is  Initiated,  it 
means  that  the  DM0  starts  processing  a  new  input.  Given  the  assumptions 
made  for  determining  the  firing  schedule,  Sx  represents  the  earliest 
instant  at  which  the  n-th  input  can  be  processed;  i.e.,  whatever  may  be. 
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the  instant  of  times  at  which  inputs  arrive  successively,  the  schedule  (s“) 
(n  *  1,2,3,...)  determines  the  earliest  instants  at  which  the  DMO  is  able 
to  handle  these  inputs.  Because  this  schedule  is  not  necessarily  strongly 
periodic  but  only  K-periodic  in  its  steady  state,  we  do  not  generally 
obtain  the  best  performance  by  regulating  inputs  at  the  entry  of  the 
system,  such  that  the  arrival  rate  be  constant.  Suppose  for  example  that 
K=2  and  that 


SN*+1  =  SN«  +  —  SNo+* 

5i  5»  2|  '  5i 


s".  ♦  £.... 


(5.11) 


as  shown  on  Figure  5.4  (Recall  that,  after  N#,  repetitions  the  process  is 
in  its  steady-state). 


Occurences 


N0+a  V*  V4 


Instants 


sN«+1 

i 


SN.«  SN.~  s».*< 


Figure  5.4  2-Periodic  Firing  Schedule 


If  the  interarrival  time  of  inputs  was  constant  and  equal  to  1/#  (f 
being  the  maximum  throughput  rate),  the  processing  of  every  two  inputs 
would  be  delayed  by  l/2i,  compared  to  what  it  could  be  if,  for  example,  two 
inputs  were  arriving  simultaneously  at  regular  interval  of  time  2/f.  In 
fact,  so  as  to  reduce  at  the  minimum  possible  the  processing  delays,  K 
inputs  should  arrive  simultaneously  at  regular  intervals  of  time  duration 
K/i.  Naturally,  this  supposes  that  the  buffer  capacity  at  the  entry  of  the 
system,  is  at  least  equal  to  K. 


5.3.2  Execution  Schedule  of  the  Output  Transition 


Let  us  denote  by: 


r11  -  Sn  +  |i 
mom 


where  pm  is  the  processing  time  of  the  output  transition  tD.  is  the 


instant  of  the  n-th  firing  termination  of  tm,  the  output  transition.  Each 


time  that  the  firing  of  t_  is  terminated,  it  means  that  the  DMO  has 


completed  the  processing  of  an  input  and  has  produced  the  corresponding 
response.  Now,  let  us  assume  that  n  inputs  arrive  simultaneously  at  t  * 
0.  Then  T^  (k  =  l,2,...,n)  represents  the  instant  of  time  at  whloh  the 
processing  of  the  k-th  input  (defined  by  the  order  in  which  inputs  are 
handled)  is  completed:  For  that  reason,  we  shall  say  that  the  schedule 


(T^)  represents  the  dynamic  response  time  of  the  DMO.  At  this  point,  we 


should  emphasize  the  difference  existing  between  the  time-delay  of  the  DMO, 
as  computed  in  [3]  and  this  measure.  In  [3],  the  time-delay  measures  the 
duration  from  input  entering  the  system  to  the  response  taking  place.  This 
time-delay  is  a  static  measure,  in  the  sense  that  the  system  is  supposed  to 
process  one  input  at  a  time.  Quite  obviously,  this  measure  would 
correspond  to  T*  =  S*  +  in  our  model,  i.e.,  the  amount  of  time  it  takes 
to  process  the  first  input,  since  the  processing  is  supposed  to  start  at 


x=0.  Likewise,  it  would  also  correspond  to  the  difference: 


d  =  T?  -  Sn  =  Sn  -  Sn 
n  m  *  m  * 


+  p 

x  'm 


(5.12) 


for  n  2  N#  (where  N,  determines  the  number  of  repetitions  above  which  the 
process  is  in  its  steady  state)  and  assuming  that  the  resource  plaoe  of  the 


overall  organization,  R0,  contains  only  one  token.  As  we  have  discussed  in 
Chapter  3,  Section  3.3.2,  it  means  indeed  that  the  DMO  can  only  handle,  in 
that  case,  one  input  at  a  time.  This  intuitive  result  can  actually  be 
proved  rigorously.  Let  us  consider  all  the  lines  (or  directed  path  of  the 


occurence  net,  as  defined  in  Section  5.1,  for  which  the  sum  of  the 
transition  firing  times  is  maximal.  By  adding  the  resource  place  R0,  we 
would  trivially  obtain  all  the  crl tical  circuits,  precisely  because  R0 
contains  only  one  token.  In  that  case,  the  maximum  average  circuit 
processing  time,  i.e.,  l/$,  is  equal  to  the  maximum  time-delay  over  all  the 
possible  directed  paths,  denoted  by  Dj^.  Now,  since  the  token  content  of 
the  critical  circuits  is  one,  the  periodicity  factor  K  equals  one  (Section 
5.2.4),  which  means  that  the  steady  state  is  strongly  periodic.  In 
particular,  for  n  2  N0: 


sn+1  =  sn  +  i 

i  i  i 


+  D, 


Max 


(5.13) 


Applying  relation  (5.7)  to  the  place  pal,  which,  following  the  notation 
used  previously,  denotes  the  resource  place  R#  of  the  overall  organization, 
yields: 


sTl  -  <sa +  +  **x) 


since  Mffll  =  1,  by  assumption. 

Now,  the  expression,  for  any  n  2  1: 


Sn  +  ji  -  (Sn  +  ji  ) 

m  m  i  i 


(5.14) 


represents  the  amount  of  time  it  takes  to  complete  the  processing  of  the 
n-th  input,  once  the  input  operation  (modeled  by  t4)  is  completed.  Hence, 
this  quantity  is  always  positive,  so  that  in  fact,  (5.14)  can  be  written: 


Sn  +  m  =  T11 
m  mm 


Using  relation  (5.13),  we  deduce: 
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s“  +  Du  *  T° 
*■  Max  in 


and  therefore,  given  the  definition  (5.12): 


T11  -  SD 
m  i 


Max 


Q.E.D. 


In  that  case,  the  amount  of  time  it  takes  for  the  DMO  to  complete  the 
processing  of  n  inputs  is  clearly  equal  to  nD^.  However,  it  should  be 
clear  that,  in  the  general  case,  where  the  DMO  can  handle  more  than  one 
input  at  a  time,  the  amount  of  time  it  takes  to  oomplete  the  processing  of 
n  inputs  may  actually  be  much  less.  With  the  assumption  that  the 


processing  starts  at  x=0,  this  amount  of  time  is  precisely  given  by  TjJ,  the 


dynamic  response  time. 


Now,  we  can  define: 


RT  =  lim 


m 


n->+« 


lim 
n— )  +• 


Sn  +  m 
m  m 


(5.15) 


as  the  average  response  time  of  the  DMO.  From  Chapter  4,  Section  4.2.1,  we 
know  that  this  limit  will  precisely  be  equal  to  a,  which  is  also  equal  to 
1/1.  Hence  the  average  response  time  of  the  DMO  corresponds  to  the  maximum 
average  oircuit  processing  time. 


At  this  point  we  should  make  a  clear  distinction  between  what  we  have 
defined  as  the  average  response  time  and  what  would  correspond  to  the 
average  processing  time  (the  equivalent  of  the  time-delay  defined  in  13]), 
that  is: 
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d  -  11a  d  -  11m  CT  -  S  ) 

.  __  n  .  __  m  » 

n->  +®  n->  +* 


(5.16) 


which  can  also  be  written  as: 


d  +  d  +. . .+  d  _ 
n  n+*  n+K-i 


for  any  n  1  N# 


because  of  the  K-periodicity  of  the  steady-state  process,  d  is  an  average 
measure  of  the  amount  of  time  It  takes  to  oomplete  the  processing  of  any 
individual  input.  In  the  general  case,  the  value  of  d  is  necessarily 

greater  or  equal  to  the  maximum  time-delay  of  the  directed  paths, 

previously  denoted  by  Dj^,  i.e.,  the  corresponding  value  when  the  DMO 

handles  only  one  input  at  a  time.  It  may  be  indeed  greater  precisely 
because  more  than  one  input  can  be  processed  at  the  same  time,  so  that 
additional  delays  occur,  due  to  the  unavailability  of  the  limited 

resources.  In  fact,  d  can  be  characterized  quite  easily.  In  the  steady- 
state  process,  inputs  are  processed  simultaneously  <c  is  the  initial 
marking  of  the  resource  place,  R# ,  of  the  overall  organization) .  Since  d 
is  the  average  processing  time  of  each  input,  every  transition  fires  at  the 


average  rate: 


Since  I  is  precisely  the  processing  rate: 


so  that,  finally: 


WWW 


Let  us  now  prove  that  we  have  necessarily:  d  2  ^Max*  43  stated  earlier,  by 
adding  the  resource  place  R#  to  the  directed  paths  (i.e.,  lines)  for  which 
the  time-delay  is  maximal  and  therefore  equal  to  PM^V,  we  obtain  a  circuit. 
The  token  content  of  such  a  circuit  is  M^1,  the  corresponding  average 
circuit  processing  rate  is: 


M 


m^ 


Max 


By  definition,  4  is  the  minimum  average  circuit  processing  rate,  therefore 


4  < 


M 


mi 


Max 


Given  (5.17),  we  deduce  immediately: 


d  2  P 


Max 


Q.E.D. 


In  brief,  we  have  defined  the  dynamic  response  time  as  the  amount  of 
time  it  takes  to  complete  the  processing  of  n  inputs,  which  would  arrive 
simultaneously  at  the  initial  instant  (i.e.,  t=0).  For  large  number  of 
Inputs  (when  n  goes  to  infinity),  dividing  by  n  yields  the  average  response 
time  of  the  organization  and  it  was  proved  that  the  average  response  time 
is  precisely  determined  by  the  maximum  average  circuit  processing  time. 
However,  the  average  processing  time,  d,  for  each  individual  input,  i.e., 
the  average  time-delay  from  one  input  entering  the  system  to  the  response 
taking  place  is  determined  by: 
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In  other  words,  d  Is  equal  to  the  total  number  of  inputs  being  processed 
simultaneously  times  the  average  response  time,  as  would  be  expected. 


5 .4  CONCLUSION 

In  this  chapter,  we  have  carried  out  the  analysis  of  the  execution 
schedule  for  the  DMO  model.  We  have  first  described  the  sequential  and 
concurrent  activities  of  the  process,  using  the  slices  of  a  Petri  Net.  A 
method  has  been  developed  to  compute  the  firing  schedule,  for  the  case 
where  all  the  transition  processing  times  are  deterministic.  The  property 
of  K-periodicity  that  characterizes  the  steady-state  process  has  also  been 
presented.  Finally,  we  have  described  different  MOPs  derived  from  the 
execution  schedule.  This  analysis  will  now  be  illustrated  in  the  next 
chapter  with  a  specific  example  of  DMO. 


CHAPTER  VI 


APPLICATION 

In  this  chapter,  an  example  of  a  five  DM  organization  is  used  to 
illustrate  the  results  presented  in  Chapter  4  and  5.  The  structure  of  the 
organization  is  fixed,  but  the  resource  and  time  constraints  vary.  In  each 
case  considered,  the  execution  schedule  and  the  maximum  throughput  rate  are 
computed.  The  measures  of  performance  obtained  are  then  analyzed  and 
compared.  The  task  processing  times  considered  in  the  examples  are  all 

deterministic. 

6.1  PRESENTATION  OF  THE  MODEL 

The  Petri  Net  representation  of  the  organizational  structure 
considered  is  shown  in  Figure  6.1.  The  five  DM  organization  models  a 
generalized  submarine  ship  control  party  performing  emergency  control 
tasks.  This  particular  organization  has  been  extensively  described  and 
analyzed  in  [22],  using  the  information  theoretic  framework  as  a  way  to 
evaluate  the  performance  of  the  organization.  More  precisely,  the 

performance  evaluation  carried  out  in  [22]  is  based  on  measures  of  the 
workload  of  the  organization  members  and  the  accuracy  of  the  response  [4], 
The  goal  here  is  to  evaluate  the  time-related  performance  of  the 
organization. 

Let  us  now  comment  on  the  -  structure  of  the  organization. 
Interestingly  enough,  the  organization  has  both  hierarchical  and  parallel 
aspects,  where  the  DMs  play  in  fact  different  roles:  DM1  acts  as  a 
supervisor,  DM2  acts  as  a  coordinator  and  DM3,  DM4  and  DM5  act  as 

subordinates.  Indeed,  DM3,  DM4  and  DM5  transmit  their  situation 

assessments  to  DM2  and  receive  commands  inputs  from  this  DM.  However, 
before  issuing  these  command,  DM2  fuses  the  information  transmitted  by 
these  three  DMs  with  his  own  assessment  and  transmits  the  resulting 


Figure  6.1  Petri  Net  Model  of  tbe  Five  DM  Organisation 


information  to  DM1,  the  suprrvisor.  DM1  combines  this  information  with  his 
own  situation  assessment  and  then  produces  a  command  to  DM2.  DM2  prooeeds 
to  select  a  response  that  he  then  transmits  back  to  the  subordinates  (DM3, 
4,  and  5)  as  a  command.  Finally,  DM3,  4  and  5  select  in  turn  a  response 
and  the  three  produce  the  output  which  constitutes  the  organization  s 
response.  At  this  point,  we  should  note  in  particular  two  characteristics 


of  the  organizational  structure.  The  first  one  is  that  DM4  and  DM5 

directly  select  a  response  from  the  command  inputs  received  by  DM2,  without 
processing  the  information  relative  to  their  situation  assessment.  The 

second  one  is  that  the  situation  assessment  transmitted  by  DM3  is  not 

combined  with  the  other  assessments  (from  DM4  and  DM5)  at  the  IF 

(Information  Fusion)  stage  of  D)C.  What  actually  occurs  is  that  the 
situation  assessment  of  DM3  is  preprooessed  by  DM2.  In  fact,  the 
transition  t4  that  models  the  SA  stage  of  DM2  corresponds  to  an  aggregate 
representation  of  the  subnet  shown  on  Figure  6.2. 


Figure  6.2  Petri  Net  Representation  of  the  SA  Stage 
with  Preprocessor  (DM2) 

In  the  represenation  of  Fig.  6.2,  the  first  transition  corresponds  to 
the  preprocessor  task  and  the  second  transition  models  the  usual  Situation 
Assessment  stage.  Now,  the  external  inputs  received  by  DM2  (in  place  pa#) 
are  first  preprocessed  together  with  the  situation  assessment  transmitted 
(in  place  pt4)  by  DM3.  The  resulting  output  from  the  preprocessor  is  then 
processed  at  the  SA  stage  of  DM2.  In  that  case,  the  preprocessor  acts  as  a 
decision  aid  to  DM2.  Note  that,  in  the  aggregated  Petri  Net  model  of  the 
DMO,  as  shown  on  Figure  6.1,  some  transitions  include  in  reality  switches, 
where  a  selection  does  occur  among  different  algorithms  possible.  The 
corresponding  transitions  are  transition  t4  (resp.  t7)  modeling  the  SA 
stage  (resp.  RS  stage)  of  DM2,  and  transition  txl,  modeling  the  RS  stage  of 
DM3. 

Finally,  the  Petri  Net  representation  of  this  DMO,  that  corresponds  to 
the  model  developed  in  Chapter  3  of  the  thesis,  is  shown  on  Figure  6.3. 


Figure  6.3  Model  of  the  Five  DM  Organization  with  Limited  Resources 


The  different  resource  places  have  been  added,  i.e..  the  resource  place  of 
each  DM,  between  the  RS  stage  and  the  SA  stage,  and  the  resource  place  of 
the  overall  organization,  between  the  output  transition  and  the  input 
transition  of  the  net.  Based  on  this  model,  the  performance  evaluation  of 


the  organization  will  now  be  conducted,  following  the  analysis  carried  out 
in  Chapter  4  and  Chapter  5. 


6.2  PERFORMANCE  EVALUATION  OF  THE  FIVE  DM  ORGANIZATION 

In  order  to  evaluate  the  performances  of  the  organization,  given  the 
resource  and  time  constraints,  it  is  first  necessary  to  describe  the 
structural  characteristics  of  the  model,  as  determined  by  the  directed 
circuits  and  the  slices  of  the  Petri  Net.  As  discussed  in  Chapter  4, 
Section  4.2.4,  the  determination  of  all  the  directed  elementary  circuits  is 
required  so  as  to  compute  the  maximum  throughput  rate.  Likewise,  the 
determination  of  the  slices  which,  as  said  in  Chapter  5,  Section  5.1, 
precisely  describe  the  concurrent  and  sequential  activities  of  the  process, 
is  necessary  to  compute  the  execution  schedule  of  the  organization. 
Accordingly,  we  will  first  describe  these  charateristics  for  the  model 
studied  here. 

6.2.1  Structural  Characteristics  of  the  Model 

Both  the  directed  circuits  and  the  slices  are  determined  from  the 
Incidence  Matrix.  Recall  that  the  Incidence  Matrix,  C,  characterizes  the 
structure  of  the  Petri  Net  in  the  following  manner:  THe  columns  of  the 
matrix  correspond  to  the  transitions  of  the  net  and  the  rows  to  the  places. 
The  element  C<  <  of  the  matrix  is  such  that : 

if  pA  is  an  input  place  of  t^ 
if  p.^  is  an  output  place  of  t^ 
otherwise 


The  Incidence  Matrix  corresponding  to  the  Petri  Net  model  of  Figure  6.3  is 
shown  on  Table  6.1.  From  the  Incidence  Matrix,  it  is  possible  to  determine 
all  the  directed  circuits  of  the  net,  using  the  algorithm  developed  by 
Alawain  and  Toudlc  [19],  that  is  presented  in  Appendix  A  of  the  thesis. 
The  list  of  the  directed  circuits  obtained  is  provided  in  Table  6.2.  The 
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0 

0 

0 

0 

0 

0 

0 

p  2  3 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

P24 

0 

0 

0 

0 

-1 

0 

0 

0 

0 

0 

p25 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

p26 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

RO 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

R 1 

0 

-1 

1 

0 

0 

0 

0 

0 

0 

0 

R2 

c 

0 

0 

-1 

0 

0 

1 

0 

0 

0 

R3 

0 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

R4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

R5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

u  u 
0  0 
0  0 
0  0 
0  0 
0  0 
1  0 
0  0 
1  0 
1  -1 


u  u 

0  0 


0  -1 

0  0 


0  -i 


1  -1 


0  -1 

0  0 


TABLE  6.2.  LIST  OF  THE  DIRECTED  CIRCUITS 


1»  p2-t3-Rl-t2-p2 

2«  p4-t5-p5-t6-p6-t7-R2-t4-p4 

3:  pl-t2-p2-t3-p3-t6-p6-t7-R2-t4-p4-t5-pl 

4 i  p4-t5-p5-t6-p6-t7-p7-t 10-p 10— t 1 1-R3— t8— pl8— t4-p4 

5:  pl-t2-p2-t3-p3-t6-p6-t7-p7-t 10-p 10-t 1 l-R3-t 8-p 18-t4-p4-t 5-p  1 

6:  p8-t9-p9-tl0-pl0-tll-R3-t8-p8 

7:  p2-t3-p3-t6-p6-t7-p7-tl0-pl0-tll-pll-t!8-R0-tl-pl9-t2-p2 
8:  p4-t5-p5-t 6-pA-t 7-p7-t  1 0-p  10-tll-pU-tl 8-RO-t  1  -p20-t4-p4 
9:  p 1 -t2-p2-t3-p3-t6-p6-t7-p7-t 10-p 10-t 1 1  — p 1 1 -t 18-RO-t l-p20-t4-p4-t 5-p J 
10 :  p4-t5-p5-t6-p6-t7-p7-tl0-pl0-tll-pll-tl8-R0-tl-p21-t8-pl8-t4-p4 
11:  pl-t2-p2-t3-p3-t6-p6-t7-p7-t 10-p 10-t 1 1-pl 1-tl 8-RO-t 1 -p 2 1 - t 8-p 1 8- t 4 -p 4- t 5-p 1 
12:  p8-t9-p9-tlO-plO-tll-pll-tl8-RO-tl-p21-tB-p8 
13:  p5-t6-p6-t7-p7-t 10-p 1 0-t 11-pll-tl 8-RO-t l-p22-tl2-p24-t5-p5 
14:  p 1  — t2-p2— 1 3-p3— t6-p6— t7-p7— 1 1 0-p 1 0-t 1 1  — p 1 1-t 18-RO-t 1 -p22-t 12-p24-t5-p 1 
15:  p9-t 10-p 10-t 1 1-pl 1-t 18-RO-t 1 -p 22- t 1 2 -p 1 2- t 9 -p 9 
16:  p5-t6-p6-t7-p25-t 1 3-pl3-t 1 4-R4-t 1 2-p24-t5-p5 
17:  pl-t2-p2-t3-p3-t6-p6-t7-p25-tl3-pl3-tl4-R4-tl2-p24-t5-pl 

18:  p4-t5-p5-t6-p6-t7-p25-t 1 3-p 1 3-t 1 4-R4-t l 2-p 1 2-t 9-p9-t 10-p 10-t 1 1 -R 3 - t 8 -p 1 8- t 4 -p 4 

19:  p 1-t 2-p2-t3-p3-t 6-p6-t7-p25-t 13-pl 3-tl4-R4-tl 2-p l2-t9-p9-t 10-p 10-t 11 -R3-t0-p 18- 
t4-p4-t 5-p 1 

20:  p2-t3-p3-t6-p6-t7-p25-tl3-pl3-tl4-pl4-tlB-R0-tl-pl9-t2-p2 

21:  p4-t5-p5— t6-p6— t7-p25-t 13-p 13— 1 14— p 14  —  1 1 8-RO-t 1 -p20-t4-p4 

22:  pl-t2-p2-t3-p3-t6-p6-t7-p25-t 13-p 1 3-t 1 4-p 1 4-t 1 8-RO-t 1 -p20-t4-p4-t5-p 1 

23:  p4-t5-p5-t6-p6-t7-p25-t 13-pl 3-t 1 4-p 1 4-t 18-RO-t l-p2 1-t 8-p lB-t4-p4 

24:  p 1 -t2-p2-t3-p3-t6-p6-t7-p25-t 13-pl 3-t i 4-p 1 4-t 1 8-RO-t 1 -p21 -t8-p 1 B-t4-p4-t5-p 1 

25:  p5-t6-p6-t7-p25-t 1 3-p 13-t 1 4-p 1 4-t 1 8-RO-t 1 -p 22 - t 1 2 -p 24 - t 5- p 5 

26:  pl-t2-p2-t3-p3-t6-p6-t7-p25-t 13-p 13-t 1 4-p 1 4-t 18-RO-t l-p22-tl2-p24-t5-pl 

27:  p4-t5-p5-t6-p6-t7-p25-t 13-pl 3-t 1 4-p 1 4-t 18-RO-t l-p22-tl 2-p 1 2-t9-p9-t 1 0-p 1 0-t 11- 
R3-t8-p 1 8 - 1 4 - p  4 

28:  p 1-t 2-p2-t3-p3-t6-p6-t 7-p25-t 13-p 13-t 1 4-p 1 4-t 10-RO-t 1 -p22-t 1 2-p 1 2-t 9-p9-t 10-p 10- 
t 1 1 -R3-t B-p 1 B-t 4-p4-t5-p 1 
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TABLE  6.2. _ (Cont'd) 


p2-t3— p3-t6-p6-t7— p25-tl3-pl3— tl4-R4-tl2-pl2— t9-p9— tlO— plO-tl 1 — p 1 1— tl8-R0-t 1- 
pl9-t2-p2 

p4-t5-p5— t6-p6-t7— p25-t 13-p 1 3-t 14-R4— 1 1 2-pl2-t9-p9-t 10-p 10— 1 1 1— pi 1-tlS-RO-tl- 
p20-t4-p4 

pl-t2-p2-t3-p3-t6-p6-t7-p25-tl3-pl3-tl4-R4-tl2-pl2-t9-p9-tl0-pl0-tll-pll-tl8- 
RO-t l-p20-t4-p4-t5-pl 

p4-t5-p5-t6-p6-t7-p25-t 1 3-p 1 3-t 14— R4— t 12— p 12-t9— p9— t 1 0-p 10-t 11-pll-tlB-RO-tl- 
p21-t8-pl8-t4-p4 

p  1  -t 2-p2-t 3-p3-t 6-p6-t 7-p25-t 13-p 1 3-t 1 4-R4-t 1 2-p 1 2-t 9-p 9- 1 10-pl0-tll-pll-tl8- 
R0-tl-p21-t8-pl8-t4-p4-t5-pl 

p5-t6-p6-t7-p7-t 10-p 10-t 1 1-pl 1-t 18-RO-t l-p23-t 15-pl5-t5-p5 

p  1  -  t2— p2-t3— p3— 16— p6— t7-p7— 1 10— p 10— 1 1 1— p 1 1 -t 18-RO-t l-p23-t 15-p 1 5-t 5-p 1 

p5-t6-p6-t7-p25-tl3-pl3-tl4-pl4-tl 8-RO-t 1 -p23-t 1 5-p 15-t5-p5 

pl-t2-p2-t3-p3-t6-p6-t7-p25-tl3-pl3-tl4-p 1 4-t 18-RO-t 1 -p23— t 15-p 15-t5— p 1 

p5-t6-p6-t 7-p25-t 13-p 13-t 14-R4-t 1 2-pl2-t9-p9-t 10-p 10-t 11 -pi 1-t 18-RO-t l-p23-t 15- 
p 1 5-t5-p5 

p 1-t 2-p2-t3-p3-t6-p6-t7-p25-t 13-p 13-t 14-R4-tl2-p 12-t9-p9-t 10-p 10-t 11 -pi 1-t 18- 
RO-t  1  -p23-t  1  5-p  1  5-t  5-p  1 

p5-t6-p6-t7-p26-tl6-pl6-tl 7-R5-t 1 5-p 1 5- t 5-p 5 

pl-t2-p2-t3-p3-t6-p6-t7-p26-tl6-pl6-tl 7-R5-t 1 5-p 1 5-t 5-p 1 

p2-t3-p3-t6-p6-t7-p26-t 16-p 1 6-t 1 7-pl 7-t 18-RO-t 1 -p 1 9— t2-p2 

p4-t 5-p 5-t 6-p6-t 7-p 26-t 1 6-p 1 6-t 17-p 17-t 18-RO-t l-p20-t4-p4 

pl-t2-p2-t3-p3-t6-p6-t7-p26-tl6-pl 6-t 1 7-pl 7-t 1 8-RO-t 1 -p20-t 4-p4-t 5-p 1 

p4-t 5-p5-t6-p6-t7-p26-t 16-p 16-t 17-p 17-t 18-RO-t 1 -p 2 1 - t 8 -p 1 8- t 4-p 4 

p  1 -t 2-p2-t3-p3-t6-p6-t7-p26-t 16-p 1 6-t 1 7-pl 7-t 1 8-RO-t l-p21-tB-pl8-t4-p4-t5-pl 

p5-t6-p6-t 7-p26-tl6-pl 6-t 1 7-pl 7-t 1 8-RO-t l-p22-t 12-p24-t5-p5 

Pl-t2-p2-t3-p3-t6-p6-t7-p26-tl6-pl 6-t 17-p 17-t 18-RO-t l-p22-tl2-p24-t5-pl 

p4-t5-p5-t6-p6-t7-p26-t 1 6-p 1 6-t 1 7-pl 7-t 18-RO-t l-p22-t 1 2 -p 1 2 - t 9-p 9- t 10-p 10-t 11- 
R3-t8-p 1 8-t4-p4 

pl-t2-p2-t3-p3-t6-p6-t7-p26-t 16-p 16-t 17-p 17-t 18-RO-t l-p22-tl 2-p 1 2-t 9-p9-t 10-p 10 
tll-R3-t8-pl8-t4-p4-t5-pl 

p5-t6-p6-t7-p26-tl6-pl6-t 17-p 17-t 18-RO-t l-p23-t 15-p 1 5-t 5-p5 
Pi-t2-p2-t3-p3-t6-p6-t7-p26-t 1 6-p 1 6-t 17-p 17-t 18-RO-t l-p23-t 15-p  1 5-t 5-p  1 
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total  number  of  circuits  amounts  to  fifty-two  (52)  and  we  can  check  that 
every  circuit  oontains  at  least  one  of  the  resource  place  (i.e., 
or  R,):  we  are  sure  therefore  that  the  organizational  form 
is  admissible,  i.e.,  the  information  structure  is  acyclical,  as  discussed 
in  Chapter  4,  Section  4.1.  Finally,  the  characterization  of  the  slices  is 
presented  in  Table  6.3.  The  procedure  described  in  Chapter  5,  Section  5.1, 
has  been  used  to  determine  the  slices.  There  are  evelen  (11)  slices  that 
model  the  different  steps  of  the  complete  process.  The  transitions 
belonging  to  the  i-th  slice  correspond  to  the  tasks  that  are  performed 
concurrently  at  the  i-th  processing  step. 


We  are  now  going  to  present  the  performance  measures  obtained,  for 
different  values  of  the  resource  and  time  constraints.  In  each  case,  we 
will  determine  the  maximum  throughput  rate  and  the  execution  schedule  of 
the  organization.  Recall  that  the  task  processing  times  considered  are  all 
deterministic. 


TABLE  6.3.  LIST  OF  THE  SLICES 


Slice  Is  tl 

Slice  2:  t0  ,  1 1 2 ,  1 15 

Slice  3:  t9  ,  t4 

Slice  4i  t5 

Slice  5i  t2 

Slice  6:  t3 

Slice  7i  1 6 

Slice  8s  t7 

Slice  9s  tlO,  1 1 3 ,  1 16 

Slice  10s  til,  1 1 4 ,  1 1 7 

Slice  lit  1 1 8 


6.2.2  Measures  of  Performance 

The  results  obtained  for  each  case  considered  will  be  presented  on  a 
Table  that  contains  the  following  information: 


-  The  transition  firing  times,  corresponding  to  the  task  processing 
times,  and  the  Initial  marking  of  the  resource  places,  corresponding 
to  the  resources  available  for  processing. 

-  The  firing  schedule  obtained  for  a  total  number  of  N  repetitions  of 
the  process  which  corresponds  to  the  number  of  Inputs  processed  (the 
choice  of  N  is  arbitrary  and  can  be  as  large  as  desired). 
Accordingly,  the  firing  schedule  will  correspond  to  the  sequence  (s“), 
for  1  =  l,2,...,m  (total  number  of  transitions  in  the  net)  and 
n-1,2, . . . ,N. 

-  The  resulting  period  n  of  the  steady-state  process  and  the  periodicity 
factor  K,  as  defined  in  Chapter  5,  Section  5.2.4. 

-  The  critical  circuits  obtained,  as  defined  in  Chapter  4,  Section 
4. 2. 4.1. 

From  this  information,  the  following  measures  of  performance  will  be 
computed : 

(a)  The  maximum  throughput  rate,  defined  as  the  ratio  (Chapter  5, 
Section  5  .2  .4 ) : 

f=-  (6.1) 

n 

It  will  be  checked  that  I  is  equal  to  the  average  cirouit  processing 
rate  of  the  critical  circuits,  as  determined  by  1/a  (Chapter  4, 
Section  4.2.3) . 

(b)  The  dynamic  response  time  corresponding  to  the  complete 
processing  of  N  inputs  as  defined  in  Chapter  5,  Section  5.3.2: 


where  tm  is  the  output  transition  of  the  net.  The  ratio 


T* 

RTn  =  —  (6.3) 

will  give  an  approximate  value  of  the  average  response  time  RT 
(defined  as  the  limit  of  RTN  when  N  goes  to  infinity).  Recall  that: 

RT  =  j  (6.4) 

(c)  The  average  processing  time  (of  any  individual  input),  defined 
as: 


d  = 


d  + 
n 


n+i 


+ . . .+  d 


n+K-x 


(6.5) 


where : 


d 

n 


Sn  +  » 

i 


Recall  that  the  following  equality  holds  (Chapter  5,  Section  5.3.2): 

M# (R* ) 

d - ^ -  (6.6) 


where  M#(R#)  denotes  the  initial  marking  of  the  resource  place  of  the 
overall  organization. 
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These  measures  will  be  analyzed  and  compared  for  the  different  cases 
considered . 


6.2.3  Numerical  Results 


Four  different  cases  of  resource  and  time  constraints  will  be 
analyzed.  The  first  example  coresponds  to  the  case  where  the  DM0  as  a 
whole  and  each  of  the  five  DMs  can  only  handle  one  input  at  a  time.  All 
the  task  processing  times  are  assumed  identical  and  equal  to  one  unit  of 
time.  In  the  second  case,  there  are  more  resources  available  for 
processing:  the  DM0  can  handle  four  external  inputs  at  the  same  time  and 
each  DM  is  able  to  process  at  most  two  inputs  simultaneously.  The  task 
processing  times  remain  the  same  as  in  the  first  case  (one  unit  of  time). 
In  the  third  case,  the  resources  available  are  assumed  to  be  the  same  as  in 
the  second  example,  but  the  time  constraints  now  differ.  The  task 
processing  times  of  DM4  are  assumed  to  be  twice  as  long  as  the  other 
processing  times,  i.e.,  these  processing  times  are  equal  to  two  units  of 
time.  Finally,  the  fourth  case  is  an  example  where  the  processing  times 
are  arbitrary.  The  resources  available  are  identical  and  set  to  three  (the 
DM  and  the  organization  can  handle  three  inputs  at  a  time). 

In  each  case,  the  firing  schedule  has  been  computed  for  a  total  of  ten 
(10)  repetitions  of  the  process,  i.e.,  N=10.  Since  the  total  number  of 
transitions  in  the  net  amounts  to  m=18,  the  firing  schedule  obtained 
corresponds  to  sequence  (S°)  for  i  =  1,2,...,  18  and  n  =  1,  2,...,  10.  The 
algorithm  described  in  Appendix  C  has  been  used  to  determine  this  sequence. 
Let  us  now  present  the  results  obtained  in  each  case. 

6. 2. 3.1  First  Case 

The  results  obtained  for  the  first  case  are  shown  on  Table  6.4.  The 
period  n  is  11  units  of  time  and  the  periodicity  factor  K  equals  one.  As 
expected,  the  steady-state  process  is  strongly  periodic,  since  the  token 
content  of  the  elementary  directed  circuits  is  one  (see  Chapter  5,  Section 
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5.2.4).  According  to  (6.1).  the  maximum  throughput  rate  is: 

}  =  -jy  =  0*091  inputs  processed  per  unit  of  time 

We  can  check  that  i  =  1/a,  where  a  is  the  maximum  average  circuit 
processing  time,  which  is  also  the  average  circuit  processing  time  of  the 
critical  circuits.  The  three  critical  circuits  obtained  have  indeed  eleven 
(11)  transitions  and  the  token  oontent  is  one.  Since  the  transition  firing 
time  equal  one  unit  of  time,  we  immediately  obtain: 

a  =  11  units  of  time 

Furthermore,  all  the  critical  circuits  have  one  unique  resource  place, 
which  is  R# .  The  first  such  critical  circuit  is  represented  in  Figure  6.4. 
By  deleting  R# ,  we  obtain  the  directed  paths  (or  lines)  for  which  the 
processing  time  (called  time-delay  in  [3])  is  maximal,  as  discussed  in 
Chapter  5,  Section  5.3.2.  Accordingly,  this  time-delay  is: 

d  =  Dmax  =  11  units  of  time 

and  the  dynamic  response  time,  corresponding  to  the  processing  of  the  ten 
inputs  is  (relation  (6.2)) 

T*°  =  S*°  +  1  =  110  units  of  time. 

6. 2. 3. 2  Second  Case 

The  results  obtained  in  the  second  case  are  shown  on  Table  6.5. 
It  turns  out  that  the  steady-state  process  is  2 -periodic  (i.e., 
the  periodicity  factor  K  equals  2)  with  a  period  n  =  9  units  of  time. 
Accordingly,  the  maximum  throughput  rate  is: 

*  2 

1  =  J  =  0.222  inputs  processed  per  unit  of  time 


TABLE  6.4.  COMPPTATION  OF  THE  EXECOTION  SCHEDULE:  CASE  1 


#  RESOURCE  and  TIME  CONSTRAINTS  * 

Marking  oi  the  Resource  Places  ;  RO  R1  R2  R3  R4  R5 

111111 

Transition  Firing  Times  i 

tl  t2  t3  t4  t5  t6  1 7  tB  t9  tlO  til  1 1 2  1 1 3  1 1 4  1 1 5  tl6  1 1 7  1 18 

1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0 


*  EXECUTION  SCHEDULE  * 


tl  i 

0.0 

11.0 

22.0 

33.0 

44.0 

tB  : 

1.0 

12.0 

23.0 

34.0 

45.0 

1 1 2: 

1.0 

12.0 

23.0 

34.0 

45.0 

1 1 5: 

1.0 

12.0 

23.0 

34.0 

45.0 

1 9  i 

2.0 

13.0 

24.0 

35.0 

46.0 

1 4  : 

2.0 

13.0 

24.0 

35.0 

46.0 

t5  : 

04 

O 

14.0 

25.0 

36.0 

47.0 

t2  : 

4.0 

15.0 

26.0 

37.0 

48.0 

t3  : 

5.0 

16.0 

27.0 

38.0 

49.0 

1 6  : 

6.0 

17.0 

28.0 

39.0 

50.0 

t7  : 

7.0 

18.0 

29.0 

40.0 

51.0 

1 1 0: 

8.0 

19.0 

30.0 

41.0 

52.0 

1 1 3 : 

8.0 

19.0 

30.0 

41.0 

52.0 

1 16: 

8.0 

19.0 

30.0 

41.0 

52.0 

til: 

9.0 

20.0 

31.0 

42.0 

53.0 

1 1 4: 

9.0 

20.0 

31.0 

42.0 

53.0 

1 1 7: 

9.0 

20.0 

31.0 

42.0 

53.0 

1 1 8: 

10.0 

21.0 

32.0 

43.0 

54.0 

Period  :  11.0  Periodicity 


55.0 

66.0 

77.0 

88.0 

99.0 

56.0 

67.0 

78.0 

89.0 

100.0 

56.0 

67.0 

78.0 

89.0 

100.0 

56.0 

67.0 

78.0 

89.0 

100.0 

57.0 

6B.0 

79.0 

90.0 

101.0 

57.0 

68.0 

79.0 

90.0 

101.0 

58.0 

69.0 

80.0 

91.0 

102.0 

59.0 

70.0 

81.0 

92.0 

103.0 

60.0 

71.0 

82.0 

93.0 

104.0 

61.0 

72.0 

83.0 

94.0 

105.0 

62.0 

73.0 

84.0 

95.0 

106.0 

63.0 

74.0 

85.0 

96.0 

107.0 

63.0 

74.0 

85.0 

96.0 

107.0 

63.0 

74.0 

85.0 

96.0 

107.0 

64.0 

75.0 

86.0 

97.0 

108.0 

64.0 

75.0 

86.0 

97.0 

108.0 

64.0 

75.0 

86.0 

97.0 

108.0 

65.0 

76.0 

87.0 

98.0 

109.0 

Factor 

:  1 

Critical  Circuits  : 


Hi  pl-t2-p2-t3-p3-t6-p6-t7-p7-tl0-p JO-t 1 1 -p 1 1 -t 1 8-RO-t 1 -p21 -t 8-p 1 8-t 4-p4-t5-p 1 
24:  pl-t2-p2-t3-p3-t6-p6-t7-p25-t 1 3-p 1 3-t 1 4-p 1 4-t 1 6-RO-t 1 -p21 -t8-p 18-t4-p4-t 5-p 1 
46:  pl-t2-p2-t3-p3-t6-p6-t7-p26-t 16-pl6-t 1 7-p 1 7-t 1 8-RO-t l-p21-t 8-p 1 8-t 4-p4-t 5-p  1 


Figure  6.4  Critical  Circuit  (Case  1) 

The  process  reaches  its  steady-state  after  the  fourth  repetition 
(N=4).  For  n  >  2  and  i  =  1.2,.... 18,  we  have: 

S*n+X  =  S*n  +8  (units  of  time) 

S*n+*  =  S*n+1  +  1  (unit  of  time) 

Inputs  are  therefore  processed  at  successive  intervals  of  eight  (8)  units 
of  time  and  then  one  (1)  unit  of  time. 

In  that  case,  the  dynamic  response-  time  is  given  by: 

X10  *  s10  +  p  =  47.0  +  1.0  =  48.0  units  of  time 

it  xt  II 


It  is  interesting  to  note  that  the  ratio: 


TABLE  6.5.  COMPUTATION  OF  THE  EXECUTION  SCHEDULE:  CASE  2 


#  RESOURCE  and  TINE  CONSTRAINTS  * 

Narking  of  the  Resource  Places  t  RO  R1  R2  R3  R4  R5 

4  2  2  2  2  2 

Transition  Firing  Times  : 

tl  t2  t3  t4  t5  t6  t7  tS  1 9  tlO  til  1 12  1 1 3  1 1 4  1 15  1 16  1 1  7  1 18 

1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0 


•  EXECUTION  SCHEDULE  • 


tl  : 

0.0 

1.0 

2.0 

3.0 

11.0 

12.0 

20.0 

21.0 

29.0 

30.0 

1 8  : 

1.0 

2.0 

10.0 

11.0 

19.0 

20.0 

28.0 

29.0 

37.0 

38.0 

1 1 2: 

1.0 

2.0 

10.0 

11.0 

19.0 

20.0 

28.0 

29.0 

37.0 

38.0 

1 15s 

1.0 

2.0 

10.0 

11.0 

19.0 

20.0 

28.0 

29.0 

37.0 

38.0 

t9  : 

2.0 

3.0 

11.0 

12.0 

20.0 

21.0 

29.0 

30.0 

38.0 

39.0 

t4  i 

2.0 

3.0 

11.0 

12.0 

20.0 

21.0 

29.0 

30.0 

38.0 

39.0 

t5  : 

3.0 

4.0 

12.0 

13.0 

21.0 

22.0 

30.0 

31.0 

39.0 

40.0 

t2  : 

4.0 

5.0 

13.0 

14.0 

22.0 

23.0 

31.0 

32.0 

40.0 

41.0 

t3  : 

5.0 

6.0 

14.0 

15.0 

23.0 

24.0 

32.0 

33.0 

41.0 

42.0 

t6  : 

6.0 

7.0 

15.0 

16.0 

24.0 

25.0 

33.0 

34.0 

42.0 

43.0 

1 7  : 

7.0 

8.0 

16.0 

17.0 

25.0 

26.0 

34.0 

35.0 

43.0 

44.0 

tlO: 

8.0 

9.0 

17.0 

IB.  0 

26.0 

27.0 

35.0 

36.0 

44.0 

45.0 

1 13: 

8.0 

9.0 

17.0 

18.0 

26.0 

27.0 

35.0 

36.0 

44.0 

45.0 

1 16: 

8.0 

9.0 

17.0 

18.0 

26.0 

27.0 

35.0 

36.0 

44.0 

45.0 

til: 

9.0 

10.0 

18.0 

19.0 

27.0 

28.0 

36.0 

37.0 

45.0 

46.0 

1 14: 

9.0 

10.0 

18.0 

19.0 

27.0 

28.0 

36.0 

37.0 

45.0 

46.0 

1 1 7: 

9.0 

10.0 

18.0 

19.0 

27.0 

28.0 

36.0 

37.0 

45.0 

46.0 

1 18: 

10.0 

11.0 

19.0 

20.0 

28.0 

29.0 

37.0 

38.0 

46.0 

47.0 

Period  :  9.0 


Periodicity  Factor  :  2 


Critical  Circuits  : 

5:  pl-t2-p2-t3-p3-t6-p6-t7-p7-tl0-pl0-tll-R3-t8-pl 8-t4-p4-t 5-p 1 
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RTi#  -  — =4.8  units  or  time  per  input 


is  already  quite  close  to  the  average  response  tine  (given  by  6.4) 


RT  =  x  *  units  of  time  per  input 


The  average  response  time,  determined  by  (6.5),  is  in  that  case: 


d,  +  dxo  <<.  -  »:>  +  <T"  -  O  18  -  18 


We  can  check  that  relation  (6.6)  is  satisfied,  since: 


M°  ( R  ) 

0 


=  4  *  4.5  =  18  units  of  time 


Finally,  there  is  a  unique  critical  circuit,  which  is  represented  on 
Figure  6.5.  This  circuit  includes  nine  (9)  transitions  and  its  token 
content  is  two  (2),  corresponding  to  the  initial  marking  of  the  resource 
place  R,.  The  maximum  average  circuit  processing  time  is  therefore: 


a  ~  Y  ~  4*5  units  of  time 


which  is  equal  to  1/i,  as  expected.  It  should  be  noted  that,  in  contrast 
to  the  first  case,  the  critical  circuit  obtained  here  is  internal  to  the 
organizational  structure.  It  is  clear,  from  the  structure  of  the  circuit, 
that  the  throughput  rate  constraints  come  from  the  interactions  between 
DM1,  DM2  and  DM3.  In  addition,  only  the  resources  of  DM3  are  critical. 
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Figure  6.5  Critical  Circuit  (Case  2) 


6. 2. 3. 3  Third  Case 

The  results  obtained  in  the  third  case  are  shown  on  Table  6.6.  The 
steady-state  process  is  2-perlodlc  with  a  period  of  n  *  11  units  of  time. 
The  maximum  throughput  rate  is  therefore: 

2 

i  =  yj-  =  0.182  inputs  processed  per  unit  of  time 

The  process  reaches  its  steady-state  after  the  fifth  repetition  (N=5).  For 
n  >  3  and  i  *  1,2,... ,18,  we  have: 


Js*n  *  S*n  +  2  (units  of  .time) 

(s*n+1  =  S*n  +  7  (units  of  time) 

The  dynamic  response  time  is,  in  that  case: 

T*°  =  Si;  +  pig  *  58.0  +  1.0  =  59.0  units  of  time 


TABLE  6.6.  COMPUTATION  OF  THE  EXBCOTION  SCHEDOLE:  CASE  3 


«  RESOURCE  and  TIME  CONSTRAINTS  * 

Narking  of  the  Resource  Places  i  RO  Ri  R2  R3  R4  R5 

4  2  2  2  2  2 

Transition  Firinq  Tiees  : 

tl  t2  t 3  t4  t 5  t6  t7  tS  t9  tlO  til  t 12  t 13  t 14  t 1 5  t 16  t 1 7  t 18 

1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0  2.0  2.0  2.0  1.0  1.0  1.0  1.0 


*  EXECUTION  SCHEDULE  • 


tl  : 

0.0 

1.0 

2.0 

3.0 

13.0 

15.0 

24.0 

26.0 

35.0 

37.0 

te  t 

1.0 

2.0 

10.0 

12.0 

21.0 

23.0 

32.0 

34.0 

43.0 

45.0 

1 12: 

1.0 

3.0 

12.0 

14.0 

23.0 

25.0 

34.0 

36.0 

45.0 

47.0 

1 15: 

1.0 

2.0 

10.0 

12.0 

21.0 

23.0 

32.0 

34.0 

43.0 

45.0 

1 9  : 

3.0 

5.0 

14.0 

16.0 

25.0 

27.0 

36.0 

38.0 

47.0 

49.0 

t4  : 

2.0 

3.0 

11.0 

13.0 

22.0 

24.0 

33.0 

35.0 

44.0 

46.0 

t5  : 

3.0 

5.0 

14.0 

16.0 

25.0 

27.0 

36.0 

38.0 

47.0 

49.0 

t2  : 

4.0 

6.0 

15.0 

17.0 

26.0 

28.0 

37.0 

39.0 

48.0 

50.0 

t3  : 

5.0 

7.0 

16.0 

18.0 

27.0 

29.0 

38.0 

40.0 

49.0 

51.0 

1 6  i 

6.0 

8.0 

17.0 

19.0 

28.0 

30.0 

39.0 

41.0 

50.0 

52.0 

1 7  : 

7.0 

9.0 

18.0 

20.0 

29.0 

31.0 

40.0 

42.0 

51.0 

53.0 

1 10: 

8.0 

10.0 

19.0 

21.0 

30.0 

32.0 

41.0 

43.0 

52.0 

54.0 

1 1 3: 

8.0 

10.0 

19.0 

21.0 

30.0 

32.0 

41.0 

43.0 

52.0 

54.0 

1 16: 

8.0 

10.0 

19.0 

21.0 

30.0 

32.0 

41.0 

43.0 

52.0 

54.0 

til: 

9.0 

11.0 

20.0 

22.0 

31.0 

33.0 

42.0 

44.0 

53.0 

55.0 

1 1 4  * 

10.0 

12.0 

21.0 

23.0 

32.0 

34.0 

43.0 

45.0 

54.0 

56.0 

1 1 7: 

9.0 

11.0 

20.0 

22.0 

31.0 

33.0 

42.0 

44.0 

53.0 

55.0 

1 18: 

12.0 

14.0 

23.0 

25.0 

34.0 

36.0 

45.0 

47.0 

56.0 

58.0 

Period  :  11.0 


Periodicity  Factor  :  2 


Critical  Circuits  : 

1 7 i  pl-t2-p2-t3-p3-t6-p6-t7-p25-tl3-pl3-tl4-R4-tl2-p24-t5-pl 


V\“-'  '  •  *• 


and  the  average  processing  tine: 


(T*  -  S*>  +  (t; 


d  +  d 


6. 2. 3. 4  Fourth  Case 


The  results  obtained  in  the  fourth  case  are  shown  in  Table  6.7.  In 
this  example,  we  have  taken  totally  arbitrary  transition  firing  times  and 
set  the  initial  marking  of  all  the  resouroe  places  to  three  (3)  tokens.  It 
turns  out  that  the  steady-state  process  is  3-perlodic  with  a  period  w-25 .3 
units  of  time.  Accordingly,  that  maximum  throughput  rate  is: 


i 


3 

23.5 


0.118 


inputs  processed  per  units  of  time 


We  note  that  the  process  reaches  its  steady-state  after  the  fourth 
repetition.  For  n  2  2  and  1  =  1,2,..., 18,  we  obtain: 


S*n  1  =  S*n_1  +3.5  (units  of  time) 
S*n  =  S*n_1  +3.5  (units  of  time) 
S*n+‘  «  S*n  +  18.5  (units  of  time) 


The  dynamic  response  time  is  in  that  case : 


T*°  *  Sl#  +u  »  100.0  +  2.0  =  102.0  units  of  time 

it  i»  'll 


and  the  average  processing  time: 


d  +  d  +  d 

•  *  i* 


(T*  -  S*)  +  (T* 


-  S*)  +  (T1#  -  S10) 


TABLE  6.7.  COMPUTATION  OF  THE  EXECUTION  SCHEDULE:  CASE  4 


»  RESOURCE  and  TIME  CONSTRAINTS  * 

Markina  of  the  Resource  Places  :  RO  R1  R2  R3  R4  R5 

3  3  3  3  3  3 

Transition  Firing  Tines  : 

tl  t2  t3  t4  t 5  t 6  t7  t 8  t9  tlO  til  t 1 2  t 13  t 14  t 15  t 16  t 1 7  t 1 8 

2.0  1.5  3.0  3.5  2.5  3.0  2.0  1.5  3.0  2.0  2.5  1.5  1.0  2.0  3.5  2.0  1.0  2.0 


*  EXECUTION  SCHEDULE  • 


tl 

0.0 

2.0 

4.0 

25.5 

29.0 

32.5 

51.0 

54.5 

58.0 

76.5 

t8 

2.0 

4.0 

6.0 

27.5 

31.0 

34.5 

53.0 

56.5 

60.0 

78.5 

1 1 2 

2.0 

4.0 

6.0 

27.5 

31.0 

34.5 

53.0 

56.5 

60.0 

78.5 

1 1 5 

2.0 

5.5 

9.0 

27.5 

31.0 

34.5 

53.0 

56.5 

60.0 

78.5 

1 9 

3.5 

6.5 

9.5 

29.0 

32.5 

36.0 

54.5 

58.0 

61.5 

80.0 

1 4 

3.5 

7.0 

10.5 

29.0 

32.5 

36.0 

54.5 

58.0 

61.5 

80.0 

t 5  :  7.0 

10.5 

14.0 

32.5 

36.0 

39.5 

58.0 

61.5 

65.0 

83.5 

t2 

9.5 

13.0 

16.5 

35.0 

38.5 

42.0 

60.5 

64.0 

67.5 

86.0 

t3 

11.0 

14.5 

18.0 

36.5 

40.0 

43.5 

62.0 

65.5 

69.0 

87.5 

1 6 

14.0 

17.5 

21.0 

39.5 

43.0 

46.5 

65.0 

68.5 

72.0 

90.5 

1 7 

17.0 

20.5 

24.0 

42.5 

46.0 

49.5 

68.0 

71.5 

75.0 

93.5 

tlO 

19.0 

22.5 

26.0 

44.5 

48.0 

51.5 

70.0 

73.5 

77.0 

95.5 

1 1 3 

19.0 

22.5 

26.0 

44.5 

48.0 

51.5 

70.0 

73.5 

77.0 

95.5 

1 1  6 

19.0 

22.5 

26.0 

44.5 

48.0 

51.5 

70.0 

73.5 

77.0 

95.5 

til 

21.0 

24.5 

28.0 

46.5 

50.0 

53.5 

72.0 

75.5 

79.0 

97.5 

1 1 4 

20.0 

23.5 

27.0 

45.5 

49.0 

52.5 

71.0 

74.5 

78.0 

96.5 

1 1  7 

21 . 0 

24.5 

28.0 

46.5 

50.0 

53.5 

72.0 

75.5 

79.0 

97.5 

1 1  8 

23.5 

27.0 

30.5 

49.0 

52.5 

56.0 

74.5 

78.0 

81.5 

100.0 

Period  :  25.5  Periodicity  Factor  :  3 

Critical  Circuits  : 

11:  p 1 -t2-p2-t 3-p3-t6-p6-t7-p7-t 10-pl0-tll-pl 1  - 1 1 8 -RO- 1 1 -p 2 1 - 1 8-p 1 8- 1 4- p 4 - 1 5- p 1 
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As  expected  (relation  (6.6): 


M*(R#) 

nr 


25.5  -  d 


There  is  also  a  unique  oritical  circuit,  which  is  the  same  as  the  first 
such  circuit  obtained  in  Case  1  (see  Figure  6.4).  The  maximum  average 
circuit  processing  time  is: 

U  +  11  +  U  +11  +  U  +U  +  11  +  11  +11  +11  +11 

*1  ^  1  Hl«  'll  *11  *1  'l  *4  K* 

a  - - 

M#(R  ) 

0 

__  1.5  +  3.0  +  3.0  +  2.0  +  2.0  +  2.5  +  2.0  +  2.0  +  1.5  +  3.5  +  2.5 

3 

»  =  8.5  units  of  time 

As  expected,  we  obtained  $  -  1/a. 

6.2.4  Summary  Analysis  of  the  Results 


6.2 .4.1  Comparison  of  the  Performance  Measures 

The  performance  measures  obtained  in  each  of  the  four  cases  analyzed 
previously  are  sunmarized  in  Table  6.8.  We  note  first  that  the  performance 
obtained  in  the  first  case  is  the  worst  one,  compared  to  the  other  oases. 
The  maximum  throughput  rate  f  is  the  lowest  and  the  dynamic  response  time 
T^  (reap,  the  average  response  time  RT)  is  the  longest.  In  the  seoond 
case,  where  there  are  more  resources  available  for  processing,  the 
performance  is  greatly  improved.  The  maximum  throughput  rate  is  more  than 
twice  the  one  obtained  in  the  first  case  while  the  response  time  is  less. 
However,  the  average  processing  time,  d,  which  represents  the  average 
amount  of  time  it  takes  to  process  any  individual  input,  is  now  larger,  as 
emphasized  in  Chapter  5,  Section  5.3.2.  In  the  third  case,  the  resources 


a 


available  for  processing  are  the  same  as  In  the  seoond  case,  but  the 
processing  times  of  DM4  are  tvlce  the  other  ones.  As  expected,  the 
performance  obtained  Is  slightly  degraded,  when  compared  to  the  second 
case,  but  Is  still  much  better  than  In  the  first  case.  Recall  that  In  the 
fourth  case,  the  task  processing  times  are  arbitrary  (but  are  greater  or 
equal  to  one  unit  of  time)  and  there  are  more  resources  than  In  the  first 
three  cases.  It  turns  out  that  the  performance  obtained  Is  worse  than  In 
case  2  and  3,  but  Is  still  better  than  In  oase  1. 


6 .2 .4 .2  Execution  Schedule  of  the  Concurrent  Tasks 

In  the  previous  section;  we  have  been  able  to  evaluate  and  oompare  the 
performance  of  the  DHO  by  computing  different  time-related  measures  from 
the  execution  schedule  obtained.  The  performance  measures  were  computed, 
given  the  time  and  resource  constraints  set  in  each  oase.  Ve  want  to 
analyze  further  how  the  constraints  can  be  modified  so  as  to  improve  the 
performance  of  the  organization.  In  Chapter  4,  this  problem  has  already 
been  addressed  for  the  maximum  throughput  rate.  It  was  said  that  the 
critical  circuits,  such  as  those  obtained  in  the  four  cases,  determine 
precisely  which  of  the  time  and  resource  constraints  bound  the  throughput 
rate.  We  want  to  investigate  now  how  the  time-delays  in  the  organization 
are  related  to  the  various  task  processing  times  and,  in  particular,  how  it 
is  possible  to  improve  the  dynamic  response  time. 

At  this  point,  the  execution  schedule  obtained  for  the  concurrent 
tasks,  i.e.,  the  instants  of  time  at  which  the  transitions  within  a  slice 
fire,  is  very  interesting  to  analyze.  In  the  four  cases  described 
previously,  we  can  note  that,  for  cases  1  and  2,  all  the  transitions  within 
a  slice  fire  at  the  same  instants  of  time  (at  each  repetition  of  the 
process),  which  is  not  true  for  case  3  and  4  (recall  that  the  list  of  the 
slices  is  provided  on  Table  6.3).  This  result  is  not  surprising:  as  stated 
in  Chapter  5,  Section  5.1.2,  when  all  the  transition  firing  times  are  equal 
to  one  unit  of  time,  the  slice  T^  (i  =  l,2,...,s)  determines  the 
transitions  that  fire  concurrently  at  the  i-th  instant.  However,  when  the 
transition  firing  times  are  arbitrary,  there  is  no  reason  why  the 
concurrent  operations  should  occur  at  the  same  instant  of  time.  This  is 
actually  what  makes  the  process  to  occur  asynchronously  in  real-time. 

From  the  execution  schedule  it  is  possible  to  identify  which 
operations  take  place  (repeatedly)  at  the  latest  instants,  oompared  to  the 
other  concurrent  operations  (i.e.,  compared  to  the  other  transitions  of  the 
same  slice).  For  instance,  in  case  3,  the  transitions  that  fire  the 
latest  are  tia  (compared  to  t#  and  t1Jt  which  are  the  other  transitions  of 


T,)  and  t14  (compared  to  t14  and  t1?,  which  are  the  other  transitions  of 
T1#).  It  turns  out  that  these  particular  transitions  are  critical  for  the 
processing  delays.  We  are  sure  that,  by  reducing  their  processing  times, 
the  time-delay  of  the  organization  will  improve.  By  determining,  from  the 
execution  schedule  obtained,  the  concurrent  tasks  that  fire  at  the  latest 
instants,  we  therefore  know  which  time  constraints  are  critical  for  the 
processing  delays.  It  becomes  easier  to  modify  the  right  constraints  that 
can  improve  the  related  performance  of  the  organization. 
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CHAPTER  VII 


CONCLUSIONS  AND  DIRECTIONS  FOR  FURTHER  RESEARCH 


7.1  CONCLUSIONS 

In  this  thesis,  the  time-related  performances  of  a  DMO  have  been 
analyzed  and  evaluated.  The  computation  of  two  measures  of  performance  has 
been  carried  out:  the  maximum  throughput  rate,  which  characterizes  the 
maximum  processing  rate  achievable  by  the  organization,  and  the  execution 
schedule,  determining  the  earliest  instants  at  which  the  various  tasks  can 
be  executed.  In  order  to  evaluate  these  MOPs,  the  DMO  has  been  modeled  as 
an  asynchronous  concurrent  system  using  Timed  Petri  Nets. 

The  modeling  addresses  two  types  of  contraints.  On  the  one  hand, 
there  are  the  time  constraints  derived  from  the  various  task  processing 
times.  On  the  other  hand,  the  resource  constraints,  which  affect  both 
individual  DMs  and  the  DMO  as  a  whole.  For  the  human  DMs  in  the 
organization,  the  resource  limitation  corresponds  to  the  limited  capacity 
of  the  human  short-term  memory,  which  bounds  the  amount  of  Information  that 
can  be  handled  at  the  same  time.  For  the  overall  organization,  the 
resource  constraint  derives  from  the  limited  capacity  of  the  system  to 
store  the  information  currently  processed.  In  both  cases,  the  resource 
limitation  is  modeled  as  a  constraint  on  the  total  number  of  inputs  that 
can  be  processed  simultaneously. 

The  maximum  throughput  rate  is  expressed  as  a  function  of  the  resource 
and  time  constraints  in  the  following  manner.  The  inclusion  of  the 
resource  constraints  in  the  Petri  Net  modeling  the  DMO  results  in  directed 
circuits  (or  loops),  which  are  characterized  by: 

-  the  circuit  processing  time,  defined  as  the  sum  of  the  different 

task  processing  times  of  the  circuit.  ^  represents  the  amount  of  time 


it  takes  for  one  input  to  complete  the  processing  operations  of  the 
circuit. 

-  The  amount  of  resources  available,  n,  which  bounds  the  total  number  of 
inputs  that  can  be  processed  at  the  same  time  in  the  circuit. 

For  a  given  circuit,  the  ratio  n/p  characterizes  the  average  circuit 
processing  rate.  As  shown  in  Chapter  4,  the  minimum  average  oircuit 
processing  rate,  taken  over  all  the  directed  circuits  of  the  net, 
determines  the  maximum  throughput  rate  of  the  deterministic  system,  i.e., 
when  all  the  task  processing  times  are  deterministic.  For  the  non- 
deterministic  system,  where  processing  times  are  discrete  random  variables 
with  a  certain  probability  distribution  (to  model  decision  switches),  we 
only  obtain  an  upper  bound  to  the  maximum  throughput  rate:  In  that  case, 
the  average  circuit  processing  time  is  computed,  taking  the  expected  (or 
mean)  processing  time  of  each  task.  The  determination  of  the  critical 
circuits,  for  which  the  corresponding  average  processing  rate  is  precisely 
minimal,  is  very  useful  for  evaluating  and  comparing  different 
organizations.  Indeed,  the  critical  circuits  characterize  the  particular 
time  and  resource  constraints  that  actually  bound  the  throughput  rate.  The 
way  with  which  the  different  constraints  affect  the  performance  of  the 
organization  can  now  be  completely  specified.  In  particular,  the  problem 
of  modifying  the  right  constraints  so  as  to  improve  the  performance  of  the 
actual  design  becomes  transparent. 

In  Chapter  5,  we  have  developed  a  method  for  obtaining  and  analyzing 
the  exact  execution  schedule  of  the  deterministic  system.  In  particular, 
we  have  described  a  representation,  as  defined  by  the  slices  of  the  net, 
that  allows  for  a  clear  characterization  of  the  causality  of  the  system. 
The  causality  refers  to  the  partial  ordering  of  the  different  operations 
that  take  place  in  the  process.  The  execution  schedule  so  obtained 
determines  the  earliest  instants  at  which  the  various  tasks  can  be  executed 
in  real  time  for  a  process  that  occurs  repetitively.  An  interesting  aspect 
of  the  steady-state  process  has  also  been  highlighted,  namely  the  condition 


of  K-periodicity,  according  to  the  results  obtained  by  Chretienne  [16]. 

Once  we  have  determined  the  maximum  rate  with  which  Inputs  can  arrive 
without  overloading  the  system,  as  specified  by  the  maximum  throughput 
rate,  the  execution  schedule  becomes  useful  for  evaluating  the  time-delays 
associated  with  the  various  processing  tasks.  Two  useful  MOPs  can  be 
computed  from  the  execution  schedule:  (a)  the  earliest  instants  at  which 
the  processing  of  the  incoming  inputs  can  start,  and  (b)  the  dynamic 
response  time  of  the  organization,  characterized  as  the  minimum  amount  of 
time  it  takes  to  complete  the  processing  of  any  large  number  of  external 
inputs.  In  our  modeling,  both  the  DMO  as  a  whole  and  individual  DMs  are 
allowed  to  handle  more  than  one  input  at  a  time,  depending  on  the  type  of 
resource  constraints.  In  such  cases,  the  response  time  necessarily  depends 
on  the  number  of  inputs  being  processed  simultaneously. 

Finally,  the  evaluation  of  the  maximum  throughput  rate,  together  with 
the  precise  execution  schedule,  completely  characterize  the  dynamic 
behavior  of  the  system.  These  measures  describe  entirely  the  time-related 
performance  of  the  DMO. 

7.2  DIRECTIONS  FOR  FURTHER  RESEARCH 

The  evaluation  of  the  maximum  throughput  rate  and  of  the  execution 
schedule  has  been  obtained,  in  this  work,  for  deterministic  systems,  where 
all  the  task  processing  times  are  assumed  deterministic.  However,  as 
discussed  in  Chapter  3,  the  Petri  Net  representation  of  the  DMO  also 
includes  switches,  which  model  the  particular  processing  stages  where 
decisions  take  place.  Because  the  algorithms  selected  by  the  decisions  may 
have  in  reality  different  processing  times,  the  processing  time  of  these 
particular  tasks  was  treated  as  a  discrete  random  variable,  with  a 
probability  distribution  that  depends  on  the  decision  rules.  The 
probability  distribution  was  assumed,  however,  to  be  independent  from  one 
input  to  the  next,  since  tokens  (which  model  the  inputs  in  the  Petri  Net 
framework)  are  indistinguishable  in  ordinary  Petri  Nets.  In  such  a  case. 


i.e.,  when  the  system  is  non-deterministic,  we  have  only  obtained  a  lower 
and  an  upper  bound  for  the  performance  measures. 

Further  research  is  needed  to  obtain  more  accurate  measures  of 
performance,  by  relaxing  some  of  the  assumptions  made  about  the  decision 
switches.  In  particular,  it  appears  promising  to  extend  the  modeling  of 
the  DMO  using  a  more  sophisticated  type  of  Petri  Net,  called  Colored  Petri 
Nets  [23].  In  contrast  to  ordinary  Petri  Nets,  the  tokens,  in  this  type  of 
Petri  Nets,  are  marked.  More  precisely,  each  token  is  identified  by  a 
certain  color,  so  that  it  is  now  possible  to  differentiate  the  tokens. 
Clearly  enough.  Colored  Petri  Nets  provide  a  useful  tool  to  account  for  the 
different  type  of  inputs  that  an  actual  organization  receives  and 
processes.  As  developed  in  [1],  external  inputs  can  be  treated  as  elements 
cf  an  alphabet  and  the  decision  that  rules  the  selection  of  a  particular 
algorithm  depends  now  on  the  specific  element  in  the  alphabet.  Therefore, 
we  can  easily  foresee  how  Colored  Petri  Nets  would  be  relevant  for 
evaluating  more  accurately  the  time  performance  of  a  DMO:  In  that  case, 
the  probability  distribution  that  rules  the  selection  among  the  possible 
algorithms  (and  consequently  the  corresponding  processing  time),  should 
depend  on  the  color  of  the  token,  rather  than  being  independent  from  one 
token  to  the  next.  The  Petri  Net  model  so  analyzed,  would  now  fully 
account  for  the  assumptions  that  are  made  in  the  information  theoretic 
framework . 
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APPENDIX  A 


ALGORITHM  TO  OBTAIN  ALL  THE  CIRCUITS  OF  AN  EVENT-GRAPH 


A.l  SUMMARY 


This  appendix  presents  the  algorithm  that  has  been  used  to  obtain  all 
the  circuits  of  an  Event-Graph.  As  discussed  In  Chapter  4,  It  Is  necessary 
to  determine  all  the  circuits  of  the  net,  so  as  to  oompute  the  maxi  mum 
throughput  rate.  The  algorithm  described  here  has  been  developed  by 
Alalwan  and  Toudlc  [19]  and  determines,  in  fact,  all  the  minimal  support 
S-invariants  of  an  ordinary  Petri  Net.  As  we  have  proved  (Theorem  2.6), 
the  directed  circuits  of  an  Event-Graph  correspond  to  the  minimal 
S-compooent  of  the  net.  This  Is  why  this  algorithm  can  be  used  directly. 


A. 2  ALGORITHM  TO  OBTAIN  ALL  THE  MINIMAL  SUPPORT  S-INVARIANTS  OF  A  GENERAL 


PETRI  NET 


The  complete  proof  of  the  algorithm  can  be  found  In  Martinez  and  Silva 
[13],  An  Improved  and  more  detailed  version  is  given  in  [19].  He  first 
present  the  description  of  the  algorithm  and  then  roughly  explain  how  it 


works . 


In  the  following  description,  C  is  assumed  to  be  the  Incidence  Matrix 


of  the  Petri  Net,  of  dimensions  n  x  m,  where  n  is  the  number  of  places  and 


m  the  number  of  transitions. 


C^j.  denotes,  as  usual,  the  element 


corresponding  to  the  ith  row  and  j-th  column  of  C  and  IQ  denotes  the  n  x  n 
identity  matrix. 


Description  of  the  Algorithm 


Step  1: 


C  ;  D 


(initialization} 


Step  2:  REPEAT  for  J  =  1  UNTIL  j 


(main  loop} 


L  *>  m  h  •  m  m  *  ■»  F  *  ■  m  ^  • *  —  -  ■  •  .  *  .  *  ■*  *  <*  »  *  «  *  .  J  *  "  »  *  *  *  ►  *  *•“»** 


Step  2,1:  Determine  the  sets  Ix  and  Ia  such  that: 


Ii  =  {1x  7  Ai  i  >  0}  and  *  {ii  7  Ai  i  <  01 
AxJ  ^xJ 


Step  2.2:  For  all  pairs  (ix,  ia)  s  Ix  x  Ia  do 


Step  2.2.1:  Append  to  the  Matrix  A  the  row  vector: 


A.  .  •  (i4-th  row  of  A)  -  A.  *  (ix-th  row  of  A) 


Step  2.2.2:  Append  to  the  Matrix  D  the  row  vector: 


L  .  *  (ia-th  row  of  D)  -  A.  •  (ia-th  row  of  D) 


Step  2.3:  Eliminate  from  A  and  D  all  the  rows  with  index  i  >  IjUl*. 

Step  2.4:  Eliminate  from  D  all  the  rows  whose  support  is  non-minimal  with 
respect  to  the  other  rows  of  D.  If  two  rows  have  the  same 
support  eliminate  one  of  them. 

Step  2.5:  Eliminate  from  A  the  rows  corresponding  to  the  rows  eliminated 
from  D. 


Step  3:  The  rows  of  D  determine  all  the  minimal  support  S-invarlant  of 

the  net. 


Remark:  In  the  algorithm,  we  call  support  of  a  n-posltlve  Integer  vector 
v,  the  set  of  elements,  i  e  (l,2,...,n),  such  that  the  1-tb  component 
of  v  is  non-null .  If  G  *  {vx,...,vk)  is  a  set  of  veotors,  vA 
(i  e  {l,2,...k})  will  be  a  vector  whose  support  is  minimal  in  G  iff  there 
does  not  exist  in  G  a  non-null  vector  Vj  such  that  its  support  is  strictly 
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included  In  the  support  of  vA 

Let  us  make  a  brief  oomment  on  how  the  algorithm  works.  The 
underlying  idea  is  to  proceed  in  m  steps  (where  m  is  the  number  of 
transitions),  by  finding  the  minimal  support  S-invariants  of  the  R*  net 
(i«=0,l, . . . ,m-l) ,  which  results  from  the  elimination  of  the  transitions 
(i  +  l,...,m).  The  Matrix  D  is  constructed  so  that  at  the  i-th  iteration, 
the  row  vectors  of  D  are  precisely  the  minimal  support  S-invariants  of  the 
net  R*.  Initially  D  =  IQ  because  R®  corresponds  to  the  net  without 
transitions  and  therefore  each  single  place  constitutes  a  minimal  support 
invariant.  Once  the  minimal  support  S-invariants  of  the  net  R*  are 
obtained,  those  of  the  net  R*+1  are  generated  in  the  following  manner:  As 
they  obviously  are  also  invariants  for  the  net  R*.  they  can  be  expressed  as 
a  positive  linear  combination  of  the  minimal  support  S-invariants  of  R* 
(Chapter  2.  Section  2.2.  Theorem  2.2),  l.e.,  as  the  rows  of  D.  However, 
the  invariants  generated  are  not  all  minimal  support  S-invariants  and  it  is 
necessary  to  eliminate  the  ones  that  are  not  by  comparing  their  supports. 

A. 3  DETERMINATION  OF  THE  DIRECTED  CIRCUITS  FOR  AN  EVENT-GRAPH 

We  have  proved  that  the  directed  circuits  of  an  Event-Graph  are  the 
minimal  S-components  (Theorem  2.6).  Let  us  recall  that  an  S-component  is  a 
subnet  constructed  as  follows: 

-  the  set  of  places  P_  is  the  support  of  the  corresponding  minimal 

9 

S-in variant, 

-  the  set  of  transitions  T_  are  all  the  transitions  of  the  PN 

9 

connected  td  the  places  of  Pg,  i.e., 

T  -  U  (p-U'  P> 

3  peP 
r  s 


Once  we  have  obtained  all  the  minimal  support  S-invariants  of  the  Event- 


Graph,  using  the  algorithm  of  Alawaln  and  Toudic  [19],  the  determination  of 
the  circuits  becomes  straightforward:  For  each  minimal  support,  we  then 
determine  the  unique  output  (or  input)  transition  of  each  place  of  the 
support  and  we  immediately  obtain  the  corresponding  S-oomponent,  which  is  a 
circuit. 

A. 4  EXAMPLE 

Let  us  illustrates  the  algorithm  with  a  simple  example.  The  example 
considered  is  shown  on  Figure  A.l. 


Figure  A.l  Example  of  an  Event-Graph 


The  Incidence  matrix  is: 


t,  t, 


.  e  4  -  , 


Let  us  now  carry  out  the  steps  of  the  algorithm  described  previously 


1 

I 


<») 

1*1 

( 

!*$•! 
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Initially  (Step  1): 


A,D  * 


For  j  =  1 

Step  2.1;  Ix  -  (2}  I,  =  {1} 

Therefore,  we  append  the  row:  row  1  +  row  2. 
Step  2.2: 


At  the  next  step,  we  eliminate  row  1  and  row  2: 


Step  2.3 : 


0  1-1 


0-11 


0  0  10 


A,0  ■  0  -1  1  .  0  0  0  1 
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•f 

9  ,'L 
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The  rows  of  D  are  obviously  minimal  support  so  no  elimination  occurs 
at  Step  2.4  and  Step  2.5. 
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Step  2.1 :  Ix  =  11}  I,  =  12(3} 

Step  2.2  and  2.3:  we  append  the  two  rows  formed  by  (row  1  +  row  2) 
and  (row  1  +  row  3)  and  eliminate  afterwards  row  1,  row  2,  and  row  3, 
so  that  the  new  matrices  are: 


0  °l.f0  0  1  M 

Lo  o  oj  Li  i  i  oj 


Once  again,  the  two  rows  of  D  are  minimal  supports.  Now  the  last 
iteration,  corresponding  to  j  =  3  is  useless,  since  the  third  column  of  A 
is  already  null.  The  two  rows  of  D  are  therefore  the  minimal  support  S- 
invariants  of  the  net,  which  correspond  to  the  set  of  places:  (pt,  p4}  and 
(px,  Pa»  p,).  The  two  directed  circuits  are  immediately  obtained: 

Px  -  <P,  t#  p4  ta  pa) 

Pa  =  (pa  fca  Pa  fca  P,  fc,  Pa} 


APPENDIX  B 


DEMONSTRATION  OF  THEOREM  5.1 


In  this  appendix,  we  present  the  proof  of  Theorem  5.1,  stated  in 
Chapter  5,  Section  5.2.2.  The  result  was  obtained  by  Ramchandani  [6],  but 
the  approach  taken  here  Is  a  different  version  of  the  proof. 


Let  us  recall  the  result  stated  in  Theorem  5.1:  A  firing  schedule 
(S^)  is  feasible  iff,  for  any  pair  of  transitions  (t^,tj)  connected  by 
place  pAj  and  for  n  =  1,2,3,...: 


(B.l) 


where : 


S"  denotes  the  instant  of  the  n-th  firing  initiation  of  transition  ti 
Ml  is  the  firing  of  time  of 

M^j  is  the  initial  marking  (at  v=0)  of  place  p±j  (i.e.,  the  unique 
place  whose  input  transition  is  tj  and  output  transition  t j ) . 

Let  us  recall  the  notation  introduced  in  Chapter  2,  Section  2.4.1: 

Ij(t)  denotes  the  number  of  Initiations  of  transition  tj  during  the 
interval  of  time  [0,t]. 

Tj(t)  denotes  the  number  of  terminations  of  transition  tj  during  the 
interval  of  time  [0,t], 

The  marking  of  place  p^j  at  any  Instant  t  is  given  by  (Equation  2.10): 
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v*>  -  Mu  +  VT)  -  V*’ 


(B.2) 


From  what  has  been  discussed  in  Chapter  5,  Section  5.2.1,  it  should  be 
clear  that  a  firing  schedule  is  feasible  if  and  only  if  the  marking  of 
each  place  remains  non- negative  at  any  Instant  of  time  (assuming  that  if  a 
firing  should  be  initiated  before  the  transition  is  enabled,  the  marking  of 
the  input  places  becomes  strictly  negative).  In  order  to  prove  the 
theorem,  we  need  therefore  to  prove  the  equivalence  between  (B.l)  and  the 
fact  that,  at  any  instant  x,  1  0. 


Proof : 


(1)  Suppose  that  the  firing  schedule  is  feasible  and  let  us  consider  the 


instant  t,  such  that: 


t  =  S 


n+Mij 


where  n  is  any  positive  integer.  This  is  the  instant  of  the 


(n+MV)-th  firing  initiation  of  transition  tj.  Accordingly: 


I  ,(t)  =  n  +  M, . 
j  ij 


Hence: 


V'> 


Mu  +  V*’ 


n-  M.  .  =  T, (t)-n 
ij  i 


from  (B.2). 


By  assumption  M^U)  1  0  (since  the  firing  schedule  is  feasible)  and 


therefore:  T^(x)  1  n.  This  implies  that,  during  the  interval  of  time 
10, t],  at  least  n  firing  terminations  of  transition  tA  have  occured. 


V. 


we 


Since  S^+(i^  is  the  instant  of  the  n-th  firing  termination  of  t^, 
deduce : 


S1  *  “i 


(2)  Let  us  assume  now  that  the  set  of  inequalities  described  by  (B.l)  are 
verified.  If  we  consider  any  instant  t  such  that: 

n+M® ,  ,  n+i+M! , 

Sj  <  .  <  Sj  « 

where  n  is  any  positive  integer,  then  the  number  of  initiation  firings 
of  transition  tj  during  the  interval  [0,t]  is  precisely: 

I ,(t)  =  n  +  M* 


By  assumption,  we  have  also: 


x  2  s”  +  (i  (from  B.l) 

which  implies  that,  during  [0,t],  there  have  been  at  least  n  firing 
terminations  of  transition  tA.  Hence: 

T^t)  >  n 

Applying  equation  (B.2)  yields: 

Vl>  *  MIj +  V” '  Mlj  2  0 

If  we  consider  now  any  instant  x  such  that 
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^  'r:  v —  r"J 


,  "ii 
0  <  x  <  S 

J 


It  means  that,  during  the  interval  £0,t],  there  have  been  leas  than 
firing  initiations  of  transition  tj.  Accordingly: 


VTl  - 


and  we  deduce  immediately  that: 


Mij(x)  =  Mij  +  Ti(t)  ~  Ij(‘c)  “  ti(t)  ~  °* 


The  firing  schedule  is  therefore  feasible. 


Q.E.D. 


APPENDIX  C 


ALGORITHM  TO  COMPOTE  THE  FIRING  SCHEDULE 


In  this  Appendix,  we  present  the  algorithm  that  has  been  developed  to 
obtain  the  execution  schedule  of  the  DMO,  according  to  the  analysis  carried 
out  in  Chapter  5.  We  first  recall  the  key  points  of  the  analysis  and  then 
describe  the  corresponding  algorithm. 

C.l  SUMMARY  ANALYSIS  OF  THE  FIRING  SCHEDULE 

The  execution  schedule  that  is  determined  here  characterizes  the 
earliest  instants  at  which  the  various  tasks  can  be  executed,  assuming  that 
the  processing  starts  at  t-0  and  occurs  repetitively.  We  denote  by: 

^  the  (deterministic)  firing  time  of  tA 

S"  the  instant  of  the  n-th  firing  initiation  of  transition  t ^ 

(corresponding  to  the  n-th  repetition  of  the  process) 

Pjj  the  unique  place  whose  input  transition  is  t^  and  output 
transition  tj 

M^j  the  initial  marking  (at  t=0)  of  the  corresponding  place  p^j 

Recall  that  M°j  =  0,  except  for  the  resource  places  of  the  organization,  in 
which  case  the  initial  marking  corresponds  to  the  resources  available  for 
processing.  This  assumption  means  that  at  x=0,  there  are  no  inputs  already 
being  processed. 

For  each  transition  tj,  we  consider  now  all  its  input  plaoes,  that  we 
denote  by  j,  pi  j , . . . ,  p^  j  and  the  input  transitions  of  each  of  these 
places,  i.e.,  t,  ,  t,  ,...,t,  ,  as  shown  on  Figure  C.l. 

■M  i2  *-T 


Figure  C.l  Input  Places  and  Input  Transitions  of  tj 


In  the  algorithm,  we  will  denote  by: 


t.  =  {t.  ,  t.  ••••it.  } 
J  X1  x2  r 


the  corresponding  set  of  transitions.  The  backward  arrow  means  that  we 
consider  all  the  transitions  that  precede  transition  tj.  Naturally,  tj  is 
uniquely  determined  by  the  structure  of  the  Petri  Net  model  of  the  DMO. 

For  each  transition  tj,  S°  is  now  determined  from  the  schedule  of  the 
transitions  tA  that  belong  to  Vj,  by  the  relation  (see  relation  (5.7)): 


S“  -  h  .  S-p  *  Mj) 


(C.l) 


In  the  relation  (C.l),  the  maximum  is  oomputed  for  all  i  such  that  t^  s  tj. 

Let  us  recall  briefly  from  where  this  relation  comes.  Let  us  oonslder 

any  transition  t^  belonging  to  tj,  l.e.,  any  of  the  transitions 

t,  , ...,t,  ,  as  shown  on  Figure  C.l.  At  the  instant  of  the  n-th  termination 
*1  ir  - 

firing  of  transition  t^,  l.e.,  at  the  instant: 


T-S“  . 
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transition  t^  has  exactly  *produoed"  n  tokens  in  plaoe  pi}  (since  t^  has 
exactly  fired  n  times  from  the  Initial  instant).  Since  there  were 
initially  M^j  tokens  in  this  place,  the  total  number  of  tokens  available  in 
pij  during  the  interval  of  time  [0,t],  is  precisely  n  +  M* j .  Now,  each 
time  that  the  firing  of  transition  tj  was  initiated,  one  token  was  consumed 
in  Pjj.  Accordingly,  the  firing  of  transition  tj  could  not  have  been 
initiated  more  than  (n  +  M°j)  times  in  the  time  interval  [0,t],  meaning 
that: 


Sn+Mlj  >  t  .  sn  + 

Sj  "  T  Si  “i 


(C.2) 


This  intuitive  result  is  proved  in  the  previous  Appendix  (Appendix  B). 
Applying  (C.2)  for  all  the  transitions  t^  belonging  to  tj  yields: 


n+M?  .= 

s:  -  (si  *  “i* 


(C.3) 


However,  in  our  modeling,  the  firing  of  a  transition  cannot  be  initiated  if 
the  transition  is  already  executing,  meaning  that  a  transition  cannot  fire 
more  than  one  token  at  a  time  (see  Chapter  2,  Section  2.4.2).  Accordingly, 
we  have  also  the  Inequality: 


SD+1  >  SD  +  ^ 

J  3 


(C.4) 


Si 

V,y| 


A  firing  schedule  (Sj)  verifying  (C.3)  and  (C.4)  for  n  =  1,2,3,...  is  said 
to  be  feasible,  as  defined  in  Chapter  5,  Section  5.2.2.  Let  us  now  assume 
that  n,  the  number  of  repetitions,  is  large  enough,  so  that: 


n  >  1  and  n  >  M.  .  for  all  1  such  that  t.  e  t . 

ij  i  j 


Then,  (C.3)  and  (C.4)  can  be  rewritten  as: 


2# 

ft 


1 


i 


>  max  (S 


n_Mij 


°r* 


(C.5) 


It  should  be  clear  from  (C.3)  that  relation  (C.l)  determines  the  earliest 
firing  schedule,  such  that  the  firing  initiations  occur  as  soon  as  the 
transitions  are  enabled. 


In  case  where  n  is  less  than  some  M^j,  we  simply  do  not  take  into 
account  the  corresponding  value  of  i  in  the  determination  of  the  maximum: 
it  means  that  the  corresponding  transition  t^  does  not  need  to  fire  in 
order  to  enable  tj,  because  the  marking  of 
still  contains  initial  tokens).  In  particular,  (C.l)  can  be  used  to 
compute  the  instant  of  the  first  occurence  of  each  transition,  l.e.,  the 
schedule  Sj  ( j-1,2, . . . ,m) ,  if  we  adopt  the  convention  that: 


S 


J 


for  j  =  1,2, ...  ,m 


(where  m  denotes,  as  usual,  the  total  number  of  transitions).  In  that 
case,  there  are  only  two  alternatives: 


-  either  the  initial  marking  of  all  the  input  places  pAj  of 
transition  tj  is  strictly  positive,  implying  that  transition  tj 
is  enabled  at  the  Initial  instant  and  therefore  starts  firing 
at  x  =  0.  Using  relation  (C.l)  for  n  *  1,  given  the  convention, 
yields  effectively: 


V° 


or  there  is  at  least  one  place  p^j  with  a  null  marking,  i.e., 
m£j=0.  Then  transition  tj  will  therefore  be  enabled  only 
after  the  corresponding  transition  tj^  has  fired,  which  is  what 
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we  obtain,  using  (C.l): 


si  *“i 


Finally,  it  turns  out  that  relation  (C.l)  provides  a  way  to  oompute  the 
firing  schedule  (Sj)  by  iteration  on  n,  the  number  of  firing  repetitions. 
However,  in  order  to  be  able  to  oompute  Sj,  we  need  to  determine  S°,  for 
all  the  transitions  t^  that  belong  to  tj  and  for  whioh  mJj  -  0  (p^j  is  any 


one  of  the  places  that  is  not  a  resouroe  plaoe).  The  order  in  whioh  (s“) 


for  j  =  l,2,...,m  (number  of  transitions)  should  be  computed  is  therefore 
defined  by  the  partial  order  between  the  transitions,  as  determined  by  the 
slices  Tx,  Ta,...,Ta  (see  Chapter  S,  Section  5.2.3).  Let  us  desoribe  now 
the  corresponding  algorithm. 


C.2  DESCRIPTION  OF  THE  ALGORITHM 


The  algorithm  is  described  by  means  of  a  flow-chart  shown  on  Figure 
C.2.  The  analysis  of  the  Flow-chart  is  straightforward,  from  the 
explanations  given  previously.  Let  us  note  that  N  represent  the  desired 
number  of  simulations,  l.e.,  of  firing  repetitions.  It  may  also  be  seen  as 
the  total  number  of  inputs  processed  by  the  organization.  It  is  emphasized 
that  there  are  in  fact  two  steps  to  oomplete,  before  using  the  algorithm: 


(1)  determine  the  slices  Tx,  Ta,...,Ts  of  the  net,  following  the 
procedure  described  in  Chapter  5,  Section  5.1.2. 


(2)  determine  tj  for  each  transition  tj,  as  defined  previously. 


These  steps  are  easily  implemented,  using  the  Incidence  Matrix,  that 
characterizes  the  structure  of  the  Petri  Net  (see  Chapter  2,  Section 


2.1.2). 
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